search for: difford

...))] of Professor Koenker's quantreg program appears to do the same thing. Question: Is there a difference between these two types of plot in the manner in which they depict the relationship between explanatory variables and the response variable ? Thank you inav for your help. Regards, Mark Difford. ------------------------------------------------------------- Mark Difford Ph.D. candidate, Botany Department, Nelson Mandela Metropolitan University, Port Elizabeth, SA.

...e, and redoing the analysis as: glm(Freq ~ as.numeric(Age) + Site, data=glm.table, family='poisson') My reasons for wanting to do this are to be able to include non-linear terms in the model, using say restricted or natural cubic splines. Thank you in advance for your help. Regards, Mark Difford. --------------------------------------------------------------- Mark Difford Ph.D. candidate, Botany Department, Nelson Mandela Metropolitan University, Port Elizabeth, SA.

...2 [3] LC_MONETARY=English_South Africa.1252 LC_NUMERIC=C [5] LC_TIME=English_South Africa.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tools_2.14.0 ----- 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.nabble.com/Graphical-option-to-update-packages-in-development-version-build-of-the-2011-07-31-r56569-for-Windowy-tp3718847p3718847.htm...

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) Mark Difford Ph.D. candidate, Botany Department, Nelson Mandela Metropolitan University, Port Elizabeth, SA.

An embedded and charset-unspecified text was scrubbed... Name: n?o dispon?vel URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080829/57df9fc7/attachment.pl>

...rcept, 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 my TideCode levels now relate to their reference level within Site. Is that correct? Thank you in advance for help. Regards, Mark Difford. Mark Difford Ph.D. candidate, Botany Department, Nelson Mandela Metropolitan University, Port Elizabeth, SA. Send instant messages to your online friends http://uk.messenger.yahoo.com

...ost users of 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 co...

Hi all, I'm trying to do model reduction for logistic regression. I have 13 predictor (4 continuous variables and 9 binary variables). Using subject matter knowledge, I selected 4 important variables. Regarding the rest 9 variables, I tried to perform data reduction by principal component analysis (PCA). However, 8 of 9 variables were binary and only one continuous. I transformed the data by

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

...rcept, 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 my TideCode levels now relate to their reference level within Site. Is that correct? Thank you in advance for help. Regards, Mark Difford. Mark Difford Ph.D. candidate, Botany Department, Nelson Mandela Metropolitan University, Port Elizabeth, SA.

...t.1 <- lme(Chla ~ PO4, random=~1|Site, data=obj) ## or simply: lme(obj) augPred(test.1) plot(augPred(test.1)) (Note that if you are using a data.frame() as your data object you will need to supply a 'primary' statement to augPred(), e.g. augPred(test.1, primary=~PO4). Regards, Mark Difford. --------------------------------------------------------------------- Ph.D. candidate, Botany Department, Nelson Mandela Metropolitan University, Port Elizabeth, SA.

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) )

Hello All, The basic premise of what I want to do is the following: I have 20 "entities" for which I have ~500 measurements each. So, I have a matrix of 20 rows by ~500 columns. The 20 entities fall into two classes: "good" and "bad." I eventually would like to derive a model that would then be able to classify new entities as being in "good

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

...r-project.org When replying, please edit your Subject line so it is more specific than "Re: Contents of R-help digest..." Today's Topics: 1. Help in error of mixed models (Bernardo Rangel Tura) 2. Where to ask general stats questions? (francogrex) 3. Re: structure vs. matrix (Mark Difford) 4. Problems with garch() function tseries package R 2.6.0 (Jos? Augusto Morais de Andrade J?nior - JAMAJ) 5. Re: Where to ask general stats questions? (Patrick Drechsler) 6. help on error message!!! (david csongor) 7. Re: structure vs. matrix (Duncan Murdoch) 8. Re: structure vs. matrix (Duncan Mu...

Hi there. I am still trying to get the RSVGTipsDevice to work, yet I can not. I have copied the first example from RSVGTipsDevice documentation: library(RSVGTipsDevice) devSVGTips("C:\\svgplot1.svg", toolTipMode=1, title="SVG example plot 1: shapes and points, tooltips are title + 1 line") plot(c(0,10),c(0,10), type="n", xlab="x", ylab="y",

Hello, I'm using the R-2.13.1 version in Windows and I'm trying to do a robust Pca with the following: x<-matrix(0.5,30,30) library("rrcov") y<-PcaCov(x) The following error occurs: Error: diff(sv)<0 ist not all TRUE The same error occurs with the other functions. What does this mean and how can I perform the robust PCA with these functions by using a quadratic

Dear all, Is there a way to plot / calculate pointwise confidence bands or interval values for a non parametric regression like sm.regression? Thank you in advance. Regards, Martin

...ultiple 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: http://www.nabble.com/Multiple-comparisons--rank-based-anova-tf4516025.html#a12881037 Sent from the R help mailing list archive at Nabble.com.

Dear Friends, require(cluster) x <- rbind(cbind(rnorm(10, 0, 0.5), rnorm(10, 0, 0.5)), cbind(rnorm(15, 5, 0.5), rnorm(15, 5, 0.5))) plot(pp <- pam(x, 2), which.plots = 1) How can I extract the coordinates used in the plot? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology Postal Address: P.O.Box 400400, Charlottesville, VA 22904-4400