Displaying 5 results from an estimated 5 matches for "rinv".
Did you mean:
ring
2009 May 11
0
Partial correlation function required
...c(x,z0,zc,method=method,na.rm=na.rm)
ryz0.zc <- pcor.rec(y,z0,zc,method=method,na.rm=na.rm)
rxy.z <- (rxy.zc - rxz0.zc*ryz0.zc)/(
sqrt(1-rxz0.zc^2)*sqrt(1-ryz0.zc^2) )
return(rxy.z)
}
}
and this piece of code also i tried to use.
Rinv <- solve <http://wiki.r-project.org/rwiki/rhelp.php?id=solve>(R)
D <- diag <http://wiki.r-project.org/rwiki/rhelp.php?id=diag>(1 / sqrt
<http://wiki.r-project.org/rwiki/rhelp.php?id=sqrt>(diag
<http://wiki.r-project.org/rwiki/rhelp.php?id=diag>(Rinv)))
P &l...
2000 Feb 25
2
partial correlation coefficients in R?
Hello,
after thorough searching of the R help files as well as S+-help, I'm coming
to the list: Is there a possibility to compute partial correlation
coefficients between multiple variables (correlation between two paired
samples with the "effects of all other variables partialled out")? All I
seem to find are the standard Pearson correlation coefficients (with cor())
and no clue
2007 Feb 22
0
Error in solve.default
...error in solve.default. The function was
originally written for an older version of SPlus. Can anyone give me some
insights into where the problem is?
Thanks
R 2.4.1 on MAC OSX 2mb ram
Mark Grant
markg at uic.edu
> attach(Aspirin.frame)
> hblm(Diff ~ 1, s = SE)
Error in solve.default(R, rinv) : 'a' is 0-diml
> traceback()
6: .Call("La_dgesv", a, b, tol, PACKAGE = "base")
5: solve.default(R, rinv)
4: solve(R, rinv)
3: summary.blm(fit)
2: eb.calc(rho[i], X, Y, s.e., df.se, corrs, prior, ...)
1: hblm(Diff ~ 1, s = SE)
>
> hblm
function(formula, s.e.,...
2001 Nov 20
0
Summary: non-negative least squares
...g")
m <- NCOL(x)
if (length(eps)==1) eps <- rep(eps,m)
x <- x * wsqrt
y <- y * wsqrt
#sometimes a rescaling of x and y helps (if solve.QP.compact fails
otherwise)
xscale <- apply(abs(x),2,mean)
yscale <- mean(abs(y))
x <- t(t(x)/xscale)
y <- y/yscale
Rinv <- backsolve(qr.R(qr(x)),diag(m))
cf <- solve.QP.compact(Dmat=Rinv,dvec=t(x)%*%y,Amat=rbind(rep(1,m)),
Aind=rbind(rep(1,m),1:m),bvec=eps*xscale/yscale,
factorized=TRUE)$sol
cf <- cf*yscale/xscale #scale back
cf
}
##- #PS: Why the scaling...
2009 Jul 23
1
[PATCH server] changes required for fedora rawhide inclusion.
...|9)JhZWq?YKbKPs+;d7)+h5ZdZ5#Cvr@<Oc-
z(uaKo80~w0T~5m0TumgvkbOY-&PiiBM4?~eP?h-dbuKN3ZUVEDFvQnM+{6^CJG_0t
zm1s;yqc?3g)N)TUP(W)gjH&-QmDd9xhx)f9{$VB2zbNV*jsW8viqq4er*=ECY?>%%
z6?ehhpGNODVwO8!XIM2K8rUjgUJMv14JVk~RMZ+j)mBP)!Ys(#G|1e*I}*fLrxTjS
zI$!7tq^(KdrBi^3j}>4Gh=VBAI#?AeEeD7LDFT81dsE(RInVs`ytdZBWRKoJzf7s|
zNm|oi+NYA~IX^GH#b?}V2m>%A8N?k9Gp6Epx==GzUhb6t!?5lI38J>|HUC?R{*Y4C
z|0EMJL7RL%Jw&;w(dFa3>Fvw>g`?xklBR&pbNKWafsRgfp-`+6$wN<N6c;DuU_27C
zRBD|d$VH8xei%4sS4APKe-`5UY$zTXp*C(Lsi8Zm=VcD}ShlWwshzfgmvFY5FOqh?
z at 2O`G$5{*+CZCUYHGhy~IjMoOmA7au4UouNzL=F_MG#h6N~)Msw1_|...