[R] asypow.noncent: how does it work?

[Uwe Ligges]

david_foreman at doctors.org.uk 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.
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html