R help - Sep 2010 - (no subject)

hi how can i plot now this function??? have to be m= 2???  because of the
dimensions?thanks for ur help


myfun <- function(n, m, alpha = .05, seeder = 1000) {
 
  set.seed(seeder)
 
  x <- matrix(rnorm(n, 0, 0.5), ncol = m)
 
  y <- matrix(rnorm(n, 0, 0.8), ncol = m)
 
  l <- diag(cor(x, y))
 
  cat("Correlations between two random variables \n", l, fill = TRUE)
 
  gute <- function(x, m, alpha) {
 
    q_1 <- qnorm(alpha, 0, 0.05)
 
    q_2 <- qnorm(1 - alpha, 0, 0.05)
 
    p <- (x^2)/sum(x^2)
 
    H <- log(m) - sum(p * log(p), na.rm = TRUE)
 
    1 - mean(q_1 <= H & H <= q_2)
 
  }
 
  dat <- seq(0, 1, length.out = 10)
 
  output <- gute(x = dat, m = m, alpha = alpha)
 
  return(output)
 
}


 		 	   		  
	[[alternative HTML version deleted]]

Did you look at:

?curve




----------------Contact
Details:-------------------------------------------------------
Contact me: Tal.Galili@gmail.com |  972-52-7275845
Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) |
www.r-statistics.com (English)
----------------------------------------------------------------------------------------------




On Sat, Sep 25, 2010 at 2:08 PM, Kartija Oompragash
<kartija@hotmail.com>wrote:
>
> hi how can i plot now this function??? have to be m= 2???  because of the
> dimensions?thanks for ur help
>
>
> myfun <- function(n, m, alpha = .05, seeder = 1000) {
>
>  set.seed(seeder)
>
>  x <- matrix(rnorm(n, 0, 0.5), ncol = m)
>
>  y <- matrix(rnorm(n, 0, 0.8), ncol = m)
>
>  l <- diag(cor(x, y))
>
>  cat("Correlations between two random variables \n", l, fill =
TRUE)
>
>  gute <- function(x, m, alpha) {
>
>    q_1 <- qnorm(alpha, 0, 0.05)
>
>    q_2 <- qnorm(1 - alpha, 0, 0.05)
>
>    p <- (x^2)/sum(x^2)
>
>    H <- log(m) - sum(p * log(p), na.rm = TRUE)
>
>    1 - mean(q_1 <= H & H <= q_2)
>
>  }
>
>  dat <- seq(0, 1, length.out = 10)
>
>  output <- gute(x = dat, m = m, alpha = alpha)
>
>  return(output)
>
> }
>
>
>
>        [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.
>
	[[alternative HTML version deleted]]