similar to: Additional output in cancor

Dear R's! I am strugling with cancor procedure in R. I cannot figure out the meaning of xcoef and of yxcoef. Are these: 1. standardized coefficients 2. structural coefficients 3. something else? I have tried to simulate canonical correlation analysis by checking the eigenstructure of the expression: Sigma_xx %*% Sigma_xy %*% Sigma_yy %*% t(Sigma_xy). The resulting eigenvalues were the same

Hi, from what I understand about the canonical correlation function 'cancor', it looks for correlations in two sets of variables, each represented in matrix form. Right? Sounds exactly like what I need. I have tried the following but I am not sure how to interpret the results. AudioPCs <- c(ArTHarF0PCA$x[,2], ArTHarF1PCA$x[,2], ArTHarF2PCA$x[,2], ArTHarF3PCA$x[,2],

Hello, I'm a beginning user of R, now I have a question about canonical correlation analysis (cca). In R,there is a function "cancor.R" used for cca; For example X(n*p1) and Y(n*p2)are the two matrix to be analyzed. In the example given by R, when n> max(p1,n2), cancor(X,Y) works; but when n<p1 or n<p2, cancor(X,Y) doesn't work well because cancor$cor == 1; how to cope

Dear all, I didn't find any graphical solution in the package "mva" to plot the canonical scores from a CCA (canonical correlation analysis). Does anybody knows how to plot or has anybody already programmed : - the map of the canonical scores, - the graph of the canonical weights, - the correlation circle i.e. the canonical loadings ? Thank you for help ...

When I use the function cancor of mva package, I found that it doesn't work when one of the matrix has only one column, or both have only one column. The function in Splus 5 with the same name works under those situations. The version of R I am using is 0.63.2 (released on Jan., 1999) on Solaris. Kenny Ye Assistant Professor Department of Applied Math and Statistics SUNY at Stony Brook

I'm trying to write a formula method for canonical correlation analysis, that could be called similarly to lm() for a multivariate response: cancor(cbind(y1,y2,y3) ~ x1+x2+x3+x4, data=, ...) or perhaps more naturally, cancor(cbind(y1,y2,y3) ~ cbind(x1,x2,x3,x4), data=, ...) I've adapted the code from lm() to my case, but in this situation, it doesn't make sense to include an

Dear R-users When running a glm polynomial model with one explanatory variable (example Y~X+X^2), with a poisson or binomial error distribution, the predicted values obtained from using the predict() function and those obtained from using the coefficients from the summary table "as is" in an equation of the form Y=INTERCEPT+ XCoef x X + XCoef x X^2, differ considerably. The former are

good morning to the R-help list!, I would like to ask if there is an alternative to the S ploting parameter "cxy" in R. I am trying to port some S code to R that makes use of this "cxy" parameter to scale the graphics, and I don't find this parameter in help(par). thank you, IOsu ------------------------------------------------------------------------------- Jes?s

Hello, I want to implement LLVM backend for a specific VLIW hardware. I am working on defining its instruction set, and assembly language. The hardware has two pipelines, int and float. Each pipeline can do 3 operations/cycle, 3 operations forms an instruction. One of the Integer Instruction looks like this: add Ri, Rj, Rk; add Rl, Rm, Rn; add Ro, Rp, Rq An int instruction and a float

I would like to draw only the ticks of an axis, but not the axis itself. I don't think this can be done using axis(), so I am trying to write a cut-down version in R, which only draws ticks. The point at which I am stuck is that the length of a tick is set by par("tcl") as a fraction of the "height of a line of text". So I would like to draw a line whose length is also

kalish at psy.uwa.edu.au wrote: > > I'm a newbie at R, and can't get libraries to really work. > I did this: > > library(help = mva) > cancor Canonical Correlations > cmdscale Classical (Metric) Multidimensional Scaling > dist Distance Matrix Computation > hclust Hierarchical Clustering

Is there an equivalent to the S-PLUS par()$cxy in R? Also, is it possible to get notice taken of some equivalent of \n in R when using text()? For what it is worth I am using RW-0.64.1. John Maindonald email : john.maindonald at anu.edu.au Statistical Consulting Unit, phone : (6249)3998 c/o CMA, SMS, fax : (6249)5549 John Dedman Mathematical

Dear all, Does anybody know the probable reason why <= gives false when it should give true? These two variables are of the same type, and everything works in the cycle but then it stops when they are equal. this is the output result > Rk[47] <= RB[21] [1] FALSE > Rk[47] [1] 0.002842007 > RB[21] [1] 0.002842007 Thanks a lot. [[alternative HTML version deleted]]

Hi Violeta, I compiled with clang -Oz and clang -Oz -mno-red-zone for comparisons against Oz and clang -O0 and clang -O0 -mno-red-zone for comparisons against a default clang. I unfortunately don’t have the clang commit I worked with on my home laptop, and don’t have access to the computer I was using at Apple, so I can’t help you there. Jessica > On Sep 28, 2016, at 11:22 AM, Violeta

I often use something like old.par <- par(set someting) on.exit(par(old.par)) but in R 1.7.0. I now get warnings: > old.par <- par() > par(old.par) Warning messages: 1: parameter "cin" can't be set in: par(args) 2: parameter "cra" can't be set in: par(args) 3: parameter "csi" can't be set in: par(args) 4:

Having searched and searched I still haven't found what's the problem with my data (I've attached the relevant file). Every time I tried to use the CANCOR-Function I got error messages. So I turned to check my two sets of variables separately by using the "colMeans"-Function. With one set no problem .. but with the other. My input : > struktur <-

Dear r-help list, I'm trying to use multiple imputation for my MSc thesis. Having good exemples using the Hmisc package, I tried the aregImpute function. But with my own dataset, I have the following error : Erreur dans matxv(X, xcof) : columns in a (51) must be <= length of b (50) De plus : Warning message: In f$xcoef[, 1] * f$xcenter :   la taille d'un objet plus long n'est pas

I've been using par() to check the graphics parameters associated with both plot(<fitted linear model>) and plot(<grouped data object>). AFAIK the only differences are in the $cxy, $usr, $xaxp, and $yaxp parameters but the background color for the grouped data plot is grey while the linear model plot has a white background. When I've tried par(bg = "white") prior to

Hi I'm using boot package for some analysis on linear regression coeficients. My problem is that I can not compute bca intervals, I get an error message > bca.ci(blm8901,index=1) Error in if (!all(rk > 1 & rk < R)) warning("Extreme Order Statistics used as Endpoints") : missing value where logical needed The problem is the empinf.reg function that is

dear R-help, i have looked carefully through the R-help archives for information on how to suppress whiskers in a bwplot. someone asked this question a while ago, but the answer he received is not available in the archives. but i did manage to get my hands on a panel function (called "my.panel") that is supposed to do this (the function is reproduced at the end of the email, below).