search for: rxz

Displaying 4 results from an estimated 4 matches for "rxz".

Did you mean: rx
2009 May 11
0
Partial correlation function required
...omit(data) } - Ignored: # recursive formula if(dim(z)[2] == 1){ tdata <- na.omit(data.frame(data[,1],data[,2])) rxy <- cor(tdata[,1],tdata[,2],m=method) tdata <- na.omit(data.frame(data[,1],data[,-c(1,2)])) rxz <- cor(tdata[,1],tdata[,2],m=method) tdata <- na.omit(data.frame(data[,2],data[,-c(1,2)])) ryz <- cor(tdata[,1],tdata[,2],m=method) rxy.z <- (rxy - rxz*ryz)/( sqrt(1-rxz^2)*sqrt(1-ryz^2) ) return(rxy.z) }else{...
2009 Jul 13
0
testing equality of two dependent correlations + normality issue
...stion. I would like to test equality of two correlation coefficients in a setting with three variables X,Y,Z, i.e. equality of r(X,Y) and r(Z,Y). I have found a formula to transform the "2 dependent correlations difference" to t-distribution with N-3 df: t = (rxy - rzy)* SQRT[{(n - 3)(1 + rxz)}/ {2(1 - rxy^2 - rxz^2 - rzy^2 + 2rxy*rxz*rzy)}] (Blalock, H., 1972. Social Statistics. NY: McGraw-Hill. Page 406-7). Am actually not sure whether this is exact or approximate (even given normality assumption, the Fisher's Z-transform which this is - I assume - based on, is approximate, right?...
2011 Jul 20
0
The C function getQ0 returns a non-positive covariance matrix and causes errors in arima()
...- ttheta[C[B]] SX <- A2%*%SU%*%t(A2) ## Now, compute correlation matrix between X and Z C1 <- matrix(0,q,q) F <- row(C1)-col(C1)+1 E <- (F>=1) & (F<=p+1) C1[E] <- tphi[F[E]] g <- matrix(0,q,1) if (q) g[1:q,1] <- ttheta[1:q] rXZ <- forwardsolve(C1,g) SXZ <- matrix(0, p, q+1) F <- col(SXZ)-row(SXZ) E <- F>=1 SXZ[E] <- rXZ[F[E]] S <- rbind(cbind(SX,SXZ),cbind(t(SXZ),diag(q+1))) } Q0 <- A%*%S%*%t(A) } The second way is to resolve brutally the equation of Gardner et al....
2009 Jul 23
1
[PATCH server] changes required for fedora rawhide inclusion.
...X}99?pyP$aW4)F{F7C6(z%yyzJ6r{)Ab? zh#Kn at d$~(1(jq1`(VFcdTb<0WfB=i+s>KW!$8|EC>Pe_4yRosSMiTna`oB53UY6iV zHF_{|SH at 9O!ttfu at dt~X=z*h4JUKh+X~Z3NxUWWdAr;&o?kB8ChCRC8{u;|J!WSq9 zryQ=1Z31muB$FY%%p3Lv(7}Geu=Kbnn8J=~?+97MQSGmhk at RTQR!&A7!E#)x?N&1* z!=alHe#{xkuw9mD*V5xUD{C~W?V8g3HRXZ23hb`Ey!FAlh!ZB`4%3}VG2Ik%_pL~e zt5w!<>Bas+($0~cFma|c7B2WzH^XAcP1rBHME6VKn|G~``0nrqaH`Ot(V2iGg>;3h zg2OH6Fss5k9O5uoRAcV$ss5lm%Gu6=bWAvQ1Yrr%5XR?bBArw}45v6+`9g6AOu?39 z1>DvZf;?Ad6i3Hk^(+M?gn-Q3fE1jMVkvjLuhcBQo}@Oigdou~?-0yB6t<?)Kg4a; z(M=`_8e!-X&(QMXsfeGxtUz!$tf4T3C1T6-)s*|P>S}9...