[EMAIL PROTECTED] wrote:
I am trying to do power calculations for the proportional odds model using the asypow library.
The code
noncenta90b10<-asypow.noncent(theta.ha=a9010,info.mat=infomatrixa90b10,constraints=constrt)
returns
Error in max(..., na.rm = na.rm) : invalid "mode" of argument.
the various arguments I've used are: a9010 [,1] [1,] -1.7357568 [2,] -0.1928619 specifying the theta.ha array as a row not a column makes no difference
infomatrixa90b10
[,1] [,2]
[1,] 0.967005807 -0.004699262
[2,] -0.004699262 0.903852346
constrt
[,1] [,2] [,3] [1,] "1" "a" "-1.92861865194525"
[2,] "1" "b" "0"
I don't think you can specify a character matrix as constraints -- for sure the package maintainer knows better.
Uwe Ligges
I'm probably missing something very simple, but I can't see what it is. Can anyone help?
_______________________________________________________________________ Most doctors use http://www.Doctors.net.uk e-mail. Move to a free professional address with spam and virus protection.
______________________________________________ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
______________________________________________ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
