adjacency matrix in WGCNA returns NA
0
0
Entering edit mode
4 months ago
Appus • 0

I have FPKM data for 6 samples (3 control, 3 treated). I tried to construct a co-expression network using WGCNA. I used the following code below however my adjacency matrix obtained from expression data in WGCNA returns NA. I don't understand why this is happening. Any suggestion will be highly appreciated. I have attached the picture of matrix below.

Thank you

library(WGCNA) library(flashClust) enableWGCNAThreads(nThreads = 64) dat=read.csv("FPKM_rep.csv") names(dat) demo=as.data.frame(t(dat[,-c(1)])) names(demo)=dat$Genes powers = c(seq(1,20,by=1)) sft = pickSoftThreshold(demo, powerVector = powers, networkType="unsigned",corFnc = cor, verbose=5) par(mfrow = c(1,2)) cex1 = 0.8 plot(sft$fitIndices[,1], -sign(sft$fitIndices[,3])sft$fitIndices[,2],xlab="Soft Threshold (power)",ylab="Scale Free Topology Model Fit,signed R^2",type="n",main = paste("Scale independence")) text(sft$fitIndices[,1], -sign(sft$fitIndices[,3])sft$fitIndices[,2],labels=powers,cex=cex1,col="red") abline(h=0.80,col="green") plot(sft$fitIndices[,1], sft$fitIndices[,5],xlab="Soft Threshold (power)",ylab="Mean Connectivity", type="n",main = paste("Mean connectivity")) text(sft$fitIndices[,1], sft$fitIndices[,5], labels=powers, cex=cex1,col="red") softPower = 8 adjacency = adjacency(demo, power=8, type="unsigned")

adjacency matrix

WGCNA • 178 views
ADD COMMENT

Login before adding your answer.

Traffic: 2693 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6