I am attemting to use glmD in order to use latex(obj) afterwards. Glm works fine, glmD throws an error: > conv.go <- > glmD(cbind(go.cnt,tot.cnt-go.cnt)~sn+rcs(relAge,4)+termfac+rate:termfac+ + I(relAge*term>(term-1.25))+I((prevbal/tot.cnt)*1e-4)+prod+newmort+ + I(dnorm(relAge*term-round(relAge*term/12)*12,sd=1.2)*as.integer((relAge*term)>3))+ + I(term %in% c(36,60,72,120)) + ,data=conv, + family="quasibinomial") Error in if (!length(fname) || !any(fname == zname)) { : missing value where TRUE/FALSE needed
> conv.go <- > glm(cbind(go.cnt,tot.cnt-go.cnt)~sn+rcs(relAge,4)+termfac+rate:termfac+ + I(relAge*term>(term-1.25))+I((prevbal/tot.cnt)*1e-4)+prod+newmort+ + I(dnorm(relAge*term-round(relAge*term/12)*12,sd=1.2)*as.integer((relAge*term)>3))+ + I(term %in% c(36,60,72,120)) + ,data=conv, + family="quasibinomial") > latex(conv.go) Error in format == "" : comparison (1) is possible only for atomic and list types > str(conv) 'data.frame': 25537 obs. of 18 variables: $ prod : Factor w/ 4 levels "CD","CL","CO",..: 2 2 2 1 2 2 2 4 2 2 ... $ strt :Class 'Date' num [1:25537] 10623 11354 11382 12084 12084 ... $ matym :Class 'Date' num [1:25537] 13180 13180 13180 13180 13180 ... $ ins : Factor w/ 2 levels "0","1": 1 1 1 1 1 1 1 1 1 1 ... $ ym :Class 'Date' num [1:25537] 13149 13149 13149 13149 13149 ... $ tot.cnt: int 24 117 25 17 25 26 102 95 52 52 ... $ prevbal: num 1435380 8991665 1652058 1056949 2240357 ... $ RAM : num 139 160 150 107 160 ... $ age : num 83.4 59.4 58.9 35.3 35.3 ... $ rate : num 6.46 7.4 7.6 6 5.44 ... $ go : num 0 0 0 0 0 0 0 0 0 0 ... $ go.cnt : num 0 0 0 0 0 0 0 0 0 0 ... $ SMM.go : num 0 0 0 0 0 0 0 0 0 0 ... $ term : num 84 60 59 36 36 24 12 6 84 84 ... $ relAge : num 0.993 0.99 0.999 0.981 0.981 ... $ sn : Factor w/ 12 levels "Jan","Feb","Mar",..: 1 1 1 1 1 1 1 1 1 2 ... $ termfac: Factor w/ 7 levels "1","2","3","4",..: 6 5 5 3 3 2 1 1 6 6 ... $ newmort: int 0 0 0 0 0 0 0 0 0 0 ... Thanks everybody. Stephen B ______________________________________________ R-help@r-project.org mailing list 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.