Error in ORDGLM function
I am trying this simple model.
Following is the code:
newrate <- c(0,1,2,0,1,2,0,1,2,0,1,2)
totalast <- c(100,200,300,120,215,280,75,145,255,125,275,310)
ordglm(newrate~totalast,link="logit",maxiter=100, weights=1)
It runs for few seconds and get the error back
Error in ordglm(newrate ~ totalast, link = "logit", maxiter = 100, weights = 1) :
Object "PearsRes" not found
That error is not from my code: I do not have a PearsRes. Your three lines work for me:
library(gnlm)
Loading required package: rmutil
newrate <- c(0,1,2,0,1,2,0,1,2,0,1,2) totalast <- c(100,200,300,120,215,280,75,145,255,125,275,310) ordglm(newrate~totalast,link="logit",maxiter=100, weights=1)
Call:
ordglm(newrate ~ totalast, link = "logit", maxiter = 100, weights = 1)
-Log likelihood 149.5272
AIC 152.5272
Iterations 5
Location coefficients
estimate s.e.
(Intercept) -51.3812 0.05031
totalast 0.1969 3.41556
Intercept contrasts
estimate s.e.
b[2] -5.513 11.46
Correlation matrix
1 2 3
1 1.0000 -0.9890 0.7268
2 -0.9890 1.0000 -0.8065
3 0.7268 -0.8065 1.0000
Jim
I tried to use the example in the package # McCullagh (1980) JRSS B42, 109-142 # tonsil size: 2x3 contingency table y <- c(0:2,0:2) carrier <- gl(2,3,6) wt <- c(19,29,24,497,560,269) ordglm(y~carrier, weights=wt) It works fine. So everything is installed correctly. may be it has something to do with the weight? Any help is welcome. Thanks a lot. Kamal Desai -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._