Question: problem of Mfuzz
0
gravatar for gaowei1160
3.4 years ago by
gaowei11600
gaowei11600 wrote:

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 • 2.1k views
ADD COMMENTlink modified 9 months ago by luishernanhuertas0 • written 3.4 years ago by gaowei11600
0
gravatar for theobroma22
3.4 years ago by
theobroma221.1k
theobroma221.1k wrote:

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

ADD COMMENTlink modified 3.4 years ago • written 3.4 years ago by theobroma221.1k
0
gravatar for luishernanhuertas
9 months ago by
luishernanhuertas0 wrote:

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 COMMENTlink written 9 months ago by luishernanhuertas0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1744 users visited in the last hour