To Mark Leeds and the others, thank you for solving my problem, and so quickly,
JM
__
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
and
,
thanks in advance,
Juan Manuel Barreneche,
Zoología de Vertebrados,
Facultad de Ciencias,
UDELAR, Uruguay.
__
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
Let's say I have the following string:
str - P = 0.0, T = 0.0, Q = 0.0
I'd like to find a function that generates the following object from
'str'.
list(P = 0.0, T = 0.0, Q = 0.0)
Thanks!
--
http://mutualism.williams.edu
__
R-help@stat.math.ethz.ch
Hello list,
Is anyone aware of a non-lattice-based alternative to xyplot()?
Thanks!
Manuel
--
http://mutualism.williams.edu
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
On Tue, 2007-07-17 at 22:18 +0100, Gavin Simpson wrote:
On Tue, 2007-07-17 at 16:39 -0400, Manuel Morales wrote:
Hello list,
Is anyone aware of a non-lattice-based alternative to xyplot()?
x - rnorm(20)
y - rnorm(20)
plot(x, y) ?
If you mean some specific aspect of xyplot(), you'll
and use it just like you would use R
command window.
*/Manuel [EMAIL PROTECTED]/* wrote:
Thanks for your answer, but i think i dont do correctly my question.
I need the command line to run Rmdr, like that:
R Rcmdr or R loadRcmdr.R where loadRcmdr has library(Rcmdr
/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
--
Manuel A. Morales
http://mutualism.williams.edu
__
R-help@stat.math.ethz.ch mailing list
https
R
*/Manuel [EMAIL PROTECTED]/* wrote:
Hi to all,
i want to know how can run Rcmdr automatically , or how to load a
library in the call of R
greetings
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch
Hi to all,
i want to know how can run Rcmdr automatically , or how to load a
library in the call of R
greetings
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
(unclass(range(atTicks)),length(atTicks)-1))
Does not help and my question simply is how to make axTicks(1) to be
unclass(atTicks).
Many thanks for any hints.
Manuel
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
of the bars, summary stats,
and CIs). This will allow you to position labels at the x-values of the
plotted bars. For example:
test - bargraph.CI(x.factor = dose, response = len, data = ToothGrowth,
xaxt=n)
axis(side=1,at=test$xvals,labels=c(bar1,bar2,bar3),las=2)
--
Manuel A. Morales
http
?
--
Manuel A. Morales
http://mutualism.williams.edu
signature.asc
Description: This is a digitally signed message part
__
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
) {
r[i] - verycomplicatedsimulation(i)
}
instead of the correct
r - numeric(1)
for(i in 1:1) {
r[i] - verycomplicatedsimulation(i)
}
Would a similar speed issue arise for the construction:
r - vector()
...
--
Manuel A. Morales
http://mutualism.williams.edu
of the
two first columns, genus and species. I am guessing this must be a simple
matter while working with R, but I am stuck at this point.
veg$genus.species - paste(veg$genus,veg$species)
will add a new column combining genus and species.
--
Manuel A. Morales
http://mutualism.williams.edu
but
there was no response... I acknowledge that programming it wouldn't be
so difficult (only needs to apply solve() over a grid), but I am not a
good programmer and I am a little bit concerned with efficiency.
Thank you for any information,
--
José-Manuel Blanco-Moreno
0 0
c 0 0 1 0
d 0 0 1 0
Any suggestions?
Thanks!
--
Manuel A. Morales
http://mutualism.williams.edu
signature.asc
Description: This is a digitally signed message part
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo
is the function lineplot.CI from the package sciplot
(a wrapper that adds SE bars to the output from interaction plots). For
example:
library(sciplot)
data(ToothGrowth)
lineplot.CI(dose, len, group = supp, data = ToothGrowth)
--
Manuel A. Morales
http://mutualism.williams.edu
. The functions in this package replicate some of the
functionality of plotmeans() from the package gplots, with differences
in the treatment of two-way and higher designs. Example graphs can be
seen at:
http://mutualism.williams.edu/sciplot
--
Manuel A. Morales
http://mutualism.williams.edu
- subset(data,subset)
x - eval(substitute(expr), envir=data)
return(mean(x))
}
eg.fn2(a, data, subset=a3)
This returns the error:
Error in eg.fn2(a, data, subset = a 3) :
object a not found
Any suggestions?
Thanks!
--
Manuel A. Morales
http
values.
HTH,
Manuel
On 10/11/06, Tony Plate [EMAIL PROTECTED] wrote:
Here's a way using apply(), and the prob= argument of sample():
df - data.frame(sample1=c(red=400,green=100,black=300),
sample2=c(300,0,1000), sample3=c(2500,200,500))
df
sample1 sample2 sample3
red
a
function for estimating P-values from MCMC samples for factors (or model
comparisons) with greater than 2 levels.
See:http://tolstoy.newcastle.edu.au/R/e2/help/06/09/1003.html
and the very long thread that accompanies it.
HTH,
Manuel
--
Manuel A. Morales
http://mutualism.williams.edu
On Wed, 2006-09-13 at 08:04 +1000, Andrew Robinson wrote:
On Tue, September 12, 2006 7:34 am, Manuel Morales wrote:
On Mon, 2006-09-11 at 11:43 -0500, Douglas Bates wrote:
Having made that offer I think I will now withdraw it. Peter's
example has convinced me that this is the wrong thing
the
significance of a three-level factor? Are there better alternatives than
simply adjusting the CI for the number of factor levels
(1-alpha/levels).
Thanks!
Manuel
=
Manuel A. Morales
Asst. Prof., Biology
Williams College
http://mutualism.williams.edu
__
R-help
Dear list,
Is there any way to load a local data file when connected to a remote
machine via ESS?
Thanks!
Manuel
__
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
=poisson)
sim.model2 - lmer(resp~pred+(1|subject), data=sim.data,
family=poisson)
chisq.sim[i] - anova(sim.model1,sim.model2)$Chisq[[2]]
}
Manuel
On Sun, 2006-08-20 at 11:22 -0700, Spencer Graves wrote:
You've raised a very interesting question about testing a
fixed
this problem. The model and data in
the example are from section 10.4 of MASS.
Thanks!
Manuel
# Load req. package (see functions to generate data at end of script)
library(lme4)
library(MASS)
# Full and reduced models - pred is a factor with 3 levels
result.full - lmer(y~pred+(1|subject), data=epil3
), {
...
})
}
On Wed, 7 Jun 2006, Manuel Morales wrote:
Dear list members,
In some of my functions, I attach the data internally to allow subset
commands or to specify a data frame. This works well except for cases
where there is a masking conflict (which returns a warning). I see
some alternative
On Thu, 2006-06-08 at 11:54 -0400, Manuel Morales wrote:
On Wed, 2006-06-07 at 21:36 +0100, Prof Brian Ripley wrote:
I suggest you investigate with().
Thanks for the suggestion! Unfortunately, it's not clear to me how to
call with() from a function. The example below fails with the error
(response=height, x.factor=voice, trace.factor=range,
data=singer, fun=median, ci.fun=conf.int)
lineplot.CI(response=height, x.factor=voice.part, data=singer,
fun=median, ci.fun=conf.int)
Manuel
On Tue, 2006-06-06 at 00:20 -0300, Mike Lawrence wrote:
Hi all,
I'm
# Function body here
output - x
return(output)
}
Thanks!
Manuel
__
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
,Manuel
Here are the arguments and options I used for comparison purposes,
working on a 1600 records * 15 variables dataset :
# the MART with R tool
lx -
mart( as.matrix(x), y, c(1,1,1,1,1,1,1,1,1,1,1,1,1,2,2)
niter=1000, tree.size=6, learn.rate=0.01,
loss.cri=2
mainly off-list on particular proposals (perhaps sharing information in
the R-Wiki) would be the ideal way to work on this. I am sorry I cannot
invest more time on this at the present moment.
Regards,
Manuel López-Ibáñez.
PS: by the way, I do use Perl, Emacs and LaTeX (almost everyday
Jan T. Kim wrote:
That's an idea I like very much too -- much better than the currently
popular idea of protecting users from the unfriendliness of
programming, anyway...
It is just my opinion that the amount of mail in R-help speaks volumes
about the current friendliness [1], or lack
,
OpenBUGS will not run natively in FC5.
Manuel
On Mon, 2006-05-01 at 13:38 -0500, Paul Johnson wrote:
Dear Jun:
How about telling us which version of Linux you use and how you make
linbugs run? As far as I can tell, the OpenBUGS people have
intentionally removed linbugs from their version 2.2
FullModel-nls(Response~lev[AorB]*exp(Ev*t)*conc/(Km[AorB]+conc),data=list(Response=yBoth,t=tBoth,conc=concBoth),start=list(lev=c(90,70),Ev=-0.5,Km=c(0.8,0)),trace=TRUE)
## but i get to a singular gradient
Any other pointers,
thanks
Manuel
--- Prof Brian Ripley [EMAIL PROTECTED]
escribió
Barry Rowlingson wrote:
The frame-based nature of the CRAN pages is slightly problematic,
since you click on a menu item and the URL doesn't change. Hence there's
no way to send someone a URL that gives them the same view as you'd get
if you go to the home page and then click on
Is it possible to include a factor in an nls formula?
I've searched the help pages without any luck so I
guess it is not feasible.
I've given it a few attempts without luck getting the
message:
+ not meaningful for factors in:
Ops.factor(independ^EE, a)
This is a toy example, my realworld case is
[a==0]^2.1)+3
library(nlme)
gnls(respo~independ^b+a,start=list(b=1.8))
Many thanks.
Manu
--- Andrew Robinson [EMAIL PROTECTED]
escribió:
Manuel,
I don't think that it works very easily. Instead,
try gnls() in the
nlme package.
Cheers
Andrew
On Thu, Apr 20, 2006 at 11:18:02AM +0200
zhang jian wrote:
* n=PIKO[status=snag,]
Error in [.data.frame(PIKO, status = snag, ) :
unused argument(s) (status ...)
= is the assignment operator, you should use == as in:
n=PIKO[status==snag,]
__
LLama Gratis a
in the first place. For example, it would be interesting to know
how many new R users don't know that there is a manual page for [.
I hope you can understand my point of view (you may disagree, though.)
Regards,
Manuel.
John Fox wrote:
Dear Manuel,
Although ?sample doesn't specifically
help(sample) does not say anything about randomly permuting the rows of
a matrix M by using M[sample(m,m),]. Perhaps it could be added as an
example of use.
John Fox wrote:
Dear Jose,
M[sample(m, m),] will randomly permute the rows of M. [You probably could
have figured this out via
and errors because the function calculates
the mean and CI for you.
HTH
Manuel
__
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
Hi Gašper
On Wed, 2006-02-22 at 14:12 +0100, Gasper Cankar wrote:
Hello everyone.
For reasons too long to explain I wanted to do plots similar to histograms
with plot(type=h).
I ran into a problem - if I set line width too high, histogram isn't accurate
anymore.
For example:
. Is there
something else I need to do? I am using a windows
based machine to do this.
Try xyplot(..., par.settings=list(background=white))
HTH,
Manuel
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
On Thu, 2006-01-26 at 21:55 -0500, Duncan Murdoch wrote:
On 1/26/2006 9:45 PM, Manuel Morales wrote:
Dear list members,
I'm looking for a way to write nested functions similar to the
function Nest or NestList in Mathematica.
E.g.,
f-function(x) x+2*x
f(f(f(2)))
might
Dear list members,
I'm looking for a way to write nested functions similar to the
function Nest or NestList in Mathematica.
E.g.,
f-function(x) x+2*x
f(f(f(2)))
might instead be written as nest(f, 2, 3)
read as, nest function f 3 times with 2 as the initial value.
Thanks!
Manuel
: 0x88a076c
But I can not find what is form.
Any help, please.
Manuel
Manuel Gutierrez wrote:
The option se.fit in predict.nls is currently
ignored.
Is there any other function available to calculate
the
error in the predictions?
Thanks,
Manuel
The option se.fit in predict.nls is currently ignored.
Is there any other function available to calculate the
error in the predictions?
Thanks,
Manuel
__
LLama Gratis a cualquier PC del Mundo.
Llamadas a fijos y móviles desde 1
to calculate the errors in the
predictions from nls?
Thanks,
Manuel
__
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
to calculate the errors in the
predictions from nls?
Thanks,
Manuel
__
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
work overload.
Thanks in advance for your patience and your good efforts!
Best wishes,
Manuel Castejón Limas
E-mail: manuel.castejon at unileon.es
Área de Proyectos de Ingeniería
Departamento de Ingeniería Eléctrica y Electrónica
Universidad de León
Spain
PS: You may read more about moodle
[names(count),1,1]-count[names(count)]
Fails with Error: indexing outside limits.
Many thanks for any help
Manuel
__
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
[as.numeric(names(count))]-count
rwname-rownames(data)
for (i in 1:dim(data)[1]) data[i,1,1]-temp[as.numeric(rwname[i])]
What works for me but I am convinced there is a far more elegant way.
Kind regards
Manuel
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf
characters sometimes seem to jump around
to other places. Funny bug, but definitely nothing to do with R.
Best regards
Manuel
-Ursprüngliche Nachricht-
Von: Henrik Bengtsson [mailto:[EMAIL PROTECTED]
Gesendet: Montag, 29. August 2005 15:43
An: Peter Dalgaard
Cc: Schneider, Manuel; r-help
and many thanks for any
suggestion.
Kind regards
Manuel
__
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
Someone knows if it is possible to perform a path analysis with sem package (or
any other) to explain a dependent *dichotomus* variable?
Thanks,
Manel Salamero
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
that usually appears in the console and
then crashes. Rterm.exe works fine.
I've now unistalled R 2.1.0 and installed R 2.1.1 with no effect, still
console is not found.
Any clues on this?
Best regards
Manuel
°°°
Manuel Schneider
Eawag
Environmental chemistry
Ueberlandstr. 133
8600 Dübendorf
Phone
-minimize.fn(result$parms)
This gives the results:
result$parms
[1] 1.184172e-05 4.878992e-03 3.045663e+00
result$SS
[1] 0.005436355
Which is very similar to what you want.
HTH,
Manuel
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman
to increase the amount of text that I can paste into an R
session?
Thanks!
Manuel
__
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
Hi Meredith,
When I've wanted to do this, I put all my files (group1.txt,
group2.txt ...) in a separate directory. Then, running R from that
directory:
files-list(files)
for (i in 1:length(files)) {
group-read.table(files[i],header=T)
do stats here
}
On Fri, 2005-04-29 at 15:17 +1000, Briggs,
On Fri, 2005-04-22 at 15:42 -0500, Marc Schwartz wrote:
On Fri, 2005-04-22 at 15:44 -0400, Manuel Morales wrote:
On Fri, 2005-04-22 at 17:48 +0200, Uwe Ligges wrote:
Waichler, Scott R wrote:
I installed R-2.1.0 from source on a Linux box running Red Hat
Enterprise Linux WS
guessing you
don't mean that X11 is not running (I use Gnome for my desktop).
Manuel
__
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 R-users
I need to convert a matrix with three columns in a new array with
multiple columns.
For example,
oldmatrix
1 4 5
1 54 52
1 9 43
2 32 5
2 54 6
2 76 6
3 54 54
3 543 7
3 54 6
and in that case any solutions to obtain a
reliable linear model.
Thanks,
Manuel
Some more detailed information:
A-lm(y~x1+x2)
summary(A)
Call:
lm(formula = y ~ x1 + x2)
Residuals:
Min1QMedian3Q Max
-4.221946 -0.484055 -0.004762 0.397508 2.542769
Hello list members,
I'm trying to get nlme to work with lsoda from odesolve. Currently, I
can use nlsList to fit a simple logistic growth model to some simulated
data, but I get the following error message with nlme: Error in
model.frame(formula, rownames, variables, varnames, extras,
for me. A toy example of
linear.hypothesis for my case would be great.
Thanks,
Manuel
--- Peter Dalgaard [EMAIL PROTECTED]
escribió:
John Fox [EMAIL PROTECTED] writes:
Dear Vito,
Since Manuel says that he wants to obtain a test
and not obtain two
tests, I assume that he's interested
and for x2 different than 0.3
Thanks
Manuel
__
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
will have rare cases or not
or where they will be located.
Thanks,
Manuel
__
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
I've a linux system with 2Gb of memory which is not
enough for reading a 446Mb netcdf file using ncdf:
library(ncdf)
ncold - open.ncdf(gridone.grd)
Error: cannot allocate vector of size 1822753 Kb
When I look at the free memory in my system I can see
that none of the Swap space is being used by
,
Manuel
--
-
Manuel Metz
Sternwarte der Universitaet Bonn
Auf dem Huegel 71 (room 3.06)
D - 53121 Bonn
E-Mail: [EMAIL PROTECTED]
Phone: (+49) 228 / 73-3660
Fax:(+49) 228 / 73-3672
__
R-help@stat.math.ethz.ch
other
solutions welcomed) to be able to do this but,
shouldn't I get a more gentle warning than the kill
message?
Thanks,
Manuel
platform i686-pc-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major2
minor
I could not find any help pages on How to test many
objects for being of equal length
Something like identical for more than two objects?
x-1:6
y-1:10
z-3:5
## For two objects I can do:
identical(length(x),length(y))
## For more than two I currently can do:
))
Error in tapply(object[[nm]], groups, FUN[[numeric]], ...) :
arguments must have same length
Thanks!
Manuel
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org
I have an openoffice spreadsheet with a column of
character strings.
Some of them contain accents.
I want to read it in R so I have saved it as a csv
file using Western Europe (ISO-8859-1) character set
(the default, I've tried other sets but it doesn't
help).
R reads it fine with
of a
UTF-8 locale, which R
does not currently support and which some Linux
distros have made their
default. However, R does issue a warning -- so did
you get one?
On Thu, 25 Nov 2004, Manuel Gutierrez wrote:
I have an openoffice spreadsheet with a column of
character strings.
Some
correlation across days. I have Pinheiro and Bate's extremely
useful book, but I can't seem to find anything that addresses this. I'm
using R 1.9.1 with version 3.1-50 of nlme.
Thanks very much for any help. I've included my sample code below.
Manuel Morales
#Below I generate a day variable (note two
?
Thank you very much.
Manuel.
[1] http://finzi.psych.upenn.edu/R/Rhelp02/archive/12926.html
[2] http://finzi.psych.upenn.edu/R/Rhelp02/archive/32849.html
inline: prueba.png__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
intervals are much bigger than the ones
calculated by STATGRAPHICS, thus I think I did something wrong.
I am not an expert in Statistics nor in R, thus if anyone has any
suggestion...
Thank you very much,
Manuel.
[[alternative HTML version deleted
?
I'm using R version 1.8.
Thanks,
Manuel
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
20 475.9776 542.7283 -217.9888 1 vs 2 44.17249 .0001
Why are my within-subject errors so large, and why doesn't intervals() work
for the autoregressive errors example in the book. I am using R version 1.8
on Windows 2000. Any insight would be greatly appreciated!
Manuel
Manuel A. Morales
to R, have the book Mixed
Effects Models in S and S-Plus, but can't seem to get this analysis to
work. Any suggestions?
Thanks!
Manuel
Data:
SUB GROUP DV TRIAL
1 1 3 1
1 1 4 2
1 1 7 3
1 1 3 4
2 1 6
Hi everyone,
does anyone know a way to calculate VIP (Variable Importance in the
Projection) from R pls.pcr package procedures?
Thanks,
--
Manuel Martin
23 rue des gâtines
75020 Paris tel: 0143664965
_
Envie de discuter
dilut 4 4 420.79.0001
sample*dilut 20 20 1.610.1486
May be possible obtain the same results with NLME function?
I would appreciate any kind of help.
Best regards,
Manuel Ato
University of Murcia
in advance.
Manuel Ato
Dpto. Psic.Básica y Metodología
Apartado 4021
30080 MURCIA (Spain)
e-mail: [EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
I would like to use Rterm but i don`t know its parameters. I have searched about this
issue but i haven´t found anything. thank you.
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Dear Rohan,
Have a look at the fdim library, it may be of interest to you as far as
fractal dimension
(or box counting if you prefer) is concerned.
Best wishes,
Manuel Castejon
__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman
86 matches
Mail list logo