Displaying 1 result from an estimated 1 matches for "siglevel".
Did you mean:
loglevel
2011 May 04
1
hurdle, simulated power
...n
#
#example power analysis
#return three power, power1 for zero part, power2 for negative binomial part
#power3 for joint test,significance level can be set, default is 0.05
#M is simulation time
#require pscl package
#library(pscl)
mypower <- function(n, beta0, beta1, alpha0, alpha1, theta,
siglevel=0.05, M=1000){
myfun <- function(n,beta0,beta1,alpha0,alpha1,theta,siglevel){
data <- mysim(n,beta0,beta1,alpha0,alpha1,theta)
require(pscl, quietly = TRUE)
res <- hurdle(y ~ trt, data = data, dist = "negbin", trace = FALSE)
est <- coef(res)#[c(2,4)]
#v<-res$vcov...