Displaying 2 results from an estimated 2 matches for "jethi".
Did you mean:
sethi
2010 Sep 25
3
3D plot
...q_1 <- qnorm(alpha, 0, 0.05)
q_2 <- qnorm(1 - alpha, 0, 0.05)
p=matrix(0,nrow=m,ncol=N)
H=matrix(0,nrow=N,ncol=1)
p[j,i]=x[j]^2/sum(x[,i]^2)
}
H[i]=log(m)-sum(p[,i]*log(p[,i]))
}
1 - mean(q_1 <= H & H <= q_2)
}
output <- gute(a = l[,i])
return(output)
}
regards
jethi
--
View this message in context: http://r.789695.n4.nabble.com/3D-plot-tp2713818p2713818.html
Sent from the R help mailing list archive at Nabble.com.
2010 Sep 26
8
the function doesn´t work
hey, my function doesn?t work. can somebody help me?
the graphic doesn?t work and also the function. thnx a lot.
N=10
n=100
p_0=c(1/5,1-1/5)
power = function(p,m) {
set.seed(1000)
H=matrix(0,nrow=N,ncol=1)
for(i in 1:N) {
x <- matrix(rnorm(n, 0, 0.5), ncol = m)
y <- matrix(rnorm(n, 0, 0.8), ncol = m)
l <- diag(cor(x, y))
q_1 = qnorm(0.05, 0, 0.05)
q_2 = qnorm(1 - 0.05, 0, 0.05)