[R] GLM HELP NEEDED!

CHIRIBOGA Xavier xavier.chiriboga at unine.ch
Wed Feb 1 14:28:18 CET 2017


Dear colleagues,


I am trying to perform a GLM. I tried again without using attach()...but still is not working.

Do you have any idea to help me?


Thank you again,


Xavier



a <- read.table(file.choose(), h<-T)
> head(a)
  time treatment transinduc
1    1   CHA0+Db     1,0768
2    1   CHA0+Db     1,0706
3    1   CHA0+Db     1,0752
4    1   CHA0+Db     1,0689
5    1   CHA0+Db     1,1829
6    1    PCL+Db     1,1423
> summary(a)
      time         treatment    transinduc
 Min.   :1.000   CHA0   :10   1,0488 : 6
 1st Qu.:1.000   CHA0+Db: 9   1,0724 : 4
 Median :1.000   Db     : 9   1,0752 : 3
 Mean   :1.433   HEALTHY:15   1,0954 : 3
 3rd Qu.:2.000   PCL    :10   1,0001 : 2
 Max.   :2.000   PCL+Db :14   1,0005 : 2
                              (Other):47
> m1<-glm(a$transinduc~a$time*a$treatment,data=a,family="poisson")
Error in if (any(y < 0)) stop("negative values not allowed for the 'Poisson' family") :
  valor ausente donde TRUE/FALSE es necesario
Además: Warning message:
In Ops.factor(y, 0) : '<' not meaningful for factors


	[[alternative HTML version deleted]]



More information about the R-help mailing list