Dear all,
I have installed the latest version of R 2.12.0 available on CRAN
(http://cran.r-project.org). When I try to load the recommended
package lattice:
library(lattice)
Error: package 'lattice' is not installed for 'arch=i386'
I am running Rgui using
C:\R\R-2.12.0\bin\i386\Rgui.exe
Dear All
I have problem with this:
IWJR.temp -IWJR.missing[sample(1:length(IWJR.missing),1),]
Error: subscript out of bounds
How I can solved this.
Thanks
IRD
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Sashi Challa wrote:
I am trying to install ncdf package on a Linux 64-bit machine.
I successfully installed netcdf using this command,
./configure --prefix=/home/challar/netcdf/ --disable-netcdf4
Hi Sashi,
Just had a similar issue today. I would suggest the --disable-netcdf4 hasn't
The current version of rgdal on CRAN is 0.6-28, not 0.3-5! The latter
was never in CRAN, but 0.3-7 dates from 2006.
Is your R version (unstated, see the posting guide) also from 2006?
On Sun, 17 Oct 2010, ciccp...@libero.it wrote:
Hi everybody,
I'm trying to install the rgdal package in R,
This is the very first item in the CHANGES file for 2.12.0
How to upgrade your R is an FAQ,
http://cran.r-project.org/bin/windows/base/rw-FAQ.html#What_0027s-the-best-way-to-upgrade_003f
On Mon, 18 Oct 2010, Renaud Lancelot wrote:
Dear all,
I have installed the latest version of R
Thank you. I have read the CHANGES file as well as the FAQ, indeed. I
have uninstalled R, removed the remaining directories and reinstalled
the new version. I have also run
update.packages(checkBuilt=TRUE, ask=FALSE)
Everything works fine except for the recommended package lattice:
I get exactly the same error with package 'RWinEdt' -- though in
http://cran.r-project.org/bin/windows/contrib/r-release/ReadMe it is
mentioned that this has probably something to do with GUI interactions??
Best,
Dimitris
On 10/18/2010 8:49 AM, Renaud Lancelot wrote:
Thank you. I have read
On Mon, 18 Oct 2010, Renaud Lancelot wrote:
Thank you. I have read the CHANGES file as well as the FAQ, indeed. I
have uninstalled R, removed the remaining directories and reinstalled
the new version. I have also run
update.packages(checkBuilt=TRUE, ask=FALSE)
Everything works fine except for
Of course you're right!!! I had a copy of lattice in C:/R/RLIBS
.libPaths()
[1] C:/R/RLIBSC:/R/R-2.12.0/library
I have removed it and this has solved the problem.
Thank you very much.
2010/10/18 Prof Brian Ripley rip...@stats.ox.ac.uk:
On Mon, 18 Oct 2010, Renaud Lancelot
Thank you, but what about dissimilarity matrix objects?
Paul
On Thu, Oct 14, 2010 at 5:30 PM, Peter Langfelder
peter.langfel...@gmail.com wrote:
On Thu, Oct 14, 2010 at 5:21 PM, Paul Rigor (ucla) pr...@ucla.edu wrote:
Hi all,
I would like to use the fpc and cluster packages for
On Mon, Oct 18, 2010 at 3:12 AM, Lorenzo Cattarino
l.cattar...@uq.edu.au wrote:
Thanks for your reply. However, I have only told you half of the story.
My intention was to create a function that outputs 10 different .RData
file. These output files need to contain similar command lines (one
Thanks for the ideas,
Just wanted to say that it was because I was using an old version of R (as
U suggested).
I have now updated to v12.0 and I can see and load reshape2.
(and I agree with Hadley that it would be nice if there was some way of
getting a more informative error message. However
Thanks for your reply. However, I have only told you half of the story.
My intention was to create a function that outputs 10 different .RData
file. These output files need to contain similar command lines (one
apply and one save command) but with slightly different arguments.
This is what I
Hi,
Thanks again for your help with this. I would like to use a variation of
this function in a similar dataset (numeric) with elements separated by a
comma e.g.
dat - read.table(tc - textConnection(
'0,1 1,3 40,10 0,0
20,5 4,2 10,40 10,0
0,11 1,2 120,10 0,0'), sep=)
to simply calculate the
I think that I have successfully installed libreadline.so.6.
But still got an error -- libreadline.so.6 needed during R installation on
SUSE Linux.
help is really appreciated!
===
# ls -lt /usr/local/lib
total 4088
-rw-r--r-- 1 root root 168858 Oct 18 07:15
On 18.10.2010 09:01, Dimitris Rizopoulos wrote:
I get exactly the same error with package 'RWinEdt' -- though in
http://cran.r-project.org/bin/windows/contrib/r-release/ReadMe it is
mentioned that this has probably something to do with GUI interactions??
In that case you upgraded R but
On my Suse Box I have
ste...@gaia:~ ll /usr/lib/libreadline.*
-rw-r--r-- 1 root root 803532 5. Jul 13:46 /usr/lib/libreadline.a
lrwxrwxrwx 1 root root 23 16. Jul 09:59 /usr/lib/libreadline.so -
/lib/libreadline.so.6.1
So libraries are in /usr/lib, not /usr/local/lib
libreadline is
I think there's something about the discrete cell versus centre value
interpretation here, and you are pushing the pixels through R's graphics
engine as well as whatever the png device has to do.
I can't enlighten you about the details of that, but by creating an image
file more directly with
I am getting problems using read.zoo
I have the following data frame
head(anlyNiftyDat[,1:10])
TIMESTAMPACC AMBUJACEM AXISBANK BAJAJ-AUTO BHARTIARTL
BHEL BPCL CAIRN CIPLA
1 2010-01-04 00:00:00 913.60106.10 992.101732.05 325.20
2426.10 650.75 285.50 337.55
2
On Sat, 2010-10-16 at 16:16 -0700, andreas wrote:
I had exactly the same problem with trying to import another .csv file. Turns
out that I was working on a german computer that instead of using a comma
when I saved it as .csv used a semicolon. Just saved it as a normal excel
file, put it on a
Dear All,
I am not an expert about time series, but I am given a time series to
analyze.
That time series stands for the list of individuals in contact with a
given individual at time t_i, where the ID of every individual is an
integer number. (let us not care right now about the meaning of
hello,
yes, thanks a lot - i noticed relevel() beeing very convinient for this
purpose.
having an authority at hand may i kindly ask, if you could reinsure me
that the contrasts below are set up correctly, supposing i want to test
the earlier mentioned hypotheses simultanously.
thanks,
kay
Dear r-help list,
I would like to run a mixed design anova to compare the results from
one population sample to another. Here my within subject variable
(stiulusID) has 45 levels and my between subject variable (group) has
two levels. In addition to my number of levels in the within subject
Hello everyone.
I need some help to understand when number precision in R is set. For this
please consider the following example
for (i in c(2:length(final))){
sizex - c(sizex,(final[i]-final[i-1],digits=2))) # round is used to remove
values that are too small like e-17.
Hi Everyone,
I am trying to install RODBC but I get the following error message
Error in library.dynam(lib, package, package.lib) :
shared library 'RODBC' not found
In addition: Warning message:
package 'RODBC' was built under R version 2.12.0
Error: package/namespace load failed for 'RODBC'
Hello everyone,
could you help me learn if there are any directives that can be used to address
the first and last element of a matrix array?
I would like to thank you in advance for your help
Best Regards
Alex
[[alternative HTML version deleted]]
Hi
r-help-boun...@r-project.org napsal dne 18.10.2010 11:58:06:
Hello everyone.
I need some help to understand when number precision in R is set. For
this
please consider the following example
for (i in c(2:length(final))){
sizex - c(sizex,(final[i]-final[i-1],digits=2))) #
X[1,1] .. first
x[NROW(x),NCOL(x)] for last element and so on
Is this what you need?
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Alaios
Sent: 18 October 2010 15:59
To: Rhelp
Subject: [R] Directive for first and last array
yep. Thanks a lot
Best
Alex
From: Santosh Srinivas santosh.srini...@gmail.com
Sent: Mon, October 18, 2010 1:14:26 PM
Subject: RE: [R] Directive for first and last array arguments
X[1,1] .. first
x[NROW(x),NCOL(x)] for last element and so on
Is this what
Dear List,
I have a table i have read into R:
NameYes/No
John0
Frank 1
Ann 0
James 1
Alex1
etc - 800 different times.
What i want to do is shuffle yes/no and randomly re-assign them to the name.
I have used sample() and permute(), however there is no way to do
On Sun, 17 Oct 2010, Andrew Yee wrote:
I noticed that if I specify the location of TMPDIR in .bashrc as follows on
a Linux 64 bit system:
export TMPDIR=/store/home/ayee/.tmp
I get the following error message when installing R
'installing' or making? Looks like the latter.
make[3]:
On Mon, Oct 18, 2010 at 4:24 AM, Santosh Srinivas
santosh.srini...@gmail.com wrote:
I am getting problems using read.zoo
I have the following data frame
head(anlyNiftyDat[,1:10])
TIMESTAMP ACC AMBUJACEM AXISBANK BAJAJ-AUTO BHARTIARTL
BHEL BPCL CAIRN CIPLA
1 2010-01-04
FAQ 7.31
?all.equal
On Mon, Oct 18, 2010 at 5:58 AM, Alaios ala...@yahoo.com wrote:
Hello everyone.
I need some help to understand when number precision in R is set. For this
please consider the following example
for (i in c(2:length(final))){
sizex -
Hi All,
For those still interested the code submitted by Phil (see below) worked a
treat and produced a vector with the optimal 'nsplit' collated from 50 runs
of the rpart model. I then produced a histogram for the vector called answer
and chose my modal number for nsplits from which I had my
Dear all,
I am trying to fit a model with crossed random effects using lme. In this
experiment, I have been measuring oxygen consumption (mlmin) in bird
nestlings, originating from three different treatments (treat), in a
respirometer with 7 different channels (ch). I have also measured body
Hello everyone.,
I am looking for a covariance function this not the first time I have this type
of problem (to find which function does something). I try in google with R
cran
covariance function but usually this ends with different results that do not
help me that much.
Could you please
Hi Alex,
There are a couple ways. You can install the SOS package, which has
some nice search features. You can also use R's built in
RSiteSearch(some relevant keywords), also, if you think you are
close to what the function should be called, you can try:
apropos(cov)
[1] ability.cov cov
Dear List,
I have a table i have read into R:
NameYes/No
John0
Frank 1
Ann 0
James 1
Alex1
etc - 800 different times.
What i want to do is shuffle yes/no and randomly re-assign them to the name.
I have used sample() and permute(), however there is no way to do
For clarification, the OP downloaded files published before rgdal was
available on CRAN from the deprecated sourceforge site. The OPs references
to the other early development packages (from roughly 5 years ago) are also
from sourceforge. The rgdal and r-spatial projects on sourceforge explicitly
Okay ... now I gotcha ... Thanks
-Original Message-
From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Sent: 18 October 2010 17:19
To: Santosh Srinivas
Cc: r-help@r-project.org
Subject: Re: [R] read.zoo issues
On Mon, Oct 18, 2010 at 4:24 AM, Santosh Srinivas
Hi
I've attached an example about something I want to do in R. This example
was done in a Fortran application called ASGL. Here's an example in
matplotlib
http://matplotlib.sourceforge.net/examples/pylab_examples/hexbin_demo.html
Basically, it's like a scatter plot, but have several
On Oct 15, 2010, at 9:21 AM, ?hagen Patrik wrote:
Dear List,
I each iteration of a simulation study, I would like to save the p-
value generated by coxph. I fail to see how to adress the p-value.
Do I have to calculate it myself from the Wald Test statistic?
No. Look at
Hi all,
I am not exactly fluent in R and I got stuck with this. I would like to
compare each elements of a vector A with any of the elements in Vector B.
For some reasons it does not work.
StartDate = as.Date(01/10/2007, %d/%m/%Y)
TimeSpan = seq(StartDate, by = 'days', length =
Hi,
On Mon, Oct 18, 2010 at 7:37 AM, Peter Francis peterfran...@me.com wrote:
Dear List,
I have a table i have read into R:
Name Yes/No
John 0
Frank 1
Ann 0
James 1
Alex 1
etc - 800 different times.
What i want to do is shuffle yes/no and randomly
On Oct 18, 2010, at 5:23 AM, Paolo Agnolucci wrote:
Hi Everyone,
I am trying to install RODBC but I get the following error message
Error in library.dynam(lib, package, package.lib) :
shared library 'RODBC' not found
In addition: Warning message:
package 'RODBC' was built under R
Hi,
You might consider ?match For example:
StartDate = as.Date(01/10/2007, %d/%m/%Y)
TimeSpan = seq(StartDate, by = 'days', length = 40)
TickLabels = as.Date(c(2007-10-01, 2007-11-01, 2007-12-01,
2008-01-01), %Y-%m-%d)
TimeSpan[1:40] == TickLabels # not designed for comparing multiple
Use %in% instead of '=='
On Mon, Oct 18, 2010 at 10:50 AM, Paolo Agnolucci
agnolucp...@googlemail.com wrote:
Hi all,
I am not exactly fluent in R and I got stuck with this. I would like to
compare each elements of a vector A with any of the elements in Vector B.
For some reasons it does
Dear Ignacio,
if you want it hexagonal (as I gather from the hexbin_demo, have a look at the
hexbin package.
Otherwise, lattice's levelplot is your friend. Or, if you prefer ggplot:
geom_tile or geom_hex.
UIf you play a bit with findFn from package sos, e.g.
findFn (plot 2d density)
findFn
Hi,
Is there a package to perform a sine function fitting to XY data?
Thx,
Ashz
--
View this message in context:
http://r.789695.n4.nabble.com/Sine-function-fitting-tp3000156p3000156.html
Sent from the R help mailing list archive at Nabble.com.
__
Dear R and hydrological community,
The first public (beta) release of two new R packages are now available
on CRAN:
# hydroTSM #
1) hydroTSM is a package for management and analysis of hydrological
time series:
http://cran.r-project.org/web/packages/hydroTSM/
IRD,
It would certainly be more helpful to us if you showed more of your code or
what the dataset IWJR.missing is. Here's a couple of suggestions though.
I can see what you are doing is picking a random sample of this dataset.
Is the IWJR.missing object a data frame or a vector or what?
Dear R users,
The OpenMx developer team takes great pride in announcing the availability
of OpenMx 1.0. The team would like to express its gratitude to the large
number of beta-testers who have helped us improve the code. Thank you.
OpenMx is a free suite of R functions and a estimation back-end
Dear all,
I'm thinking of organising a tutorial on Emacs Speaks Statistics (ESS)
for next year's useR meeting.
http://www.warwick.ac.uk/statsdept/useR-2011/
Tony (Rossini) organised one a few years ago which covered the following
topics:
\begin{enumerate}
\item Introduction (now, 15
Margaret,
I'm not so sure about the pseudo or quasi random value, and without more
information about what you are trying to do I can't provide any more than
the following.
in your R console,
? rnorm
rnorm is the function for pulling random variables from a normal
distribution, and defaults to
Yes - that worked. Thank you. Incidentally I was also comparing a date to a
string which surely didn't help
Thanks again
Ron
On Mon, Oct 18, 2010 at 2:01 PM, Henrique Dallazuanna www...@gmail.comwrote:
Use %in% instead of '=='
On Mon, Oct 18, 2010 at 10:50 AM, Paolo Agnolucci
As announced, the R user conference
useR! 2011
is scheduled for August 16-18, 2011, and will take place at the
University of Warwick, Coventry, UK.
Before the official program, half-day tutorials will be offered on
Monday, August 15.
We invite R users to submit proposals for three hour
Gabor,
You are suggesting some very advanced usage that I do not understand, but it
seems this is not what I meant when I said loop.
I have a df with 47k rows and each of these is fed to a 'predict' which will
output about 62 rows, so the number of groups is very large and I implied that
I
Try this (I think your result in [2,2] is incorrect):
dat - read.table(tc - textConnection(
+ '0,1 1,3 40,10 0,0
+ 20,5 4,2 10,40 10,0
+ 0,11 1,2 120,10 0,0'), as.is = TRUE)
closeAllConnections()
# split the data and create new matrix
newDat - lapply(dat, function(.col){
+ # split by
On Mon, Oct 18, 2010 at 9:55 AM, Bond, Stephen stephen.b...@cibc.com wrote:
Gabor,
You are suggesting some very advanced usage that I do not understand, but it
seems this is not what I meant when I said loop.
I have a df with 47k rows and each of these is fed to a 'predict' which will
Hi!
Using R, I plotted a log-log plot of the frequencies in the Brown Corpus
using
plot(sort(file.tfl$f, decreasing=TRUE), xlab=rank, ylab=frequency,
log=x,y)
However, I would also like to add lines showing the curves for a Zipfian
distribution and for Zipf-Mandelbrot.
It's fairly
Folks,
I have the following dataframe:
x - structure(list(name = c(EU B, EU B, EU B, EU B, EU B,
EU B, AU A, AU A, AU A, AU A, AU A, AU A), date = c(2010-10-11,
2010-10-12, 2010-10-13, 2010-10-14, 2010-10-15, 2010-10-18,
2010-10-11, 2010-10-12, 2010-10-13, 2010-10-14, 2010-10-15,
Strange problem with download.file . for non existent URL an empty file is
created but I am not able to delete the without shutting down R
Example:
download.file(http://test.com/test.txt,test.txt;)
trying URL 'http://test.com/test.txt'
Error in download.file(http://test.com/test.txt;,
Greetings
I want to remove numbers from a string of characters that identify
sites so that I can merge two data frames. For example, a site in one
frame is called 001a Frozen Niagara Entrance whereas the same site in
the other data frame is called Frozen Niagara Entrance. It seems to me
Hi,
I have a list of n items and the ith element has m_i elements within it.
I want to do something like:
predicted.values- lapply(all.predicted.values,'[[',max.growth[[i]])
Where max.growth[[i]] is the element I want to extract from each of the ith
predicted elements. Thus, for example, I
See agrep function:
agrep(Frozen Niagara Entrance, 001a Frozen Niagara Entrance) 0
To remove the numbers:
gsub(\\d, , 001a Frozen Niagara Entrance)
On Mon, Oct 18, 2010 at 12:58 PM, kurt_h...@nps.gov wrote:
Greetings
I want to remove numbers from a string of characters that identify
Try this:
data.frame(split(x$Jem, factor(x$name, levels = unique(x$name))),
check.names = FALSE)
third question:
unstack(x, as.formula(sprintf('%s ~ name', wantedName)))
On Mon, Oct 18, 2010 at 12:47 PM, murali.me...@avivainvestors.com wrote:
Folks,
I have the following dataframe:
x -
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Gabor Grothendieck
Sent: Monday, October 18, 2010 7:03 AM
To: Bond, Stephen
Cc: r-help@r-project.org
Subject: Re: [R]
Try this:
diag(sapply(all.predicted.values, '[[', 'max.growth'))
On Mon, Oct 18, 2010 at 12:59 PM, Gregory Ryslik rsa...@comcast.net wrote:
Hi,
I have a list of n items and the ith element has m_i elements within it.
I want to do something like:
predicted.values-
Unfortunately, that gives me null everywhere. Here's the data I have for
all.predicted.values and max.growth. Perhaps this will help. Thus I want
all.predicted.values[[1]][[4]] then all.predicted.values[[2]][3]] and then
all.predicted.values[[3]][[4]].
I've attached what your statement outputs
On 18-Oct-10 14:58:05, kurt_h...@nps.gov wrote:
Greetings
I want to remove numbers from a string of characters that
identify sites so that I can merge two data frames.
For example, a site in one frame is called
001a Frozen Niagara Entrance whereas the same site
in the other data frame is
Try posting your data using 'dput' so it is easily read for testing.
On Mon, Oct 18, 2010 at 11:17 AM, Gregory Ryslik rsa...@comcast.net wrote:
Unfortunately, that gives me null everywhere. Here's the data I have for
all.predicted.values and max.growth. Perhaps this will help. Thus I want
On 18-Oct-10 15:03:22, Henrique Dallazuanna wrote:
See agrep function:
agrep(Frozen Niagara Entrance, 001a Frozen Niagara Entrance) 0
To remove the numbers:
gsub(\\d, , 001a Frozen Niagara Entrance)
That results in a Frozen Niagara Entrance, which is not
what he said he wants (his
I see that both which(condition) and subset(data,condition)
both treat NA's in condition that same as FALSE's. This
leads people to use those functions for their NA-treating
properties instead of for their main functionality (which
may not be the best way to get things done).
I wonder how much
I want to specify a 2-level mixed model using the lme function in order to
test an a priori hypothesis about the between-group values of the slopes but
don't know how to do this . Here is the problem.
Consider first the case of a single group. The model is: Y_i= a +bX_i +
error where I indexes
Sashi Challa wrote:
I am trying to install ncdf package on a Linux 64-bit machine.
[...]
Hi Sashi,
Just had a similar issue today. I would suggest the --disable-netcdf4
hasn't
been picked up.
Try installing an earlier version and using that path instead.
See:
On Mon, Oct 18, 2010 at 4:38 AM, John Haart anothe...@me.com wrote:
Dear List,
I have a table i have read into R:
Name Yes/No
John 0
Frank 1
Ann 0
James 1
Alex 1
etc - 800 different times.
What i want to do is shuffle yes/no and randomly re-assign them to
Date: Mon, 18 Oct 2010 05:46:14 -0700
From: a...@walla.co.il
To: r-help@r-project.org
Subject: [R] Sine function fitting
Hi,
Is there a package to perform a sine function fitting to XY data?
Since no one replied AFAIK, are you asking about
Hello,
I'm trying to track down more information on hdiffplot than what is supplied
in ?hdiffplot. More specifically, the example code found at the bottom
?hdiffplot (##Compare *three* of them:) is something I'm very interested
in using for my own analysis. However, I don't see how to add a
See this email:
http://www.mail-archive.com/r-h...@stat.math.ethz.ch/msg10849.html
http://www.mail-archive.com/r-h...@stat.math.ethz.ch/msg10849.htmlKevin
On Mon, Oct 18, 2010 at 7:06 AM, anord andreas.n...@zooekol.lu.se wrote:
Dear all,
I am trying to fit a model with crossed random
You probably need mapply since you have 2 list of arguments which you want to
use in sync
mapply(function(x1,x2)x1[[x2]],all.predicted.values,max.growth)
might be what you want.
On Oct 18, 2010, at 5:17 PM, Gregory Ryslik wrote:
Unfortunately, that gives me null everywhere. Here's the data
Hi,
Well I am not completely sure of the R gurus reasons for what they do,
but one explanation is that data is not thrown away unless you ask it
to be. Factors are categorical variables and each level could have
meaning even when there are no cases in it (or particularly when there
are no cases
Hello!
The code below works - if you run it you'll see a stacked area chart
generated based on the data example.
I only have one understanding question about the legend location (the
very last snippet of code):
legend(par()$usr[2],
mean(par()$usr[3:4]),
rev(order.of.vars),
xpd=T,
bty=n,
If you know the period of the sine that you want to fit (just fitting the
amplitude, phase shift, and offset) and are willing to assume normal errors (or
at least normal enough for the CLT) then you can just use the lm function. If
you need to find the period as well (but still willing to
Dear R users,
x-values (EI) = Adw, EG1, LA1, Ad1, LA2, LA3...(14 levels, insect
stages)
y-valus = antpop
within the boxplot function x-values are ordered alphabetically
Idea: x-values ranked by list order (insect stage: Egg stage 1 is
followed by Larvae 1 and not by Egg stage 2 as it
Hello Everyone,
Trying to help someone recover the contents of an SPSS.sav file using
read.spss. This seemed to work well but produced a warning message. My code and
the warning are displayed below.
Spent some time looking for previous questions about this warning. Found a lot
of questions
Sibylle -
Since you didn't give a reproducible example, I'll try
to make one that will illustrate what you need to know to
answer your question.
Suppose we have a variable with levels four,five,
and six.
ff = factor(rep(c('four','five','six'),c(5,7,2)))
Since the table() function will
Hi Everyone,
This is closer to what I need but this returns me a matrix where each element
is a factor. Instead I would want a list of lists. The first entry of the list
should equal the first column of the matrix that mapply makes, the second entry
to the second column etc...
I've attached
Hi,
I was able to reproduce your problem (I changed the date sequence to
starting 01-Oct-2010, but that does not really matter). The
interesting thing is that I could delete some files, and they were not
random. Whenever Windows showed one of the zip folders' size as 0KB,
then I could not
Hello!
If I run this example:
library(nlme)
fm1 - lme(distance ~ age+Sex, Orthodont, random = ~ age + Sex| Subject)
If I run:
summary(fm1)
then I can see the fixed effects for age and sex (17.7 for intercept,
0.66 for age, and -1.66 for SexFemale)
If I run:
ranef(fm1)
Then it looks like it's
Hi,
It seems that the files did not make it through the mailer. Perhaps it didn't
like my extensions. I have now attached the files as .txt's as well as copied
in the contents of each file:
list(list(structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
Michael Sumner-2 wrote:
I think there's something about the discrete cell versus centre value
interpretation here, and you are pushing the pixels through R's graphics
engine as well as whatever the png device has to do.
I can't enlighten you about the details of that, but by creating an
Dmitri:
Not quite sure what you mean by easier ... fixef() and ranef() will
both give coefficients which can be easily manipulated to produce the
results for all subjects.
However, note that there are numerous built-in lme
functions(especially for graphics) that do this internally to produce,
Thank you very much, but not I am not sure now - does ranef(fm1) give
the (total) slope and
intercept values directly for each group or not?
Thanks a lot for clarifying - because I might well have been wrong.
Dimitri
On Mon, Oct 18, 2010 at 5:57 PM, Bert Gunter gunter.ber...@gene.com wrote:
Oh -- I get your question (I think). Not the total, just the random
effects. You have to add them to the fixed effects.
See e.g. p. 39 of Bates and Pinheiro.
-- Bert
On Mon, Oct 18, 2010 at 3:00 PM, Dimitri Liakhovitski
dimitri.liakhovit...@gmail.com wrote:
Thank you very much, but not I am
Yes, sorry for the confusion. Maybe I should have used a different term.
So, I guess, I was right - it gives only the random effects that I
have to add to the fixed effects.
And there is no way to get it done by R (not that I can't do it myself)?
Dimitri
On Mon, Oct 18, 2010 at 6:24 PM, Bert
If you do a search of the archives you will find several similar
questions. As far as I know the answer has generally been unavailable.
It is after all only a warning and you should have gotten uesful
material in the returned value.
--
David.
On Oct 18, 2010, at 2:16 PM, Paul Miller wrote:
Does this do what you want:
x - lapply(seq_along(MaxGrowth), function(.num){
+ AllPredictedValues[[.num]][[MaxGrowth[[.num
+ })
x
[[1]]
[1] 2 1 2 2 2 2 0 2 2 0 0 2 2 0 0 2 2 1 2 0 1 1 0 0 0 2 0 0 0 2 2 0
0 1 0 0 2 0 1 0 2 0 0 2 1 0 0 0 2 1 0 2 2
[54] 2 2 0 2 0 1 0 2 0 1 0 0 0 1 0 0 2
Dear R-helpers ... any thoughts on the below issue ... will help me complete
a small project!
Strange problem with download.file . for non existent URL an empty file is
created but I am not able to delete the without shutting down R
Example:
download.file(http://test.com/test.txt,test.txt;)
On 18/10/2010 8:34 PM, Santosh Srinivas wrote:
Dear R-helpers ... any thoughts on the below issue ... will help me complete
a small project!
Strange problem with download.file . for non existent URL an empty file is
created but I am not able to delete the without shutting down R
Example:
Hi:
One way to permute your sample 1000 times is to use the r*ply() function
from the plyr package. Let s denote your original vector, randomly generated
as follows:
s - rbinom(800, 0.6) # simulates your Yes/No vector
library(plyr)
u - raply(1000, sample(s))# generates a 1000 x
1 - 100 of 105 matches
Mail list logo