Search Archives

Guy, For a partial least squares approach look at packages plspm and pathmox. Also look at sem.additions. Regards, Mark. -- View this message in context: http://r.789695.n4.nabble.com/path-analysis-tp2528558p2530207.html Sent from the R help...

Apologies to all for the multiple posting. Don't know what caused it. Maybe it _is_ time to stop using Nabble after all... Regards, Mark. -- View this message in context: http://r.789695.n4.nabble.com/lmm-WITHOUT-random-factor...

Bob, >> Does anybody know how to eliminate the double quotes so that I can use >> the >> variable name (generated with the paste function) further in the code... ?noquote should do it. ## > "varName" [1] "varName" > noquote("varName") [1] varName Regards, Mark...

Full_Name: Mark Difford Version: R version 2.7.2 Patched (2008-08-26 r46442) OS: Windows -- Vista Submission from: (NULL) (198.54.202.146) legend(graphics) does not reset xpd to its pre-call state when it has finished...

Hi Raquel, >> routine in R to compute polychoric matrix to more than 2 categorical >> variables? I tried polycor >> package, but it seems to be suited only to 2-dimensional problems. Bur surely ?hetcor (in package polycor) does it. Regards, Mark...

Dear R-help, Has anyone implemented sparklines in the strips of a lattice plot? What I have in mind is, say, highlighting that part of a time series that one is examining in more detail in a set of lattice...

On Mar 09, 2011; 11:09am Mark Seto wrote: >> How can I extract the adjusted R^2 value from an ols object (using rms >> package)? >> library(rms) >> x <- rnorm(10) >> y <- x + rnorm(10) >> ols1 <- ols(y ~ x) ## ols1$stats...

In fact, though Chapter 12 still exists in the?TOC of the document, the relevant chapter (Statistical Models in R)?has mysteriously become chapter 11 (it was a long time ago that I last looked at it). The relevant section...

Dear List-Members, Is the application of multiple comparison procedures (using the multcomp package) to the output from a rank-based ANOVA straightforward, or do I need to take heed ? That is, is it as simple as: glht( aov(rank...

Thanks for the feedback, Uwe. Will try the most recent development version. Best, Mark Difford. ----- Mark Difford (Ph.D.) Research Associate Botany Department Nelson Mandela Metropolitan University Port Elizabeth, South Africa -- View this message in context: http://r.789695.n4...

On Mar 17, 2011; 04:29pm Thierry Onkelinx wrote: >> You cannot compare lm() with lme() because the likelihoods are not the >> same. Use gls() instead of lm() And perhaps I should have added the following: First para on page 155...

John, There is a good example of one way of doing this in "multcomp-examples.pdf" of package multcomp. See pages 8 to 10. Regards, Mark. ----- Mark Difford (Ph.D.) Research Associate Botany Department Nelson Mandela Metropolitan University Port Elizabeth...

Jane, >> Does someone know how to do fa and cfa with strong skewed data? Your best option might be to use a robustly estimated covariance matrix as input (see packages robust/robustbase). Or you could turn to packages FAiR or...

On Mar 30, 2011; 11:41am Mikhail wrote: >> I'm wondering if there's any way to do the same in R (lme can't deal >> with this, as far as I'm aware). You can do this using the...

On Jul 26, 2011 Lao Meng wrote: > glht(result, linfct = mcp(f_GROUP="Tukey") ) > Error in `[.data.frame`(mf, nhypo[checknm]) : undefined columns selected It is almost certainly the underscore in the name ("_") that is causing the problem. Try putting...

On Dec 28, 2011 at 3:47am T.M. Rajkumar wrote: > I need a way to get at the Variance Extracted information. Is there a > simple way to do the calculation. Lavaan > does not seem to output this. It does...

On Nov 04, 2011 at 6:55pm Katherine Stewart wrote: > Is there a way to determine r2 values for an SEM in the SEM package or > another way to get > these values in R? Katherine, rsquare.sem() in package sem...

Jim, >> In the glm object I can find the contrasts of the main treats vs the >> first i.e. 2v1, 3v1 and >> 4v1 ... however I would like to get the complete set including 3v2, 4v2, >> and 4v3 ... along with >> the...

Dear List, mboost:::plot.glmboost uses a function called coefpath that I wish to access (I want to change the way that the labels on axis 4 are implemented in plot.glmboost). Can anyone tell me how to "get my...

On Feb 28, 2011; 10:33pm Gregory Macfarlane wrote: >> It seems as though the mlogit.data command tries to reassign my >> row.names, >> and doesn't do it right. Is this accurate? How do I move forward? Take the time...