problem of Mfuzz
3
0
Entering edit mode
7.4 years ago
gaowei1160 • 0

I'm using the software of Mfuzz. And it's a very useful method for clustering. But I got some problems when using it. I have an expression matrix stored as table in the file like this. Same format with the examples "yeast.table".

#GENE.ID    S1_1    S2_1    S3_1    S4_1    S5_1    S6_1
#TIME   1   2   3   4   5   6
#ENSSG13342 3.563908454 3.636868307 3.757429854 3.487813157 3.924786191 4.050614183
#ENSSG01130 0.816725847 0.233591149 0.045866491 0.081422805 0.078299677 0.063410441
#ENSSG18225 1.586500075 1.961881101 0.428823101 1.447323694 1.463457149 2.587627581
#ENSSG11910 0.649597066 0.882924714 0.001441974 1.22588446  1.161404661 0.550070261

But when I begin to run, it return some errors as below.

#> yeast<-table2eset(file="test2.txt")
#Read 7 items
#Read 7 items
#Read 7 items
#> yeast.r <- filter.NA(yeast, thres=0.25)
#0 genes excluded.
#> yeast.f <- fill.NA(yeast.r,mode="mean")
#> yeast.s <- standardise(yeast.f)
#> m1 <- mestimate(yeast.s)
#> m1 # xx
#[1] 1.71159
#> cl <- mfuzz(yeast.s,c=16,m=1.71)
#***Error in cmeans(exprs(eset), centers = centers, method = "cmeans", m = m, :
#  NA/NaN/Inf in foreign function call (arg 1)***
#> yeast.s
#ExpressionSet (storageMode: lockedEnvironment)
#assayData: 19189 features, 6 samples
#  element names: exprs
#protocolData: none
#phenoData
#  sampleNames: S1_1 S2_1 ... S6_1 (6 total)
#  varLabels: time
#  varMetadata: labelDescription
#featureData: none
#experimentData: use 'experimentData(object)'
#Annotation:
#> yeast
#ExpressionSet (storageMode: lockedEnvironment)
#assayData: 19189 features, 6 samples
#  element names: exprs
#protocolData: none
#phenoData
#  sampleNames: S1_1 S2_1 ... S6_1 (6 total)
#  varLabels: time
#  varMetadata: labelDescription
#featureData: none
#experimentData: use 'experimentData(object)'
#Annotation:
#> traceback()
#3: .C(R_cmeans, as.double(x), as.integer(xrows), as.integer(xcols),
#       centers = as.double(centers), as.integer(ncenters), as.double(weights),
#       as.double(m), as.integer(dist - 1), as.integer(iter.max),
#       as.double(reltol), as.integer(verbose), u = double(xrows *
#           ncenters), ermin = double(1), iter = integer(1))
#2: cmeans(exprs(eset), centers = centers, method = "cmeans", m = m,
#       ...)
#1: mfuzz(yeast.s, c = 16, m = 1.71)

I don't know what's the problem of it. Can someone help me to solve this problem? Thanks a lot.

RNA-Seq software error • 4.2k views
ADD COMMENT
0
Entering edit mode
7.4 years ago
theobroma22 ★ 1.2k

When you use the fill.NA function use the 'knn' mode rather than 'mean'.

ADD COMMENT
0
Entering edit mode
4.8 years ago

Hi, i have the same problem with the mfuzz software:

zga <-table2eset(file="rpkmglobal.txt") Read 6 items Read 6 items zga.f <- fill.NA(zga, mode = "kmn") zga.r <- filter.NA(zga.f, thres=0.25) 0 genes excluded. zga.s <-standardise(zga.r) m1 <- mestimate(zga.s) m1 [1] 1.994004 cl <- mfuzz(zga.s,c=16,m=1.99) Error in cmeans(exprs(eset), centers = centers, method = "cmeans", m = m, : NA/NaN/Inf en llamada a una funciĆ³n externa (arg 1) zga.s ExpressionSet (storageMode: lockedEnvironment) assayData: 15206 features, 5 samples element names: exprs protocolData: none phenoData sampleNames: rpf_2h rpf_5h ... rpf_48h (5 total) varLabels: time varMetadata: labelDescription featureData: none experimentData: use 'experimentData(object)' Annotation:

Could someone help me with this error?

ADD COMMENT

Login before adding your answer.

Traffic: 1693 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