hi, the part with the power function does not work. this is the " gute "
function in my programm. because i´m not shure on what dependet my power
function gute. and ofcourse the graphik does not work. thnx a lot:)
regards
Kaja
Date: Fri, 24 Sep 2010 19:21:17 -0700
From: ml-node+2653667-1721823656-167...@n4.nabble.com
To: kart...@hotmail.com
Subject: Re: Power Function
Hi,
What part of your code does not work or does not do what you want? I
am guessing that you want gute to return more than a single value
(which is what it does in your code). Can you show us an example of
the results you would like?
Best regards,
Josh
I reworked some of your code to simplify, clean, and wrap it all in a function:
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)
}
This results in:
> myfun(100, 5)
Correlations between two random variables
-0.5887829 0.07042411 -0.06082638
0.2395193 -0.1038213
[1] 1
On Fri, Sep 24, 2010 at 5:26 PM, jethi <[hidden email]> wrote:
>
> Hi, at first, i´m from germany, so sorry for my bad english. but i need ur
> help in R to programm a power function and to make at last a graphik of it.
>
> i have already tried my best. but it doesn´t work.the topic is: the
> homogeneity test of correlation based entropy.
>
> so it means, that i have to check if all correlations of a bivariate random
> vectors are same or not. for that i saperate the n bivariate random vectors
> (x1,y1),...,(xn,yn) in blocks m so, so that i at first calculate the
> correlation in a block. n=m*k. the numbers of the blocks m are
> user-defined. the test value is the entropy. pls help me!!!
>
>
> set.seed(1000)
> n=100
> m=5
> k=n/m
>
> x=rnorm(n,0,0.5)
> y=rnorm(n,0,0.8)
>
>
> #alpha/2 Quantil
> q_1=qnorm(0.05,0,0.05)
>
> #1-alpha/2 Quantile
> q_2=qnorm(0.95,0,0.05)
>
>
> l=matrix(0,nrow=m,ncol=1)
> for(i in 1:m){
>
>
> l[i]=print(cor((x[(((i-1)*k)+1):(((i-1)*k)+k)]),
> (y[(((i-1)*k)+1):(((i-1)*k)+k)])))
>
>
> }
> güte=function(l){
> p=matrix(0,nrow=m,ncol=1)
> for(i in 1:m){
> p[i]=l[i]^2/sum(l^2)
>
> }
> H=log(m)-sum(p*log(p))
> 1-mean(q_1<=H & H <=q_2)
> }
>
>
> l=seq(0,1,len=10)
>
>
>
>
>
> plot(l,güte, type="o",pch=20,ylim=c(0,1),col="red")
> --
> View this message in context:
> http://r.789695.n4.nabble.com/Power-Function-tp2631929p2631929.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> [hidden email] 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.
--
Joshua Wiley
Ph.D. Student, Health Psychology
University of California, Los Angeles
http://www.joshuawiley.com/
______________________________________________
[hidden email] mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
View message @
http://r.789695.n4.nabble.com/Power-Function-tp2631929p2653667.html
To unsubscribe from Power Function, click here.
--
View this message in context:
http://r.789695.n4.nabble.com/Power-Function-tp2631929p2694493.html
Sent from the R help mailing list archive at Nabble.com.
[[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.
