Hi Jethi,
On Sat, Sep 25, 2010 at 3:42 PM, jethi <kartija at hotmail.com>
wrote:>
> hey, how can i plot this function??? thanks for ur help
I don't really have the time (or inclination, sorry) to look through
your code, but if it relates to your subject line, then there are
several ways to plot 3d objects in R.
You can look at the scatterplot3d package, or more interestingly (I
think) the rgl package, and its plot3d, lines3d, etc. functions.
-steve
>
> n=1000
> m=2
> k=n/m
> N=100
> myfun <- function(n, m, alpha = .05, seeder = 1000) {
> l=matrix(0,nrow=m,ncol=N)
> for(i in 1:N){
> set.seed(i)
> for(j in 1:m){
> x=rnorm(n,0,0.5)
> y=rnorm(n,0,0.8)
> l[j,i]=cor((x[(((j-1)*k)+1):(((j-1)*k)+k)]),
> (y[(((j-1)*k)+1):(((j-1)*k)+k)]))
> }
> }
>
> for(i in 1:N){
> for (j in 1:m){
> gute <- function() {
> ? ?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.
>
> ______________________________________________
> R-help at 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.
>
--
Steve Lianoglou
Graduate Student: Computational Systems Biology
?| Memorial Sloan-Kettering Cancer Center
?| Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact