想利用GSE数据集做WGCNA分析,但到了设置软阈值就卡住了,报错我一直解决不了,请各路同志们支支招呀!
load("WGCNA_datExpr")
powers = c(c(1:10), seq(12,20, by = 2))
dim( datExpr )
datExpr[1:5,1:5]
sft = pickSoftThreshold(datExpr, powerVector=powers,networkType=type, verbose=0) #计算软阈值
load("WGCNA_datExpr")
powers = c(c(1:10), seq(12,20, by = 2)) #经验取值
dim( datExpr )
[1] 16 5000
datExpr[1:5,1:5]
DDX3Y EIF1AY RPS4Y1 UTY CRH
Normotensive control pregnancy_1 9.403308 7.715867 8.209864 8.529364 8.430410
Normotensive control pregnancy_2 5.633431 5.086674 5.544075 5.169188 8.138067
Normotensive control pregnancy_3 5.445700 4.937016 5.601075 5.157427 8.483299
Normotensive control pregnancy_4 9.502624 8.358172 8.468212 8.464839 7.324035
Normotensive control pregnancy_5 9.317775 7.524146 8.193160 8.289145 10.092900
sft = pickSoftThreshold(datExpr, powerVector=powers,networkType=type, verbose=0) #计算软阈值
Error in as.vector(x, "character") :
cannot coerce type 'closure' to vector of type 'character'运行结果及报错内容