Search Archives

On Nov 07, 2011 at 9:04pm Mark Difford wrote: > So here the intercept represents the estimated counts... Perhaps I should have added (though surely unnecessary in your case) that exponentiation gives the predicted/estimated counts, viz 21 (compared to...

On Oct 20, 2011; 10:07am J?rg Reuter wrote: > I want compare a Matrix row by row and at the end I want to a Matrix with > the Levenshtein-Distance. J?rg, To begin with, try the following at...

Hi Frank, >> I believe that glmnet scales variables by their standard deviations. >> This would not be appropriate for categorical predictors. That's an excellent point, which many are likely to forget (including me) since one is using a model matrix...

Nov 08, 2011; 4:58am Rolf Turner wrote: >(in response to >>> Professor Colin Aitken, >>> Professor of Forensic Statistics, >!!!) ><SNIP> >> >> Do you suppose you could provide a data-corpse for us to dissect? >Fortune nomination!!! I think Sherlock would have said...

On Aug 16, 2011 Eduardo M. A. M.Mendes wrote: > Can I gather from your email that there is nothing available on R that > deals > with dynamic models (k-step ahead and free-run)? Eduardo, There may be others, but...

Hi He Zhang, >> Is the following right for extracting the scores? >> ... >> pca$loadings >> pca$score Yes. But you should be aware that the function principal() in package psych is standardizing your data internally, which you might not want. That is...

Dear R-help, I am sorry if this is more of a stats question than an R-question, but I have found it difficult to get a clear answer by other means. Q. Would it be "wrong" to specify a...

On 2011-02-20 20:02, Karmatose wrote: >> I'm trying to include multiple variables in a non-parametric analysis >> (hah!). So far what I've managed to >> figure out is that the NPMC package from CRAN MIGHT be able...

On Jun 28, 2012 at 9:08pm Brianna Wright wrote: > I'm trying to graph some data in a boxplot-like style, but I want to set > the box and whisker limits myself... See ?bxp and ?boxplot for the format...

On Aug 27, 2011 franco salerno wrote: > ...problem with ordered probit model -polr function (library MASS). > a) how to calculate marginal effects > b) how to calculate the relative importance of each independent variables Hi Franco, Have a look at John...

Hi Angelo, Look carefully at package vcd; and at log-linear models (+ glm(..., family=poisson)). For overdispersion there are more advanced methods. HTH, Mark. Angelo Scozzarella wrote: > > Hi, > > how can I treat data organised in classes and frequencies? > > Ex. > > class...

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() Hi Thierry, Of course, I stand subject to correction, but unless something dramatic...

On Feb 27, 2012 at 9:30pm Joyous Fisher wrote: > Q: is there a way to do princomp or another method where every row has at > least one missing column? You have several options. Try function nipals in packages ade4...

Marcel, Here is one way: spplot(meuse.grid, zcol = "part.a", par.settings = list(panel.background=list(col="grey"))) ## trellis.par.get() trellis.par.get()$panel.background Regards, Mark. > On 03/05/2011 01:06 PM, Marcel J. wrote: >> Hi...

Dear List, I would greatly appreciate help on the following matter: The RuleFit program of Professor Friedman uses partial dependence plots to explore the effect of an explanatory variable on the response variable, after accounting for the average effects of...

Hi Pavel, First, annonations should have the same cex-size on each axis. That said, the way that this is implemented is not too cexy (ouch!). You need to plot without axes, e.g. plot(obj, axes=F), then you...

Hi Emma, >> R gives you the tools to work this out. ## Example set.seed(7) TDat <- data.frame(response = c(rnorm(100, 5, 2), rnorm(100, 20, 2))) TDat$group <- gl(2, 100, labels=c("A","B")) with(TDat, boxplot...

On May 04, 2011 at 8:26pm Lucia Ca?as wrote: > I would like to combine lattice plot (xyplot) and standard R plot (plot > and plotCI) in an unique figure. Hi Lucia, Combining the two systems can be done. See...

Hi J?rg, >> I haven't found anything in par()... No? Well don't bet your bottom $ on it (almost never in R). ?par (sub mgp). Mark. J?rg Gro? wrote: > > Hi, > > How can I make the distance between an...

>> I'd prefer to stick with JPEG, TIFF, PNG, or the like. I'm not sure EPS would fly. Preferring to stick with bitmap formats (like JPEG, TIFF, PNG) is likely to give you the jagged lines and other distortions...