Jose == Jose Quesada [EMAIL PROTECTED]
on Sat, 27 Jan 2007 23:42:34 +0100 writes:
Jose Hi Martin, Thanks for your detailed answer.
Jose x - Matrix(1:12, 3,4, sparse = TRUE)
I hope that you are aware of the fact that it's not
efficient at all to store a dense matrix (it
I apologize if this isn't the right forum for this question. I'm a rookie R
user trying to save a source file and i keep getting this error:
2007-01-27 18:02:12.195 R[175] *** -[NSBigMutableString
writeToFile:options:error:]: selector not recognized [self = 0x68c8c20]
2007-01-27 18:02:12.211
On 27 January 2007 at 16:19, Dirk Eddelbuettel wrote:
| In any event, the bug here is with CGIwithR as it assumes that unzip is the
| real thing. That may be true, but isn't guaranteed.
Sorry, a correction: R2HTML is the one making the assumption about
options(unzip) being unzip, not CGIwithR.
Hello everybody!
I've been working now for quite a while with my R envoirment. However,
today I tried to load it as usal, but I only get the error message
Error in methods:::mlistMetaName(mi, ns) :
the methods object name for 'plot' must include the name of the package
that contains the
Hello,
Is there exist a R Package (or R Code) of nonparametric density
estimation with adaptive methods (k-nearest neighbors,...) for
multivariate data ? I have found the package Locfit but it is only
for the univariate case.
Thank you.
--
Florent BONNEU
LSP
Université Paul Sabatier
118
Does R support grid or parallel computing?
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
Dear Xiaopeng,
There is certainly support for, among others, MPI and PVM; check
packages Rmpi, rpvm, snow, and papply, in CRAN.
Best,
R.
On 1/29/07, xiaopeng hu [EMAIL PROTECTED] wrote:
Does R support grid or parallel computing?
[[alternative HTML version deleted]]
At 21:13 28/01/2007, Bob Green wrote:
Michael,
Thanks. Yes, clearly the volume number for the Schanda paper I cited is wrong.
Where things are a bit perplexing, is that I used the same method as
Peter suggested on two papers by Eronen (referenced below). I can
reproduce in R a similar odds
Dear all,
I have a directory with my research project, containing files
.RData
and
inflow.RData
I am just curious, is there any way to explore contents of inflow.RData from
command line without affecting .RData and without copying inflow.RData to
another location?
I can see names and
You can try copy the file into another location and in the R:
load(file.choose())
ls()
choose the file .RData
On 29/01/07, Vladimir Eremeev [EMAIL PROTECTED] wrote:
Dear all,
I have a directory with my research project, containing files
.RData
and
inflow.RData
I am just curious, is
On Jan 28, 2007, at 9:39 PM, Benilton Carvalho wrote:
This seems to be due to the fact that you didn't have enough memory
when running lmer2.
I might be wrong, but I think Calloc tries to get contiguous
memory, so this might the problem.
If you are positive that you have enough memory,
Thank you, I have used this way already.
I would like to avoid copying files, as I have asked before.
I have found one more solution in Windows.
Just typing inflow.RData (the file name) in command line and pressing enter
runs Rgui, and it loads the file.
Then ls() and str() will give the
Please read the last line of every message to r-help.
On 1/29/07, Shubha Vishwanath Karanth [EMAIL PROTECTED] wrote:
Hi R,
The error I get if I try to merge the zoo objects, intra1, intra2, and
intra3 are as follows:
z=merge(intra1,intra2,intra3)
Error in dimnames(x) - dn : length of
On Mon, 29 Jan 2007, Vladimir Eremeev wrote:
Thank you, I have used this way already.
I would like to avoid copying files, as I have asked before.
attach(inflow.RData)
ls(2)
...
detach(2)
I have found one more solution in Windows.
Just typing inflow.RData (the file name) in command line
Dear All,
I am using R for my research and I have two questions about it:
1) is it possible to create a loop using a string, instead of a numeric vector?
I have in mind a specific problem:
Suppose you have 2 countries: UK, and USA, one dependent (y) and one
independent variable (y) for each
Hi R,
What package of R can I use for Bayesian States Space Modeling? And any
other supporting packages?
Thanks in advance,
Shubha
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
On Mon, 29 Jan 2007 [EMAIL PROTECTED] wrote:
Dear All,
I am using R for my research and I have two questions about it:
1) is it possible to create a loop using a string, instead of a numeric
vector? I have in mind a specific problem:
Suppose you have 2 countries: UK, and USA, one
So, I decided to give it a try (and just now noticed that this is the
example in lmer2)
I just gave it a try on a PPC G4 and it worked as expected. I'm
copying R-sig-mac (sorry for the crosspost) as the experts there
might give you a better suggestion.
fm1 - lmer2(Reaction ~ Days +
Carlo,
try something like:
for (i in c(UK,USA))
{
summ-summary(lm(y ~ x), subset = (country = i))
assign(paste('output', i, sep = ''), summ);
}
(note: it is untested, sorry).
On 1/29/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Dear All,
I am using R for my research and I have two
C.Rosa wrote:
Dear All,
I am using R for my research and I have two questions about it:
1) is it possible to create a loop using a string, instead of a numeric
vector? I have in mind a specific problem:
for (i in c(UK,USA))
output{i}-summary(lm(y{i} ~ x{i}))
In other words, at
That is
C.Rosa wrote:
for (i in c(UK,USA))
output{i}-summary(lm(y{i} ~ x{i}))
for (i in c(UK,USA)) {
lm.txt-paste(output,i,-,lm(,y,i,x,i,),sep=) # 1. produce a
character string containing needed expression
eval(parse(text=lm.txt)) #
2.
Vladimir Eremeev wrote:
That is
C.Rosa wrote:
for (i in c(UK,USA))
output{i}-summary(lm(y{i} ~ x{i}))
for (i in c(UK,USA)) {
lm.txt-paste(output,i,-,lm(,y,i,~,x,i,),sep=) # 1.
produce a character string containing needed expression
eval(parse(text=lm.txt))
Dear All,
Thank you very much for your help!
Carlo
-Original Message-
From: Wensui Liu [mailto:[EMAIL PROTECTED]
Sent: Mon 29/01/2007 15:39
To: Rosa,C
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Loop with string variable AND customizable summary output
Carlo,
try something like:
for
I am returning this to the R-help list. Please keep followups on the list.
Yes, it can be done. It is not currently easy because multcomp doesn't
have the syntax yet. Making this easy is on Torsten's to-do list for the
multcomp package.
See the MMC.WoodEnergy example in the HH package. The
Often you will find that if you arrange your data in a
desirable way in the first place everything becomes
easier. What you really want is a data frame such
as the last three columns of the builtin data frame
CO2 where Treatment corresponds to country and
the two numeric variables correspond to
In thinking about this a bit more here is an even shorter one yet it
does show the level in the Call output. See ?bquote
lapply(levels(CO2$Treatment), function(lev)
eval(bquote(lm(uptake ~ conc, CO2, subset = Treatment == .(lev)
On 1/29/07, Gabor Grothendieck [EMAIL PROTECTED] wrote:
On Mon, 29 Jan 2007, Michael Kubovy wrote:
On Jan 28, 2007, at 9:39 PM, Benilton Carvalho wrote:
This seems to be due to the fact that you didn't have enough memory
when running lmer2.
I might be wrong, but I think Calloc tries to get contiguous
memory, so this might the problem.
If you
And yet one more. This one does not use eval but uses do.call, quote
and bquote instead:
lapply(levels(CO2$Treatment), function(lev) do.call(lm,
list(uptake ~ conc, quote(CO2), subset = bquote(Treatment == .(lev)
On 1/29/07, Gabor Grothendieck [EMAIL PROTECTED] wrote:
In thinking
Hi, Mark
Thanks for the examples- this is great, and has helped me understand alot
more what's going on in the plotting functions.
Now that I'm trying to work error bars into this, I was curious if someone
might give me a hand indexing this properly so that I can format my error
bars to
Look at the 'fame' package I recently put up. You don't need to have the FAME
database installed to use it. Among other things, the package defines a class
tis (Time Indexed Series) that can handle weekly time series.
Wensui Liu [EMAIL PROTECTED] writes:
Monthly and Quarterly ts obj. is easy
Mike,
Using interaction.plot():
with(ex.dat, interaction.plot(x1, x2, response = y1,
type = b, pch = 21,
col = c(red, blue),
ylim = range(ex.dat$y1)))
arrows(1:2, xbar + sem, 1:2, xbar - sem,
In the model:
lm.1 - lm(variable ~ BLOC + TIL * YEAR , data=selvanera)
I found TIL*YEAR interaction significant. Then I am trying to compare
means of the different levels of TIL inside every YEAR using:
mc.2 - glht(lm.1, linfct = mcp(TIL*YEAR=Tukey))
summary(mc.2, test = univariate())
Prior answers are certainly correct, but this is where lists and lapply
shine:
result-lapply(list(UK,USA),function(z)summary(lm(y~x,data=z)))
As in (nearly) all else, simplicity is a virtue.
If you prefer to keep the data sources as a character vector,dataNames,
Hmisc 3.1-2 has been uploaded to the CRAN incoming directory.
Change log
3.2-1 1/25/2007:
Hmisc function 'ecdf' has been renamed 'Ecdf' to deconflict it
with the existing 'ecdf' function in base.
Fixed Bug in format.df that would create numbers with many
trailing
Dear friends,
In configuring R 2.4.1 for Solaris 11, using SunStudio 11 compilers, I
get the following error.
checking readline/history.h usability... no
checking readline/history.h presence... no
checking for readline/history.h... no
checking readline/readline.h usability... no
checking
Or, to throw yet another couple of possibilities into the mix:
lapply(split(YourDF, YourDF$country),
function(x) summary(lm(y ~ x, data = x))
and:
library(nlme)
summary(lmList(y ~ x | country, YourDF))
See ?split and help(lmList, package = nlme)
HTH,
Marc Schwartz
On Mon,
Dear all,
on a practical level an alpha 0 can be found, when a scale is
constructed / evaluated consisting only a few items (say 5) and one
of the items is coded in the wrong direction (values that should
represent a high score wrongfully represent a low score).
Rense
On Jan 24, 2007,
Note that the nlme solution seems to give the same coefficients
but appears to use a single error term rather than one error
term per level of the conditioning variable and that would change various
other statistics relative to the other solutions should that matter.
summary(lmList(uptake ~ conc
I'm not sure - what is the question here? It works for me on a both
PowerPC G5 and Mac Pro (R 2.4.1 CRAN binary):
fm1 - lmer(Reaction ~ Days + (Days|Subject), sleepstudy)
fm1
Linear mixed-effects model fit by REML
Formula: Reaction ~ Days + (Days | Subject)
Data: sleepstudy
AIC BIC
Hi everybody,
I want to do a contourplot in lattice with my raw data overlaid on it(xyplot)
which seemed to be a very easy thing to do.
I've tried it in many ways, but I haven't succeeded at obtaining it.
let's say x, y, and z as the variables that comes from the data frame ''ex'':
Hi all,
I have a file with a dozen weather stations with a dozen weather
variables. I am trying to get a percentile of each variable for each
station. I imported the file into a data frame Data.
To get the percentile for Temperature for Baxley weather station, I used
the following 2
qing wrote:
Dear All,
I am a beginner in learning the package of Biostrings currently and for
practice doing an example.
library(Biostrings);
dnaAlph-new(BioPatternAlphabet,DNAAlphabet(),c(N=AGCT,
+ B=CGT,D=AGT,H=ACT,K=GT,M=AC,R=AG,S=CG,
+ V=ACG,W=AT,Y=CT));
Error in
I'm trying to fit a simple linear regression of just Y ~ X, but both X
and Y are noisy. Thus instead of fitting a standard linear model
minimizing vertical residuals, I would like to minimize
orthogonal/perpendicular residuals. I have tried searching the
R-packages, but have not found
On 1/28/07, Dylan Beaudette [EMAIL PROTECTED] wrote:
Greetings,
I have a groupedData (nmGroupedData) object created with the following syntax:
Soil - groupedData(
ksat ~ conc | soil_id/sar/rep,
data=soil.data,
labels=list(x='Solution Concentration', y='Saturated Hydraulic
On 1/29/07, Kathy-Andrée Laplante-Albert
[EMAIL PROTECTED] wrote:
Hi everybody,
I want to do a contourplot in lattice with my raw data overlaid on it(xyplot)
which seemed to be a very easy thing to do.
I've tried it in many ways, but I haven't succeeded at obtaining it.
let's say x, y,
On Mon, 2007-01-29 at 14:30 -0500, Gabor Grothendieck wrote:
Note that the nlme solution seems to give the same coefficients
but appears to use a single error term rather than one error
term per level of the conditioning variable and that would change various
other statistics relative to the
[Originally sent this to [EMAIL PROTECTED], but in case that's
the wrong list I'm re-posting. Apologies if this becomes a re-post]
---BeginMessage---
I'm trying to fit a simple linear regression of just Y ~ X, but both X
and Y are noisy. Thus instead of fitting a standard linear model
Hi Jonothan,
try the smatr package.
I hope that this helps,
Andrew
On Mon, Jan 29, 2007 at 03:32:35PM -0500, Jonathon Kopecky wrote:
[Originally sent this to [EMAIL PROTECTED], but in case that's
the wrong list I'm re-posting. Apologies if this becomes a re-post]
Date: Mon, 29 Jan 2007
If you are interested in linear Gaussian State Space models, then
package dlm may be of interest.
Giovanni
Date: Mon, 29 Jan 2007 20:50:49 +0530
From: Shubha Vishwanath Karanth [EMAIL PROTECTED]
Sender: [EMAIL PROTECTED]
Precedence: list
Thread-topic: Bayesian States Space Modeling
Dear all,
given I have a data.frame in a format like this
mydf - data.frame(age=rep(1:3,5),
year=c(rep(1996,3), rep(1997,3), rep(1998,3),
rep(1999,3), rep(2000,3)),
income=1:15)
mydf
Now I convert it to some 2D-frequency table like
Please consult the R-admin manual, as the INSTALL file asked you to.
It explains this.
On Mon, 29 Jan 2007, Octavio Tourinho wrote:
Dear friends,
In configuring R 2.4.1 for Solaris 11, using SunStudio 11 compilers, I
get the following error.
checking readline/history.h usability... no
Greetings,
I am pleased to announce the availability of the document, R for SAS
and SPSS Users, at
http://oit.utk.edu/scc/RforSASSPSSusers.doc . It presents an
introductory view of R for people who already know SAS and/or SPSS.
Included are 27 programs written in all three languages (i.e. 81
On Mon, 29 Jan 2007, Roland Rau wrote:
Dear all,
given I have a data.frame in a format like this
mydf - data.frame(age=rep(1:3,5),
year=c(rep(1996,3), rep(1997,3), rep(1998,3),
rep(1999,3), rep(2000,3)),
income=1:15)
mydf
This is what I was hoping for!
Thanks,
Roland
On 1/29/07, Achim Zeileis [EMAIL PROTECTED] wrote:
On Mon, 29 Jan 2007, Roland Rau wrote:
Dear all,
given I have a data.frame in a format like this
mydf - data.frame(age=rep(1:3,5),
year=c(rep(1996,3),
Dear R Community,
Due to maintenance work CRAN will not be available at 30th January around 10
o'clock in the morning CET.
Best regards,
Stefan Theussl
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
Or equivalently:
as.data.frame.table(mymatrix)
On 1/29/07, Achim Zeileis [EMAIL PROTECTED] wrote:
On Mon, 29 Jan 2007, Roland Rau wrote:
Dear all,
given I have a data.frame in a format like this
mydf - data.frame(age=rep(1:3,5),
year=c(rep(1996,3),
Hi,
I have a contour plot and I want to shade a polygon (the area below a line)
but the polygon shading wipes out the contour lines. Does anybody know how
to shade the polygon and still see the contour lines? Thanks.
Jeff
__
R-help@stat.math.ethz.ch
Thank you very much, Weiwei and Jim!
Yeah, I did read the post by Andy, the contributor of this package. It seems
that classwt is not implemented yet. For Weiwei's options, I have a few more
questions. Thanks!
1. try to use rf in fortran by following the linky below
If you look at the examples in ?xyplot.zoo in the zoo package
there is an example of placing a rectangle behind a lattice plot.
Maybe that applies here too?
For classic graphics there is an example of displaying a
rectangle behind a plot here:
http://www.mayin.org/ajayshah/KB/R/html/g5.html
On
Hi, Betty:
1. Fortan code
(http://www.stat.berkeley.edu/~breiman/RandomForests/cc_examples/prog.f)
if(jclasswt.eq.0) then
do j=1,nclass
classwt(j)=1
enddo
endif
if(jclasswt.eq.1) then
c fill in
The fifth option:
actually it might be the easiest way:
you boost your minority by like 10 fold (just repeat each minority
record 10 times). Then run rf on the boosted sample. The learning
process does not exactly behave like using classwt (setting classwt[2]
= 10 will exactly gives weight=10 in
hi,
i have a list of data.frame that has same structure. i would like to know a
efficient way of rbind-ing it.
right now, i write:
n = length(temp) # 'temp' is a list of data.frames
temp2 = data.frame()
for (i in 1:n) temp2 = rbind( temp2, temp[[i]])
return(temp2)
but this is not an
do.call(rbind, temp)
On 1/29/07, jiho.han [EMAIL PROTECTED] wrote:
hi,
i have a list of data.frame that has same structure. i would like to know a
efficient way of rbind-ing it.
right now, i write:
n = length(temp) # 'temp' is a list of data.frames
temp2 = data.frame()
for (i in
Hi Weiwei, thanks a lot for the detailed help!! I tried the option 2 in R.
It works pretty well! You mention that you also implemented RF. Could you
plz share your code with me? Thanks!
Betty
On 1/29/07, Weiwei Shi [EMAIL PROTECTED] wrote:
Hi, Betty:
1. Fortan code (
I'm trying to build a package. The machine is PowerPC G4 with Mac OS 10.4.8,
and I'm using R2.4.1.
I get R CMD build roots working, and it created roots.tar.gz. But I get
the following message when I run R CMD INSTALL -l ../myrlibrary
roots.tar.gz
Hi
here is another approach based on Gabor Grothendieck's idea to add
lines to existing lattice plot.
# based on Gabor Grothendieck's code suggestion
# adds straight lines to panels in lattice plots
addLine- function(...) {
tcL - trellis.currentLayout()
for(i in 1:nrow(tcL))
for(j in
Please update your Xcode tools. According to
http://cran.r-project.org/bin/macosx/RMacOSX-FAQ.html
2.2 or later is needed, and 2.4.1 is current.
On Tue, 30 Jan 2007, C. Lillian Yau wrote:
I'm trying to build a package. The machine is PowerPC G4 with Mac OS 10.4.8,
and I'm using R2.4.1.
I
Hi everybody,
How do we append (note: it's not merging) two zoo objects with the same
column names?
Thanks,
Shubha
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
68 matches
Mail list logo