Thank you for providing such a complete, self contained example.
I found that 'predict.nlme' does not like a factor in the 'fixed'
argument as you used it, fixed=list(Asym~x1, R0+lrc~1). To see this,
I added 'x1.' as a numeric version of the factor 'x1' and reran it
successfully:
Hello,
I am trying to merge two very large data sets, via
pubbounds.prof -
merge(x=pubbounds,y=prof,by.x=user,by.y=userid,all=TRUE,sort=FALSE)
which gives me an error of
Error: cannot allocate vector of size 2962 Kb
I am reasonably sure that this is correct syntax.
The trouble is that
hi all
i have date and the return series like below, but the dates are not in
uniform intervals. Please show me the way how to create a time series in
'R' so that dates are also associated with the returns.
thanks in advance
Sayonara With Smile With Warm Regards :-)
G a u r a v Y a
Version 0.4.4 of odfWeave is available from CRAN. A Windows binary
should be available shortly.
This version requires base R version 2.3.1 or greater.
Changes from the last version include
- Non-English character sets are handled better. For example, Chinese
characters can be included in R
Which version of R?
Please try 2.4.0 alpha, as it has a different and more efficient
algorithm for the case of 1-1 matches.
On Wed, 6 Sep 2006, Adam D. I. Kramer wrote:
Hello,
I am trying to merge two very large data sets, via
pubbounds.prof -
This is a very common computation in finance.
On the public domain page of the Burns Statistics website
in the financial part is the code and R help file for
'factor.model.stat'. Most of the complication of the code
is to deal with missing values.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525
Dear all.
R 2.3.1, W2k.
I am working with a field trial series where, for the moment, I do
regressions using more than one covariate to explain the protein levels
in malting barley.
To do this I use lme() and a mixed call, structured by both experiment
(trial) and repetition in each
Hi CG,
I think that the best pair of summary plots are
1) the fitted values without random effects against the observed
response variable, and
2) fitted values with random effects against the observed response
variable.
The first plot gives a summary of the overall quality of the fixed
Hi!
Im sorry to bother you but I cant fix this.
I use the lattice function levelplot and I want the colorkey at the
bottom, how do I get it there? I have tried changing colorkey.space and
changing in legend but I cant get it right, plz help
btw I'd like to speceify strings to appear at the tick
Hi,
I am using barplot and would like to know if it is possible to have bars
filled with one color while use a different color for the shading lines.
The following code colors the shading lines, leaving the bars in white:
barplot(1:5, col=c(1:5), density=c(1:5)*5)
while the colors are
Ernst O Ahlberg Helgee wrote:
Hi!
Im sorry to bother you but I cant fix this.
I use the lattice function levelplot and I want the colorkey at the
bottom, how do I get it there? I have tried changing colorkey.space and
changing in legend but I cant get it right, plz help
btw I'd like to
Dear lmer-ers,
My thanks for all of you who are sharing your trials and tribulations
publicly.
I was hoping to elicit some feedback on my thoughts on denominator
degrees of freedom for F ratios in mixed models. These thoughts and
practices result from my reading of previous postings by Doug
Dear list,
I have a list of data.frames (generated by by), that I want to stack into a
single data.frame.
I can do this by cbind, but only by subsetting the list explicitly like this:
cbind(l[[1]],l[[2]],l[[3]],l[[4]])
I find this ugly and not very general.
I tried
cbind(l)
try this:
do.call(cbind, l)
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
Hi everyone!
I have 100 tables of the form:
XCOORD,YCOORD,OBSERVATION
27.47500,42.52641,177
27.48788,42.52641,177
27.50075,42.52641,179
27.51362,42.52641,178
27.52650,42.52641,180
27.53937,42.52641,178
27.55225,42.52641,181
27.56512,42.52641,177
27.57800,42.52641,181
27.59087,42.52641,181
On Thu, 7 Sep 2006, isidora k wrote:
Hi everyone!
I have 100 tables of the form:
XCOORD,YCOORD,OBSERVATION
27.47500,42.52641,177
27.48788,42.52641,177
27.50075,42.52641,179
27.51362,42.52641,178
27.52650,42.52641,180
27.53937,42.52641,178
27.55225,42.52641,181
27.56512,42.52641,177
You can use the 'zoo' or 'its' packages. For 'zoo' see the
documents listed at the end of:
http://cran.r-project.org/src/contrib/Descriptions/zoo.html
On 9/6/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
hi all
i have date and the return series like below, but the dates are not in
uniform
Hello everyone,
I would be glad if you could help out an R-beginner here... I have a
vector of categorial data like this
v - c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4)
When I do
hist(v)
I get the x-axis of the histogram with floating point labels: 1.0, 1.5,
2.0, etc. Is it possible to tell R that the
probably you're looking for a barplot, e.g.,
v - c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4)
plot(factor(v))
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel:
Some of the coordinates might not match and also I do
not have the same number of observations in every
table but I want to get only the common ones back.
This is where it gets tricky!I have tried merge, scan
and every joining function I could find but nothing
seems to do what I want.
the
On Thu, 2006-09-07 at 14:35 +0200, Rotkiv, Rehceb wrote:
Hello everyone,
I would be glad if you could help out an R-beginner here... I have a
vector of categorial data like this
v - c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4)
When I do
hist(v)
I get the x-axis of the histogram with floating
On Thu, 2006-09-07 at 06:18 -0500, Hao Chen wrote:
Hi,
I am using barplot and would like to know if it is possible to have bars
filled with one color while use a different color for the shading lines.
The following code colors the shading lines, leaving the bars in white:
barplot(1:5,
Thanks for your summary, Hank.
On 9/7/06, Martin Henry H. Stevens [EMAIL PROTECTED] wrote:
Dear lmer-ers,
My thanks for all of you who are sharing your trials and tribulations
publicly.
I was hoping to elicit some feedback on my thoughts on denominator
degrees of freedom for F ratios in
I would like to join repeated measures for patients across two visits using
a line. The program below uses symbols to represent each patient. Basically,
I would like to join each pair of symbols.
This is easy in ggplot:
install.packages(ggplot)
library(ggplot)
qplot(visit, var, id=patient,
Hi,
I am currently analysising a counting process form of a cox model allowing for
the inclusion of time dependent covariates. An example model I have fitted is
modlqol-coxph(Surv(Tstart,Tstop,cens.time)~tmt.first+risk
+lqol+cluster(id),data=cat)
summary(modlqol)
My question is quick. I am
Dear Douglas,
I would be happy to re-institute p-values for fixed effects in the
summary and anova methods for lmer objects using a denominator
degrees of freedom based on the trace of the hat matrix or the rank
of Z:X
Please do!
if others will volunteer to respond to the these answers are
Joris De Wolf wrote:
Are your sure your second solution does not work? Try again...
Turns out the second approach did work - but only once I stopped
cutting-and-pasting between two different operating systems (Linux and
Windows under Linux). Apparently, some of the cut-and-paste things I
Dear list,
I have been trying to run the function qvalue under the package qvalue
on a vector with about 20 million values.
asso_p.qvalue-qvalue(asso_p.vector)
Error: cannot allocate vector of size 156513 Kb
sessionInfo()
Version 2.3.1 (2006-06-01)
i686-pc-linux-gnu
attached base packages:
Dear Dr Bates,
Many thanks for such a useful response to my problem.
Regarding Variance Components .
The VarCorr function runs fine for lmer objects once the nlme package is
removed.
Regarding the format of the nested random effects for an lmer object, you
said:
In recent
On Thu, 7 Sep 2006, alex lam (RI) wrote:
Dear list,
I have been trying to run the function qvalue under the package qvalue
on a vector with about 20 million values.
asso_p.qvalue-qvalue(asso_p.vector)
Error: cannot allocate vector of size 156513 Kb
sessionInfo()
Version 2.3.1
DB == Douglas Bates [EMAIL PROTECTED]
on Thu, 7 Sep 2006 07:59:58 -0500 writes:
DB Thanks for your summary, Hank.
DB On 9/7/06, Martin Henry H. Stevens [EMAIL PROTECTED] wrote:
Dear lmer-ers,
My thanks for all of you who are sharing your trials and tribulations
On 9/7/06, Martin Maechler [EMAIL PROTECTED] wrote:
DB == Douglas Bates [EMAIL PROTECTED]
on Thu, 7 Sep 2006 07:59:58 -0500 writes:
DB Thanks for your summary, Hank.
DB On 9/7/06, Martin Henry H. Stevens [EMAIL PROTECTED] wrote:
Dear lmer-ers,
My thanks for all of
Martin Maechler [EMAIL PROTECTED] writes:
DB == Douglas Bates [EMAIL PROTECTED]
on Thu, 7 Sep 2006 07:59:58 -0500 writes:
DB Thanks for your summary, Hank.
DB On 9/7/06, Martin Henry H. Stevens [EMAIL PROTECTED] wrote:
Dear lmer-ers,
My thanks for all of you who
Hi, Patrick: Thanks very much. I'll try it. Spencer Graves
Patrick Burns wrote:
This is a very common computation in finance.
On the public domain page of the Burns Statistics website
in the financial part is the code and R help file for
'factor.model.stat'. Most of the complication of
[EMAIL PROTECTED] asked:
I am trying to divide the columns of a matrix by the first row in the
matrix.
Dividing columns of a matrix by a vector is a pretty fundamental
operation, and the query resulted in a large number of suggestions:
x/matrix(v, nrow(x), ncol(x), byrow = TRUE))
sweep(x, 2, v,
Hello,
I want to colorize the area between two curves, but one of these
curves isn't continuous.
The best solution I found is the 2nd example in the help of polygon,
but how can I get no area filling for the missing data in the 2nd curve.
example:
x1 = c(1:8)
x2 = c(1:8)
y1 =
Dear R People:
Way Off Topic:
Is anyone aware of a website that contains important dates
in statistics history, please?
Maybe a sort of This Day in Statistics, please?
I thought that my students might get a kick out of that.
(actually I will probably enjoy it more than them!)
Thanks for any
Hello Marc Schwartz
On Thu, Sep 07, 2006 at 07:54:05AM -0500, Marc Schwartz wrote:
On Thu, 2006-09-07 at 06:18 -0500, Hao Chen wrote:
Hi,
I am using barplot and would like to know if it is possible to have bars
filled with one color while use a different color for the shading lines.
On Thu, 2006-09-07 at 12:14 -0500, Hao Chen wrote:
Hello Marc Schwartz
On Thu, Sep 07, 2006 at 07:54:05AM -0500, Marc Schwartz wrote:
On Thu, 2006-09-07 at 06:18 -0500, Hao Chen wrote:
Hi,
I am using barplot and would like to know if it is possible to have bars
filled with one
I am trying to install the rgdal package on my Mac OS X 3.9. I
downloaded and installed the GDAL libraries from Fink and then tried to
install rgdal and got the following message. I tried to determine if
the GDAL libraries were in my path but I'm not sure how to do that.
Any ideas? Thanks.
On Thu, 2006-09-07 at 11:57 -0500, Erin Hodgess wrote:
Dear R People:
Way Off Topic:
Is anyone aware of a website that contains important dates
in statistics history, please?
Maybe a sort of This Day in Statistics, please?
I thought that my students might get a kick out of that.
I'm working with an historical image that may be (one of?) the first
uses of gray-scale shading to show the pattern of values in a
matrix/table, later used by Bertin in his 'reorderable matrix'
and sometimes called a scalogram.
The image is at
On Thu, 7 Sep 2006, Prof Brian Ripley wrote:
Which version of R?
Previously, 2.3.1.
Please try 2.4.0 alpha, as it has a different and more efficient
algorithm for the case of 1-1 matches.
I downloaded and installed R-latest, but got the same error message:
Error: cannot allocate vector of
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to
announce our 2-day October 2006 R/S-plus Fundamentals and Programming
Techniques : www.xlsolutions-corp.com/Rfund.htm
*** Washington DC / October 12-13, 2006
*** Seattle Wa / October 19-20
*** San Francisco / October 26-27
One obvious alternative is an SQL join, which you could do directly in
a DBMS, or from R via RMySQL / RSQLite /... Keep in mind that creating
indexes on user/userid before the join may save a lot of time.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf
All:
Is there any way to run an R script file (i.e., *.R) from the command
prompt in the console window. Ultimately, I'm looking to put such code
in a script file so that it can set off other R scripts/programs as
needed.
Thanks.
Marc
[[alternative HTML version deleted]]
On 7 September 2006 at 14:39, Zodet, Marc W. (AHRQ) wrote:
| Is there any way to run an R script file (i.e., *.R) from the command
| prompt in the console window. Ultimately, I'm looking to put such code
| in a script file so that it can set off other R scripts/programs as
| needed.
Which
Using R in batch mode should work on both Windows and Linux:
R CMD BATCH (assuming that R.exe is in your path)
Even without R's location in your path, you could issue the following
command at the prompt (in windows):
c:\Program Files\R\R-2.3.1\bin\R.exe CMD BATCH --vanilla --slave
Dear all,
I'm a newbie to R and I would really apperciate any help with the following:
I have two lists, l1 and l2:
l1:
$A*0101
[1] 0.076 0.109 0.155 0.077 0.09 0 0 0.073
[9] 0.33 0.0034 0.0053
$A*0247
[1] 0 0 0.5 .004 0 0 0
$A*0248
[1] 0 0 0.3 0 0.06
l2:
$A*1101
[1] 0.17 0.24
Hi Andrew,
Thanks a lot, That would give me what I want.
But using my own data and models resulted in this:
plot(fitted(tcos31.c.cp, level=1), FCR.c$g.cp)
Error in xy.coords(x, y, xlabel, ylabel, log) :
'x' and 'y' lengths differ
This is quite correct, as there are some missing values
try something like the following:
lis1 - c(lapply(1:10, rnorm, n = 10))
lis2 - c(lapply(1:10, rnorm, n = 10))
lis1[[5]] - lis2[[8]] - numeric(0)
ind - sapply(lis1, length) 0 sapply(lis2, length) 0
lis1 - lis1[ind]
lis2 - lis2[ind]
Hi,
If I set the p.adjust=none, does it meant that the output p values from the
pairwise.t.test will be the same as those from individual t.tests (set
var.equal=T, alternative=t)?
I actually got different p values from the two tests. See below. Is it supposed
to be this way?
Thanks
Johnny
Are there functions to read image files in jpg, gif or even a pdf file?
Thanks ../Murli
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
I am trying to add text at specific location on my graph. I know this
can be done in R but I can't recollect.
I was trying using locator() to identify the position and use identify()
but I can get it to work. Can someone jog my memory?
Thanks ../Murli
[[alternative HTML version
All,
I'm trying to create an augPred plot in the nlme library, similar to the
plot on
p.43 of Pinheiro Bates (Mixed Effects Models in S and S-Plus) for
their Pixel data.
My data structure is the same as the example but I still get the error
msg below.
comp.adj.UKV - groupedData(adj.UKV ~
no, because the formula for the test statistics ( even assuming that
variances are equal ) of the two different tests are different. in the
pairwise t test, the pairwise differences are
viewed as one sample so it turns into a one sample test. any intro stat book
will have the formulas.
MARK LEEDS wrote:
no, because the formula for the test statistics ( even assuming that
variances are equal ) of the two different tests are different. in the
pairwise t test, the pairwise differences are
viewed as one sample so it turns into a one sample test. any intro stat book
will have
On Thu, 7 Sep 2006, Jonathan Boyd Thayn wrote:
I am trying to install the rgdal package on my Mac OS X 3.9. I
downloaded and installed the GDAL libraries from Fink and then tried to
install rgdal and got the following message. I tried to determine if
the GDAL libraries were in my path
thanks. i assumed we we were talking about the standard textbook difference
between the t test and pairwise t test.
my bad.
- Original Message -
From: Chuck Cleland [EMAIL PROTECTED]
To: MARK LEEDS [EMAIL PROTECTED]
Cc: Li,Qinghong,ST.LOUIS,Molecular Biology [EMAIL PROTECTED];
MARK LEEDS [EMAIL PROTECTED] writes:
thanks. i assumed we we were talking about the standard textbook difference
between the t test and pairwise t test.
my bad.
Notice the difference between paired and pairwise...
- Original Message -
From: Chuck Cleland [EMAIL PROTECTED]
To:
Hello, dear useR,
Is there anyways to get the posterior probabilites for each group by using
knn() instead of only get the proportions of winning class?
Example knn(train=Train[,-c(1:3)], test=Test, cl=group.id.train,k=K, prob=True)
will give you the proportions for winning votes, but I also
Issue this command and then click anywhere on the plot.
loc - locator(1); do.call(text, c(loc, abc))
On 9/7/06, Nair, Murlidharan T [EMAIL PROTECTED] wrote:
I am trying to add text at specific location on my graph. I know this
can be done in R but I can't recollect.
I was trying using
If you don't need borders on the polygons then it can be simply done two points
at a time checking that neither point is an NA:
# data
x1 - x2 - 1:8
y1 - c(1,5,6,1,4,5,5,5)
y2 - c(0,3,3,NA,NA,1,3,4)
# plot
plot(x1,y1,type=l)
lines(x2,y2)
# fill in area between curves with green two points at a
Could somebody program this kind of plot type to R, if none exists,
based on mds or correlation tables or some more suitable method? What
do you think about idea? Does it work? None similar or better exists?
http://weightedassociationmap.blogspot.com/
Atte Tenkanen
University of Turku,
64 matches
Mail list logo