Hi Greg,
Since you haven't yet had a response, you could distill this. It uses the
pixel dataset from nlme() as an example.
## To get separate files, do this
postscript("c:\MyGraph%03.ps", onefile=F)
plot(Pixel, display = "Dog", inner = ~Side, layout=c(4,1))
dev.off()
## To get your layout int
Hi J. Andres,
You are probably better off using the ade4 package, which has two functions
that will do exactly what you want, i.e. a PCA using mixed quantitative and
categorical variables:
## You will need to download ade4 first
library(ade4)
?dudi.hillsmith
?dudi.mix
Regards,
Mark Difford
at 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 Departmen
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
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. ca
uthors have done in package vegan.
There are also some interesting implementations of in package psy, and in
the compositions package.
Best Regards,
Mark Difford.
Jose Claudio Faria wrote:
>
> Dears,
>
> I've been learning biplot (Gabriel, 1971) and I found the function
>
Hi Martin,
Do please, at least, read the documentation for the package you are using!:
?sm.options ## sub: display
## Example
with(iris, sm.regression(Sepal.Length, Sepal.Width, display="se"))
Regards,
Mark Difford.
M. P. Papadatos wrote:
>
> Dear all,
>
Hi Mister_Bluesman,
Sorry to sadden you further, but your example svg works perfectly on my
machine (R 2.5 running on Vista). All the "ToolTips" display. Perhaps it's
your browser [setup]: I'm using Opera 9.20 (Build 8771).
HTH,
Mark Difford.
mister_bluesman wrote:
>
Mister_Bluesman,
Perhaps I should have been more precise: your included svgplot1.svg displays
fine...!
mister_bluesman wrote:
>
> 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:
>
> libr
d lots of "null"
arguments to various functions. It's likely that the problem lies there.
HTH,
Mark Difford.
mister_bluesman wrote:
>
> Ah. Now that's intersting. It works in Opera. But do you get an annoying
> 'Null' label by the cursor when you place it over th
Michael,
Why should that confuse you? Have you tried reading some of the literature
on these methods? There's plenty about them on the Net (Wiki's often a
goodish place to start)---and even in R, if you're prepared to look ;).
BestR,
Mark.
michael watson (IAH-C) wrote:
>
> I'm looking for s
and MDS. If you know of any, of course I would be happy to read them.
>
> Many thanks
> Mick
>
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Mark Difford
> Sent: 14 June 2007 12:49
> To: r-help@stat.math.ethz.ch
> S
Hi Sarah,
You will kick yourself (perhaps?). Try:
?par ## sub: font.main (& font, just above)
## Example
par(mfrow=c(2,1))
plot(1:10, 1:10, main="Figure A", font.main=1)
plot(1:10, 1:10, main=expression("Figure B"))
Regards,
Mark.
Sarah Goslee wrote:
>
> Hello,
>
> I imagine that I'm mi
Hi Arne,
You might also take a look at Prof. Harrell's function:
require(Hmisc) ## req. library [note: also needs lattice]
?smedian.hilow
Then look at:
?xYplot
And, in particular, at: panel.xYplot(), sub "Usage:"
This does what you want.
## Example; also look at Prof Harrell's e
Hi Pablo,
> DF <- data.frame(x=rnorm(20), y=rnorm(20), g1=rep(letters[1:2], 10),
>g2=rep(LETTERS[1:2], each=10),
> g3=rep(rep(letters[3:4],each=5),2))
>
> xyplot(y ~ x | g1 + g2, groups=g3, data=DF)
...
I remember findling with this some time ago and getting most of the way
there. If
Hi Deepayan,
I, and probably quite a few others, will find this very useful until you
find the time to wrap up a proper implementation.
Many thanks,
BestR,
Mark.
Deepayan Sarkar wrote:
>
> On 6/19/07, Juan Pablo Lewinger <[EMAIL PROTECTED]> wrote:
>> I've searched the archives and read the x
Hi Birgit,
Just to add to what Gavin has said. There are two other very powerful
packages in R that handle this kind of thing: ade4 and vegan. Have a
thorough look at both of them. You should be looking at Principal
Coordinate Analysis (Classical Scaling) and Non Metric Multidimensional
Scalin
chool french but anyway, I would like to have a
> look at this lecture notes, but I can not open the link you wrote me.
> Ah sorry now the link works.
>
> Greetings and thanks a lot.
>
> Birgit
>
> Hey Gavin,
> I found it funny that your institute is in the Gower S
Hi James,
Have a look at Cadima et al.'s subselect package [Cadima worked with/was a
student of Prof Jolliffe, one of _the_ experts on PCA; Jolliffe devotes part
of a Chapter to this question in his text (Principal Component Analysis,
pub. Springer)]. Then you should look at psychometric stuff:
andhealth/People/Faculty/Varadhan.html
>
>
>
>
>
>
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Mark Difford
> Sent: Monday, July 02,
EMAIL PROTECTED] On Behalf Of Ravi Varadhan
> Sent: Monday, 2 July 2007 1:29 PM
> To: 'Patrick Connolly'
> Cc: r-help@stat.math.ethz.ch; 'Mark Difford'
> Subject: Re: [R] Question about PCA with prcomp
>
> The PCs that are associated with the smaller eigenvalu
In the absence of a data set, it may help to read the help file carefully:
?hist
Note, in particular, that the argument freq defaults to TRUE "if and only if
breaks are equidistant (and probability is not specified)."
Regards,
Mark.
Sarah Goslee wrote:
>
> Well, how about an example of what
Indeed! And, apropos of the expression, "to be Ripleyed" (and so be
condemned to eating cookies for a long, long time), what about being
"Billasted"?
BestR,
Mark.
Simon Blomberg-4 wrote:
>
> I second the nomination!
>
> Simon.
>
> On Tue, 2007-07-10 at 10:02 -0600, Greg Snow wrote:
>> I nom
Hi Massimo,
Professor Ripley has given you your answer.
It may help you further to know that factor levels aren't automatically
dropped when you subset a data set; you have to do it manually. Some time
ago I scrounged the following command from Andy Liaw's randomForest package:
it removes all e
Hi Alan,
There is a good, ready made, way of doing this, with additional options, in
Professor Fox's car package. See:---
require(car)
?scatterplot.matrix
Cheers,
Mark.
Alan S Barnett-2 wrote:
>
> I'm using pairs() to generate a scatterplot matrix;
>
> pairs(~
> Fuzzy.gray.white.ratio+Fuz
Dear List Members,
I would very much appreciate any pointers you could give me on the following
matter:
Main Question:
To what extent does the "rule" that it is unreasonable to talk about main
effects if there are significant interactions in a model depend upon effect
size [of the significant in
Hi Robert,
Thank you for your reply, and for what appears to be some good/sound
practical advice on managing this kind of problem.
Best Regards,
Mark.
Robert A LaBudde wrote:
>
> Statistical significance is "detectability", and depends upon the
> size of the sample as well as the effect. A
Hi Jean,
You haven't yet had a reply from an authoratitive source, so here is my
tuppence worth to part of your enquiry.
It's almost certain that the "receiving box" is a receptacle into which tags
were placed after they had been drawn and the inscribed measurement noted
down. Measurements on t
Hi GĂ©rald,
I can't help you directly, but you haven't yet had a reply, so...
Googling, as you have found, will waste your time if you know more that you
Google for. Clementine's quite unusual --- in the field of statistical
methods --- so target that. Take the main Stats methods journals and s
ute your lme() statement &c on the object, e.g.:
test.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' statemen
Jorn,
For your model,
model<-lme(Biomass~Age,random=~1|Age/Stand)
think about nesting age in stand (doesn't that makes more sense, anyway ?). If
you're lucky the NaN will zip. So, do
model <- lme( Biomass ~ Age, random = ~ 1 | Stand/Age
I've had a similar problem with unbalanced data when I
iable ?
Thank you inav for your help.
Regards,
Mark Difford.
-----
Mark Difford
Ph.D. candidate, Botany Department,
Nelson Mandela Metropolitan University,
Port Elizabeth, SA.
__
R-help@stat.math.et
Dear Professors Koenker and Varadhan,
Thank you for your detailed and engaging replies. The (very) muddy waters
clear slowly, but only if I keep moving my hands!
Kind regards,
Mark Difford.
Mark Difford
Ph.D. candidate, Botany Department,
Nelson Mandela Metropolitan University,
Port
y 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.
_
34 matches
Mail list logo