On Sat, 26 Mar 2005, Montse Rue wrote:
I am estimating the following model:
so2.lme-lme(so2~1+I(alcadakm^2)+dia,data=subjectes2,na.action=na.omit)
And when I try to plot the random effects of the quadratic term with
respect to a covariate (mam) I get an error:
so2.lmeRE-ranef(so2.lme,augFrame=T)
I have 100 cases. So i think the dimension is (100, 1).
The PC has a pentium 4 CPU with 512M memory. I don't know if it is enough?
From: Christian Schulz [EMAIL PROTECTED]
To: zhihua li [EMAIL PROTECTED]
CC: r-help@stat.math.ethz.ch
Subject: Re: [R] learning networks with a large number of
i know we can carry out principle component analysis by prcomp() and
princomp().and the result is similar but not identical.
i know ratation in prcomp is the loadings form princomp,sdev from both command
is the standard deviations of the principal components.
and my question is: is the x from
Dear Larry,
dat - data.frame(CAT=sample(c(a,b,c,d),100,rep=T), x=rnorm(100))
tapply(dat[,2],dat[,1, drop = FALSE], mean)
tapply(dat[,2],dat[,1, drop = FALSE], sum)
I hope this helps,
Samuel.
Larry White [EMAIL PROTECTED] wrote:
At the risk of being wacked for asking what should be
I still have 40 invitations! Do not hesitate ;)
This is last remainder on r-help.
--
Lep pozdrav / With regards,
Gregor Gorjanc
University of Ljubljana
Biotechnical Faculty URI:
Hi all:
I am trying to do PCA on the following matrix.
N1 N2A1 A2 B1 B2
gene_a 90110190210290310
gene_b 190210390410590610
gene_c 90110110 90120 80
gene_d 200100400 90
Heck, if anyone cares, I've got 50 or so. Just drop an email. It's
not a bad service.
On Sat, 26 Mar 2005 15:19:04 +0100, Gorjanc Gregor
[EMAIL PROTECTED] wrote:
I still have 40 invitations! Do not hesitate ;)
This is last remainder on r-help.
--
Lep pozdrav / With regards,
Gregor
On Sat, 26 Mar 2005, Srinivas Iyyer wrote:
Hi all:
I am trying to do PCA on the following matrix.
N1 N2A1 A2 B1 B2
gene_a 90110190210290310
gene_b 190210390410590610
gene_c 90110110 90120 80
Dear friends,
R 2.0.1 on windows XP:
TOT - seq(0.01,1,by=0.01)
which(TOT==0.06)
# numeric(0)
which(TOT0.06)
# [1] 6 7 8 etc
which(TOT0.06)
#[1] 1 2 3 4 5
TOT[6]
#[1] 0.06
TOT[6]==0.06
#[1] FALSE
but
TOT[5]==0.05
#[1] TRUE
and
which(TOT==0.05)
# [1] 5
TOT looks as expected when
Troels Ring wrote:
Dear friends,
R 2.0.1 on windows XP:
TOT - seq(0.01,1,by=0.01)
which(TOT==0.06)
# numeric(0)
which(TOT0.06)
# [1] 6 7 8 etc
which(TOT0.06)
#[1] 1 2 3 4 5
TOT[6]
#[1] 0.06
TOT[6]==0.06
#[1] FALSE
but
TOT[5]==0.05
#[1] TRUE
and
which(TOT==0.05)
# [1] 5
TOT looks as
Thanks! I already tried it, but then I have problems with the augPred
values. I get straight lines instead of quadratic lines when plotting
the augPred values.
What can I do?
Montse
On Sat, 26 Mar 2005 09:39:07 + (GMT), Prof Brian Ripley
[EMAIL PROTECTED] wrote:
On Sat, 26 Mar 2005,
i want to know how the predict function to caculate the principal
component,especially when has the newdata argument,i try to find out,but i
fail.anyone knows?thank you!
__
R-help@stat.math.ethz.ch mailing list
On Sat, 26 Mar 2005, Montse Rue wrote:
Thanks! I already tried it, but then I have problems with the augPred
values. I get straight lines instead of quadratic lines when plotting
the augPred values.
What can I do?
WHy are you expecting a non-linear plot from a linear fit?
Montse
On Sat, 26 Mar
Dear all,
I have 5 plots that I would like to include in a single figure, spread over two
rows. If I use mfrow=c(2,3), and produce my plots one after the other, I will
end up with three plots in the first row, and 2 in the second row, which is
what I want. However, I would like the two plots
On Sat, 2005-03-26 at 17:44 +, Hess, Stephane wrote:
Dear all,
I have 5 plots that I would like to include in a single figure, spread over
two rows. If I use mfrow=c(2,3), and produce my plots one after the other, I
will end up with three plots in the first row, and 2 in the second
Hi
Hess, Stephane wrote:
Dear all,
I have 5 plots that I would like to include in a single figure,
spread over two rows. If I use mfrow=c(2,3), and produce my plots one
after the other, I will end up with three plots in the first row, and
2 in the second row, which is what I want. However, I
Does R have a function similiar to the SPlus functoin lmRobMM?
R 2.0.1
Linux
Thanks
John
John Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
Baltimore VA Medical Center GRECC and
University of Maryland School of Medicine Claude Pepper OAIC
University of Maryland School of
Hello,
I am trying to wrap some code that I repeatedly use into a function for
efficiency. The following is a toy example simply to illustrate the problem.
foobar.fun-function(data,idvar,dv){
id.list-unique(idvar)
result-numeric(0)
for (i in id.list){
tmp1-subset(data, idvar == i)
On Sat, 2005-03-26 at 15:27 -0500, John Sorkin wrote:
Does R have a function similiar to the SPlus functoin lmRobMM?
R 2.0.1
Linux
John,
Page 161 of VR's MASS4 suggests using the rlm() function in the MASS
package, which has a 'method = MM' argument.
So:
library(MASS)
?rlm
HTH,
Marc
On Sat, 26 Mar 2005, Hess, Stephane wrote:
HS Dear all,
HS
HS I have 5 plots that I would like to include in a single figure, spread
over two rows. If I use mfrow=c(2,3), and produce my plots one after the other,
I will end up with three plots in the first row, and 2 in the second row,
I tried to replicate in R the example that Pinheiro and Bates present
on section 1.5 of their mixed-effects book, the pixel intensity
example.
fm1Pixel-lme(pixel~day+day^2,data=Pixel,random=list(Dog=~day,Side=~1))
When I run it in R the second order term does not appear in the
estimated model.
On Sat, 2005-03-26 at 15:43 -0500, Doran, Harold wrote:
Hello,
I am trying to wrap some code that I repeatedly use into a function
for efficiency. The following is a toy example simply to illustrate
the problem.
foobar.fun-function(data,idvar,dv){
id.list-unique(idvar)
Apologies for cross posting
-
Salford Systems Data Mining 2005
New York, March 29-30, 2005
Focusing on the Contributions of Data Mining to Solving Real World
Challenges
Two Full
Vicky Landsman wrote:
Thanks to Prof. Ripley, Kjetil and Spencer Graves for help.
I will be more specific.
I have to simulate a bivariate lognormal pair (Y1,Y0) where E(Y1)=X'b,
E(Y0)=X'd, Var(Y1)=c1, Var(Y0)=c0,
X is a data matrix, and b and d are vectors of parameters.
Vicky.
You did'nt specify
On Sun, 2005-03-27 at 14:36 +1200, Brett Stansfield wrote:
Dear R
Is there a way to apply row names as labels to a scatter plot matrix ,
I tried
pairs(dogs, labels=row.names)
Error in strwidth(labels, user) : cannot coerce type closure to character
vector
I'm not sure what this means,
On Sat, 2005-03-26 at 20:36 -0500, John Sorkin wrote:
R 2.0.1
Linux
I am using rlm() to fit a model, e.g. fit1-rlm(y~x). My model is more
complex than the one shown.
When I enter summary(fit1)
I get estimates for the model's coefficients along with their SEs, and
t values, but no p
Hello,dear all:
I want to install the package ordinal,but I don't see the package listed
under package sources.
I try to search it by google,then I found this:
http://euridice.tue.nl/~plindsey/rlibs.html
but the connect does not work.
Where can I found the package ordinal ? Is it still
On Sun, Mar 27, 2005 at 03:26:39AM +0200, Gorjanc Gregor wrote:
Hello!
Has anyone successfully compiled F90 sources in R-package? I found the
same question on r-devel list from 2002 and I wonder if there is any
progress. I heard that g95 and gfortran can be usable.
And one more thing. I
28 matches
Mail list logo