Thanks for all your help. I will try to fix this.
Very helpful.
Best,
Daofeng
On Fri, Jun 7, 2013 at 11:25 AM, Sarah Goslee wrote:
> As Marc already pointed out, take a close look at this part of your loop:
>
> R> i <- 6
> R>
> R> y <- as.numeric(data[i,-1])
> R> y
> [1] 3 3 3 3 4 4 4 4
> R> g
Thanks Marc...still didn't get it...
Tried this...
> str(group)
num [1:8] 1 1 1 1 0 0 0 0
> fit=glm.nb(y~group)
Error in while ((it <- it + 1) < limit && abs(del) > eps) { :
missing value where TRUE/FALSE needed
> group <- c('1','1','1','1','0','0','0','0')
> str(group)
chr [1:8] "1" "1" "1" "
Thank you Sarah and Marc for your fast and nice response.
Apology for didn't include all information.
I have a input file like following:
gene1 18 15 13 13 16 9 20 24
gene2 15 8 8 7 0 12 18 4
gene3 10 9 8 12 9 11 12 12
gene4 4 0 4 3 0 5 0 0
gene5 0 1 0 0 1 5 1 0
gene6 3 3 3 3 4 4 4 4
gene7 0 4 0
Sorry Sarah.
> dput(dat)
structure(list(gene = structure(c(1L, 3L, 4L, 5L, 6L, 7L, 8L,
9L, 10L, 2L), .Label = c("gene1", "gene10", "gene2", "gene3",
"gene4", "gene5", "gene6", "gene7", "gene8", "gene9"), class = "factor"),
b1 = c(18L, 15L, 10L, 4L, 0L, 3L, 0L, 4L, 11L, 6L), b2 = c(15L,
8L,
As Marc already pointed out, take a close look at this part of your loop:
R> i <- 6
R>
R> y <- as.numeric(data[i,-1])
R> y
[1] 3 3 3 3 4 4 4 4
R> group
[1] 1 1 1 1 0 0 0 0
R> fit <- glm.nb(y~group)
Error in while ((it <- it + 1) < limit && abs(del) > eps) { :
missing value where TRUE/FALSE neede
Hi,
On Fri, Jun 7, 2013 at 11:36 AM, Daofeng Li wrote:
> Thank you Sarah and Marc for your fast and nice response.
> Apology for didn't include all information.
>
> I have a input file like following:
>
> gene1 18 15 13 13 16 9 20 24
> gene2 15 8 8 7 0 12 18 4
> gene3 10 9 8 12 9 11 12 12
> gene4
On Jun 7, 2013, at 10:36 AM, Daofeng Li wrote:
> Thank you Sarah and Marc for your fast and nice response.
> Apology for didn't include all information.
>
> I have a input file like following:
>
> gene1 18 15 13 13 16 9 20 24
> gene2 15 8 8
On Jun 7, 2013, at 9:44 AM, Daofeng Li wrote:
> Dear R Community,
>
> I have encountered a problem while using the R function glm.nb.
> The code that produce the error was following two lines:
>
> group=c(1,1,1,1,0,0,0,0)
> fit=glm.nb(y~group)
>
> While the y contains 8 sets of number like:
>
Hi,
On Fri, Jun 7, 2013 at 10:44 AM, Daofeng Li wrote:
> Dear R Community,
>
> I have encountered a problem while using the R function glm.nb.
> The code that produce the error was following two lines:
>
> group=c(1,1,1,1,0,0,0,0)
> fit=glm.nb(y~group)
>
> While the y contains 8 sets of number li
Dear R Community,
I have encountered a problem while using the R function glm.nb.
The code that produce the error was following two lines:
group=c(1,1,1,1,0,0,0,0)
fit=glm.nb(y~group)
While the y contains 8 sets of number like:
gene2750 1 0 0 1 5 1
10 matches
Mail list logo
