HELP with GLM
Dear Xavier: It sounds like you have a right mess! Perhaps others cleverer and more diligent than I can sort through it and diagnose the problem. However, it really *does* sound like you need to step back, take a deep breath, and spend some time with an R tutorial or two (there are many good ones on teh web) to learn how R works before proceeding further. You seem to be groping around in the dark. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Jan 31, 2017 at 6:35 PM, CHIRIBOGA Xavier
<xavier.chiriboga at unine.ch> wrote:
Dear colleagues, I am trying to perform GLM ..but I got some objects masked: and an error message below 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
attach(a)
The following objects are masked from a (pos = 12):
time, transinduc, treatment
The following objects are masked from a (pos = 13):
time, treatment
The following objects are masked from a (pos = 14):
time, treatment
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(transinduc~time*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
I DO NOT HAVE NEGATIVE VALUES IN MY DATASET, Do you know what is going wrong?
Thank you for you help,
Xavier
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.