Hi R people!
I have a doubt in lme().
I use this model:
m2.lme-lme(log(cmort)~idade+ano,random=~idade|ano,data=dados)
If i use summary I recive this output:
summary(m2.lme)
Linear mixed-effects model fit by REML
Data: dados
AIC BIClogLik
1139.313 1170.554 -562.6563
Random
On Thu, 31 Mar 2005, Martin Maechler wrote:
BertG == Berton Gunter [EMAIL PROTECTED]
on Wed, 30 Mar 2005 10:09:38 -0800 writes:
BertG If you are on Windows and want to go GUI, see
BertG ?choose.files, ?winMenuAdd, ?winDialog, ?select.list
with the big drawback that it will
you can get the estimated covariance matrix of the random-effects
using:
library(nlme)
m - lme(Orthodont)
# scaled by the residuals variance
pdMatrix(m$modelStruct$reStruct)
# raw
lapply(pdMatrix(m$modelStruct$reStruct), *, m$sigma^2)
However, it seems that in your model you use as a
Is there a better way than:
par(mar=c(6,6,6,6))
plot(1:10,yaxt=n,ylab=)
axis(4)
text(12,5.5,'y-label',xpd=T,srt=90)
to get the y-ticks _and_ the y-label to the rhs of the plot? I did not
find anything in the 'par', 'plot', 'axis' and 'title' manpages to
solve the problem. (the above is ugly,
Hello,
Please check if this is acceptable for you..
par(mar=c(6,6,6,6))
plot(1:10,yaxt=n,ylab=)
axis(4)
mtext('y-label',4,line=2)
Regards,
Carlos.
On Thu, 31 Mar 2005 12:31:49 +0200, joerg van den hoff
[EMAIL PROTECTED] wrote:
Is there a better way than:
par(mar=c(6,6,6,6))
I am working with a matrix (x),
where rows are non-overlapping geographical areas (postcodes; here
denoted with letter row names),
and columns are income bands with a central value (e.g. £7.5k, here
denoted x1,x2, etc.).
The data are counts (how many households with 0-5k pa, 5-10k pa?, etc.)
Dear all:
I am trying to fit a multinomial log linear model to the following data:
worms- data.frame(year= rep(2000:2004, c(3,3,3,3,3)),age=rep(1:3,5),
mud=c(2,5,0,8,7,7,5,9,14,12,8,7,5,13,11),sand=c(4,7,13,4,14,13,20,17,15,23,20,9,35,27,18),
Sir,
I found your description of the dataset about nodal involvement in prostate
cancer. It comes from the book biostatistics casebook. I like to use the
dataset for doing logistics regression. Can you tell me where I can find the
dataset.
Thanks and greetings
Wim van Baarle
[EMAIL PROTECTED]
You can use merge but to do so you will need to define the common key
first. This can be a rowname in the case of a matrix or names in the
case of a vector.
v1 - 1:10
names(v1) - LETTERS[1:10]
v2 - 101:105
names(v2) - sample( LETTERS[1:10], 5 )
merge( v1, v2, by=0, all=TRUE )
Row.names x
emma pilgrim (IGER-NW) wrote:
Hello
I am trying to fit a REML to some soil mineral data which has been
collected over the time period 1999 - 2004. I want to know if the 19
different treatments imposed, differ in terms of their soil mineral
content. A tree model of the data has shown differences
wim van baarle wrote:
Sir,
I found your description of the dataset about nodal involvement in prostate
cancer. It comes from the book biostatistics casebook. I like to use the
dataset for doing logistics regression. Can you tell me where I can find the
dataset.
Thanks and greetings
Wim van
i think it is right!
i rearrage your data in the following form and carry multinom angain,it get the
same result as yours.
multinom(k~year+age,data=worms)
Call:
multinom(formula = k ~ year + age, data = worms)
Coefficients:
(Intercept) year age
sand-106.3405
Dear useRs,
I want to simulate a time series (stationary; the distribution of
values is skewed to the right; quite a few ARMA absolute standardized
residuals above 2 - about 8% of them). Is this the right way to do it?
#
load(rdtb)#the time series
Dear R-experts,
my goal is to visualize the following polynomial regression as a 3D-surface:
Z = b0 + b1*X + b2*Y + b3*XY + b4*X^2 + b5*Y^2
I believe that a solution to this problem may be of interest to a wider
range of scientists because the problem is a derivative of a more general
problem,
Bogdan Romocea wrote:
I want to simulate a time series (stationary; ... snip ...
values is skewed to the right; quite a few ARMA absolute standardized
snip
sim - list(NULL) #simulated
for (i in 1:5) {
sim[[i]] - as.vector(arima.sim(list(ar=c(farma[[coef]][1]),
May be a simple question, but not understanding why in princomp I get different
results for loadings and summary for my eigenvectors and eigenvalues.
When I use pc.cr$loadings using the USArrests dataset the proportion of
variance is equal for each of the components, but when summary(pc.cr)
On Thu, 31 Mar 2005, bogdan romocea wrote:
Dear useRs,
I want to simulate a time series (stationary; the distribution of
values is skewed to the right; quite a few ARMA absolute standardized
residuals above 2 - about 8% of them). Is this the right way to do it?
#
Please note there is no `scatter3d' function in R.
There is one in John Fox's package Rcmdr: please give credit where it is
due.
However, I think you have overlooked functions like persp, image, contour,
cloud wireframe and levelplot (lattice), all of which can plot any
function of two
Dear Johannes and Brian,
scatter3d() in the Rcmdr packages does indeed fit a full quadratic surface
including the product term, which can you verify by setting the argument
model.summary=TRUE in the call to the function. Perhaps the grid over which
scatter3d evaluates the fitted surface is too
On Wed, 30 Mar 2005, Chao Zhu wrote:
But it only gives me the curve with combined treatment. How do I get a
plot with two curves for different treatments. I don't want to do
new-data.frame(treat=c(0,1)) and plot(survfit(fit,newdata=new))
cause it gives you predicted curves for two treatments.
R-folks,
Is there a function, like aggregate, that allows users to bin values?
I've got to break down a data frame into classes of 5cm (or something like
it), and I only know how to do it using code like,
signif - symnum( stems$dbh,
corr = FALSE,
na = FALSE,
Fatal error: cannot find unused tempdir names
I am having the same problem as Michael Epstein (100% of the time) even after
clearing out all the Rtmp* files from the temporary directory. I have tried
reinstalling R, to no avail. Is there some profile setting that is trying to
create temp
On Thu, 2005-03-31 at 09:17 -0800, Jeff D. Hamann wrote:
R-folks,
Is there a function, like aggregate, that allows users to bin values?
I've got to break down a data frame into classes of 5cm (or something like
it), and I only know how to do it using code like,
signif - symnum(
(R 2.0.1)
summary() gives the correct results.
The print.loadings() print method is just giving the squared length of each
of the eigenvectors, which is 1 by definition, of course. I don't know
whether this is a bug or intentional, but it certainly seems silly.
In general, it is good practice
Is the pictex graphics device known not to support color?
In R 2.0.1 Patched (2004-11-17), it produced very pretty output:
## pictex(file = modern-metatheonomy.tex, bg = transparent);
## plotData(data);
## dev.off()
... but it appears to have ignored all color information.
__
Steve
I would like to inquire if there is a function in R that performs the partial
mantel test that controls for the effects of TWO matrices.
The function mantel.partial controls for the effect of one matrix only.
Thank you very much.
-
- Original Message -
From: Thomas Lumley [EMAIL PROTECTED]
To: Chao Zhu [EMAIL PROTECTED]
Cc: r-help@stat.math.ethz.ch
Sent: Thursday, March 31, 2005 11:09 AM
Subject: Re: [R] Cox model qustion in R
On Wed, 30 Mar 2005, Chao Zhu wrote:
But it only gives me the curve with combined
Hi.
I'm looking for a way to plot autocorrelation, but in a little bit
different way than
plot.acf does. Instead of plotting NxN graphs (assuming N is ht enumber of
variables) like plot.acf does, I'd like to have one graph of sum of
all autocorrelations
vs. lag. Is there any function that
Hi
[EMAIL PROTECTED] wrote:
Is the pictex graphics device known not to support color?
In R 2.0.1 Patched (2004-11-17), it produced very pretty output:
## pictex(file = modern-metatheonomy.tex, bg = transparent);
## plotData(data);
## dev.off()
... but it appears to have ignored all color
On Thu, 31 Mar 2005 [EMAIL PROTECTED] wrote:
Is the pictex graphics device known not to support color?
Yes, amongst other settings like lwd and metric info for character output.
If you read the code you will see that the arguments fg and bg are
ignored. In any case, there is no obligation for
Together with Enrique's running start and Prasad's work, we figured out
how to get tick data and bulk data from Bloomberg into R. Here is a code
snippet which builds on Enrique's.
require(RDCOMClient)
blCon - try(blCon - COMCreate(Bloomberg.Data.1), silent=TRUE)
#
Dear all,
I am trying to use hexbin and read the very interesting article on grid
( http://www.ci.tuwien.ac.at/Conferences/useR-2004/Keynotes/Murrell.pdf ) and
am hoping for some advice from more experienced users of hexbin.
I am trying to visualise a data and fit a straight line trough it. For
We would like to announce the availability on CRAN of a new version (1.5)
of our package distr .
-
Changes from 1.4 to 1.5
-package is now using /lazy loading/
-minor changes in the help pages
-minor
Hey guys,
I have problem in dealing with R code for testing the equality of two
independent Binomial Populations. The null hypo. is p1=p2=p (Notice: not only
p1=p2 !). I don't know how to write the test code which includes the value of
'p'. Your help is appreciated!
Shuang
--- jiaotuo niu [EMAIL PROTECTED] wrote:
Hey guys,
I have problem in dealing with R code for testing the equality of two
independent Binomial Populations. The null hypo. is p1=p2=p (Notice:
not only p1=p2 !). I don't know how to write the test code which
includes the value of 'p'. Your
Hi
Adaikalavan Ramasamy wrote:
Dear all,
I am trying to use hexbin and read the very interesting article on grid
( http://www.ci.tuwien.ac.at/Conferences/useR-2004/Keynotes/Murrell.pdf ) and
am hoping for some advice from more experienced users of hexbin.
I am trying to visualise a data and fit a
I'm trying to make my R scripts more readable by others, so I searched for
some R programming conventions and found the following two:
- http://www.ci.tuwien.ac.at/Conferences/useR-2004/Keynotes/Maechler.pdf
- http://www.maths.lth.se/help/R/RCC
the latter being quite extensive. Are there some
Dear all R-helper,
if I want to use stdout from other language as my R program input ,which is
the best way for design the API,using Pipe function or produce a temporary
file,using scan function to read the file ?
thanks
Michael
__
Dear Qian,
Yes, when bootstrap sampling by subject, when a subject is included
in a bootstrap dataset multiple times, you must give that subject
different IDs, or some statistics would be incorrect.
Here's what S+Resample does (from help(bootstrap.args)):
If subject is the name of a
39 matches
Mail list logo