Displaying 1 result from an estimated 1 matches for "icg_rcpivot_group2".
2007 Jul 16
2
Error while fitting Partial Proportional Odds model using vglm
...et an
error while executing the code. Here is the data, code, and error:
Data: first row is the column names. a = age, and 1,2,3, 4 and 5 are
condition grades.
a 1 2 3 4 5
1 1 0 0 0 0
2 84 2 7 10 2
3 16 0 6 6 2
4 13 0 3 4 0
5 0 0 0 1 0
Library(VGAM)
rc13<-read.table("icg_rcPivot_group2.txt",header=F)
names(rc13)<-c("a","1","2","3","4","5")
ppo<-vglm(cbind(rc13[,2],rc13[,3],rc13[,4],rc13[,5],rc13[,6])~a,family =
cumulative(link = logit, parallel = F , reverse = F),na.action=na.pass,
data=rc13)
summary(ppo)
I...