search for: mark_difford

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 plots. Regads,. Mark Difford. PS: (Andreas Loffler has implemented a simple but functional version for TeX/LaTeX: http://www.tug.org/tex-archive/help/Catalogue/entries/sparklines.html)

Dear List, I would appreciate help on the following matter: I am aware that higher dimensional contingency tables can be analysed using either log-linear models or as a poisson regression using a generalized linear model: log-linear: loglm(~Age+Site, data=xtabs(~Age+Site, data=SSites.Rev, drop.unused.levels=T)) GLM: glm.table <- as.data.frame(xtabs(~Age+Site, data=SSites.Rev,

...softwares don't really care about who gets credit ( Who wrote Windows Vista ,or Mac OS or Ubuntu Linux), and the NYT is a newspaper not a journal. Does any student, or teacher for that matter care whether Newton or Leibntiz invented calculas. On Thu, Feb 5, 2009 at 11:35 AM, Mark Difford <mark_difford@yahoo.co.uk>wrote: > > >> I think that all appeared on January 8 in Vance's blog posting, with a > >> comment on it by David M Smith on Jan 9. So those people have -27 days > > Then there was no need for vituperative comments (not from you, of course): > simpl...

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 the other variables. The plot method [plot(summary(rq(y ~ x1 + x2, t=seq(.1,.9,.05))))] of Professor Koenker's quantreg program

Hi: I have a large dataset mydata, of 1000 rows and 1000 columns. The rows have gene names and columns have condition names (cond1, cond2, cond3, etc). mydata<- read.table(file="c:/file1.mtx", header=TRUE, sep="") I applied PCA as follows: data_after_pca<- prcomp(mydata, retx=TRUE, center=TRUE, scale.=TRUE); Now i get 1000 PCs and i choose first three PCs and make a

He everybody, Well I have a biplot CCA-like origined from plot.cca (vegan package). I need to rotate on y axis the lines and symbols of constrained and sites representation. If I do that on an image editor, I rotate everything, including titles, axes labels and positions. I just need to rotate the inner par and keep the variables names (constrained) and symbols in the new positions but with the

Hi, I'm a beginner with R, but I'm getting excellent results with it. Well I've got an capscale object (vegan package) and I want to made a biplot with symbols representing six groups of areas. With the plot.cca function and some par attributes (like 'labels') I was able to substitute the samples names for keyboard symbols ('x','o','#',ect), but it's

Dear List-members, Hopefully someone will help through my confusion: In order to get the same coefficients as we get from the following ## require (MASS) summary ( lm(Gas ~ Insul/Temp - 1, data = whiteside) ) ...................... we need to do the following (if we use model.matrix to specify the model) ## summary ( lm(Gas ~ model.matrix(~ Insul/Temp - 1) - 1, data = whiteside) )

...___________________________________________ AOL Email goes Mobile! You can now read your AOL Emails whilst on the move. Sign up for a free AOL Email account with unlimited storage today. ------------------------------ Message: 5 Date: Sun, 3 Aug 2008 06:55:41 -0700 (PDT) From: Mark Difford <mark_difford at yahoo.co.uk> Subject: Re: [R] How do I get an inset graph (i.e. graph within a graph)? To: r-help at r-project.org Message-ID: <18798809.post at talk.nabble.com> Content-Type: text/plain; charset=us-ascii Hi Arthur, This can be done quite easily using the appropriate arguments liste...

...and "Reference Cards". Further, on the "Getting Help" tab, there are also links to the ESS section of the R- Wiki and the Emacs - Wiki. Martin Maechler, ETH Zurich ------------------------------ Message: 77 Date: Mon, 31 Aug 2009 02:30:31 -0700 (PDT) From: Mark Difford <mark_difford@yahoo.co.uk> Subject: Re: [R] test for bimodality&In-Reply-To= To: r-help@r-project.org Message-ID: <25220627.post@talk.nabble.com> Content-Type: text/plain; charset=us-ascii Hi John, >> Has a test for bimodality been implemented in R? You may find the code at the URL below u...

...tures */*Trafic -- Socio-économie */2, rue de Bruxelles, BP 275 59019 Lille cedex FRANCE Tél: +33 (0) 3 20 49 60 70 Fax: +33 (0) 3 20 49 63 69 [[alternative HTML version deleted]] ------------------------------ Message: 53 Date: Sun, 21 Sep 2008 23:48:47 -0700 (PDT) From: Mark Difford <mark_difford@yahoo.co.uk> Subject: Re: [R] Variable Selection for data reduction and discriminant anlaysis To: r-help@r-project.org Message-ID: <19602702.post@talk.nabble.com> Content-Type: text/plain; charset=us-ascii Hi Gareth, >> My data is transformed to the clr or alr under Aitchison...

..._____ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > --Forwarded Message Attachment-- From: mark_difford at yahoo.co.uk To: r-help at r-project.org Date: Mon, 23 Aug 2010 13:08:01 -0700 Subject: Re: [R] Coinertia randtest Hi Petar, >> I dunno why, but I cannot make randtes[t].coinertia() from ade4 package >> working. I have two nice distance matrices (Euclidean): >> Could anyone...

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 nested model and retain a common intercept, e.g. lm(NH4 ~ Site/TideCode + 1) I am aware (?) that my Site-coefficients are now calculated relative to my reference Site (treatment.contrasts), *but* that

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(NH4) ~ Site, data=mydat), linfct=mcp(Site="Tukey") ) Thanks in advance for your help. Regards, Mark Difford. -- View this message in context:

On May 24, 2011; 11:06pm Belle wrote: > Does anyone know how to run Multinomial logistical Model in R in order to > get predicted probability? Yes. I could stop there but you shouldn't. The author of the package provides plenty of examples (and two good vignettes) showing you how to do this. Suggest you do some work in that area. Look especially at how model formulas are

On Nov 12, 2011 at 8:29pm Alex wrote: > Has anyone used SIAR package add on? I posted a reply to an earlier question from you on this subject. See http://r.789695.n4.nabble.com/Errors-in-SIAR-td4029804.html. In it I note that there are problems with the function from siar (not SIAR) you are using, but that this may not be "your" problem, that the function calls for matrices (you

Dear R-help, I posted this on 4 Oct but got no response (I wasn't even told to go away and do some more background reading ;) ). I am reposting it in the, perhaps, vain hope that someone with knowledge of the subject will reply, if only to point me in a different direction to which I am now facing. Earlier Posting:--- I am sorry if this is more of a stats question than an R-question, but I

Dear R-core/development-team, The problem noted in the subject-line has been a problem in the last three development versions of R for Windows that I have downloaded and tested, the most recent of them being a version I downloaded this morning. Update.packages() using the graphical option, i.e. called as update.packages(ask='graphics', checkBuilt=TRUE) does not work as it should, but

Dear R-users, I am a newbie to this site and a relative new-comer to S/R, so please tread lightly, for you tread... There have been several posting relating to problems with augPred() from the nlme library. Here is a "fix" for one of these problems which may lie at the root of others. In my case the problem with augPred() lay in gsummary(), which augPred() uses, causing it to fail.