Hi Jethi, On Sat, Sep 25, 2010 at 3:42 PM, jethi <kart...@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@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 ______________________________________________ R-help@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.