# -----
Com.ave.gd<-read.csv("Com.SCA.gd.csv")
hete<-as.matrix(Com.ave.gd[,-c(1:3)])
a<-rcorr(hete, type="pearson")
write.csv(a$r,"corr_SCA_gd.csv")
write.csv(a$P,"corr_SCA_gd_pvalue.csv")
# ----匹配heterosis SCA 和gd
gd<-read.csv("gd.nei between parents.csv")
gd$com1<-paste0(gd$p1,gd$p2)
gd$com2<-paste0(gd$p2,gd$p1)
# 1 heterosis 和gd
hete<-read.csv("heterosis_growth.csv")
names(hete)
hete<-hete[,-11]
for(i in 1:dim(hete)[1]){
for(j in 1:dim(gd)[1]){
if(hete$X[i]==gd$com1[j])(hete$GD[i]<-gd$GD[j])
}
}
for(i in 1:dim(hete)[1]){
for(j in 1:dim(gd)[1]){
if(hete$X[i]==gd$com2[j])(hete$GD[i]<-gd$GD[j])
}
}
# 2 SCA 和gd
sca<-read.csv("Com.SCA.gd.csv")
sca$com2<-paste0(sca$Mum,sca$Dad)
for(i in 1:dim(sca)[1]){
for(j in 1:dim(gd)[1]){
if(sca$com2[i]==gd$com1[j])(sca$GD[i]<-gd$GD[j])
}
}
for(i in 1:dim(sca)[1]){
for(j in 1:dim(gd)[1]){
if(sca$com2[i]==gd$com2[j])(sca$GD[i]<-gd$GD[j])
}
}
# -----cor
library(Hmisc) #相關(guān)顯著性檢驗
# 1 heterosis 和gd
hete.cor<-as.matrix(hete[,-1])
a<-rcorr(hete.cor, type="pearson")
# a<-cor(hete,method = c("pearson"), use = "complete.obs")
write.csv(a$r,"corr_hete_gd.csv")
write.csv(a$P,"corr_hete_gd_pvalue.csv")
# 1 sca 和gd
names(sca)
sca.cor<-as.matrix(sca[,-c(1:3,13)])
b<-rcorr(sca.cor, type="pearson")
# a<-cor(hete,method = c("pearson"), use = "complete.obs")
write.csv(b$r,"corr_sca_gd.csv")
write.csv(b$P,"corr_sca_gd_pvalue.csv")
heterosis SCA 和gd 關(guān)系
最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
- 文/潘曉璐 我一進店門锌仅，熙熙樓的掌柜王于貴愁眉苦臉地迎上來，“玉大人墙贱，你說我怎么就攤上這事热芹。” “怎么了惨撇？”我有些...
- 文/不壞的土叔 我叫張陵伊脓，是天一觀的道長。 經(jīng)常有香客問我魁衙，道長报腔，這世上最難降的妖魔是什么株搔？ 我笑而不...
- 正文 為了忘掉前任，我火速辦了婚禮纯蛾，結(jié)果婚禮上邪狞，老公的妹妹穿的比我還像新娘。我一直安慰自己茅撞，他們只是感情好帆卓，可當我...
- 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著米丘，像睡著了一般剑令。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上拄查，一...
- 文/蒼蘭香墨 我猛地睜開眼，長吁一口氣：“原來是場噩夢啊……” “哼糊探！你這毒婦竟也來了钾埂？” 一聲冷哼從身側(cè)響起，我...
- 正文 年R本政府宣布奶是，位于F島的核電站，受9級特大地震影響，放射性物質(zhì)發(fā)生泄漏聂沙。R本人自食惡果不足惜秆麸，卻給世界環(huán)境...
- 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望及汉。 院中可真熱鬧沮趣，春花似錦、人聲如沸坷随。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽温眉。三九已至缸匪，卻和暖如春，著一層夾襖步出監(jiān)牢的瞬間类溢，已是汗流浹背凌蔬。 一陣腳步聲響...
推薦閱讀更多精彩內(nèi)容
- 小學(xué)語文修改病句的方法 修改病句是小學(xué)語文考試中常見的題型殖妇，在修改病句之前，我們應(yīng)該清晰的了解有哪些病句現(xiàn)象破花，下面...
