On Mon, 21 Feb 2005, Scholz, Fritz wrote:
I am wondering whether there is a bug in rnorm.
I am wondering if there is a bug in your procedures, or if you used a
very old version of R.
If you call `bug', do have the courtesy to present your reasons. See the
posting guide for other missing
se.contrast computes standard errors of contrasts. It does not compute
the contrast itself, which is fairly simple to do directly.
c1 - c(1,-1)[A]*c(1,-1,0)[B]
sum(c1*score)
[1] -18
se.contrast(fit, as.matrix(c1))
Contrast 1
14.24547
(18/14.24547)^2
[1] 1.596583
so the contrast has value -18,
Corey Bradshaw [EMAIL PROTECTED] writes:
Hello colleagues,
I am attempting to determine the nonlinear least-squares estimates of
the nonlinear model parameters using nls. I have come across a common
problem that R users have reported when I attempt to fit a particular
3-parameter
Is it possible that you've mistyped sum(x-4) instead of sum(x -4)?
Oscar
-Mensaje original-
De: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] nombre de Adaikalavan
Ramasamy
Enviado el: martes, 22 de febrero de 2005 2:24
Para: Scholz, Fritz
CC: R-help
Asunto: Re: [R] rnorm??
1) Why do you
Hi Sebastian. Another option you may want to entertain is the Rdbi
packages from GRASS, http://grass.itc.it/statsgrass/r_and_dbms.html. This
page begins with Postgresql but does have an option for MySQL. I'm running
FreeBSD and unixODBC and never got RODBC running as it should. Also, the
Hi,
I'm still trying to figure out that GLM procedure in SAS.
Let's start with the simple example:
PROC GLM;
MODEL col1 col3 col5 col7 col9 col11 col13 col15 col17 col19 col21 col23
=/nouni;
repeated roi 6, ord 2/nom mean;
TITLE 'ABDERUS lat ACC 300-500';
That's the same setup that I had in my
Bela Bauer wrote:
snip
A quick addition:
I've read https://stat.ethz.ch/pipermail/r-help/2000-November/007457.html
but I really can't get around using those mechanisms because it will
already be quite a transition for users to go from SAS to R...and I
can't give them different results, too!
I've
Hello!
Those of you, who use Sweave a lot, will probably find my shell script
usable. You can get it at:
http://www.bfro.uni-lj.si/MR/ggorjan/programs/shell/Sweave.sh
No warranty, however don't hesitate to contact me if you find an error or
have a patch!
--
Lep pozdrav / With regards,
Bela Bauer [EMAIL PROTECTED] writes:
Hi,
I'm still trying to figure out that GLM procedure in SAS.
Let's start with the simple example:
PROC GLM;
MODEL col1 col3 col5 col7 col9 col11 col13 col15 col17 col19 col21
col23 =/nouni;
repeated roi 6, ord 2/nom mean;
TITLE 'ABDERUS lat ACC
It seems that you want to do is in a Anova with the within factors
factCond and factRoi, right?
If this is correct, then you should try:
summary(aov(vectdata~factCond*factRoi+Error(factCond*factRoi)))
Christophe Pallier
www.pallier.org
Bela Bauer wrote:
Hi,
I'm still trying to figure out that
On Tue, 22 Feb 2005 14:29:00 +0100,
Gregor GORJANC (GG) wrote:
Hello!
Those of you, who use Sweave a lot, will probably find my shell script
usable. You can get it at:
http://www.bfro.uni-lj.si/MR/ggorjan/programs/shell/Sweave.sh
No warranty, however don't hesitate to contact
Hi everybody,
I have the source of a C program that includes some archives .c and some
libraries .h. I'm developing a program using R and I want it to load the C
program I told before. How can I do that? I was looking for some function in
R to do that and I found the .C() but I can't understand
Dear R-list,
*When comparing two logistic regression models with the anova CHi test, I
obtain the following error: (there are no NA's in the time series). How can
this be solved such that I can compare two models on the same dataset were
different explanatory variables are used?
l.KBDI -
R-help,
I am tring to create a plot with two y-axis.
I found an example which is fine but the problem is that the range of
the second y-axes appears in the first y-axes causing confusion.
The example I refer to is :
##
dev.off() ## start with a new graphics
Peter Dalgaard wrote:
However, the real crucial thing here is that SAS is de facto fitting a
mixed-effects model, with random effects being everything with Subject
inside. So, except for the GG/HF business, you should get something
similar if you try
summary(aov(vecData ~ facCond*facRoi +
You need the models to be NESTED for this test to be valid. Look up
`analysis of deviance' in your textbooks.
You get the error because you have two unnested models with the same
number of parameters, so the F numerator is divided by 0. (But this is
merely a symptom of incorrect usage.)
On
these models are not nested and thus LRT should not be used. You get
no p-value because Df=0 in your case. As an alternative you could look
at the deviances.
I hope it helps.
Best,
Dimitris
p.s., you referred to `non-linearity of the data' if you mean that the
probability of positive responses
I am tring to create a plot with two y-axis.
I found an example which is fine but the problem is that the range of
the second y-axes appears in the first y-axes causing confusion. [...]
The key is the ann = FALSE in the second plot():
n = 100
x = rnorm(n)
y = rnorm(n)
z = rnorm(n) * 250
par(mar =
This is what you are trying to get at I think:
dev.off() ## start with a new graphics device
# X11() or postscript()
#par(mar = c(3,3,1,3), oma = c(0,0,0,0), mgp = c(2, 1, 0), bg = white)
par(mar = rep(5,4))
plot(x-rnorm(100),y-rnorm(100), ylab = This is y, xlab = This is x)
z-rnorm(100)*250
I would like to use the 'nice' command when running CPU-intensive R
functions on a server. Will starting R with 'nice R' do the trick?
--
Robert Burrows, PhD
New England Biometrics
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
On a Unix like system you can do `nice +19 R' or perhaps `nice +19 R
CMD BATCH commands.R'.
-roger
Robert Burrows wrote:
I would like to use the 'nice' command when running CPU-intensive R
functions on a server. Will starting R with 'nice R' do the trick?
--
Roger D. Peng
I was not able to reproduce the behavior reported
in my message below. I was trying it on the same laptop
in the same workspace using again R 2.0.0, but the
low counts of 4 to 7 (witnesses by some of my colleagues)
did not show up again. Please do not spend any more
effort on this.
I am
Dear Prof Ripley, Dear Prof Dalgaard
Thank you both for your help. I tried it with helmert contrasts
and got a result that is consistent with lme. I didn't realize
that the parameterization of the model has an influence on the
contrasts that I tried to test.
It seems that I should read a little
Dear colleges,
Does anyone know, how to insert a color bar as used with
filled.contour when using image?
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
On Tue, 22 Feb 2005, Christoph Buser wrote:
Dear Prof Ripley, Dear Prof Dalgaard
Thank you both for your help. I tried it with helmert contrasts
and got a result that is consistent with lme. I didn't realize
that the parameterization of the model has an influence on the
contrasts that I tried to
Hi,
i got a little question about the NA's in vectors and matrices. When I
want to do some operations on a matrix with some NA's it is possible
to stripp them before the computation proceed (by na.rm=TRUE).
But how can I stripp NA's when I want to proceed a 'if', 'for'
or 'while' string?
useRs,
I'm using version 2.0.1 on Windows XP. I am a bit of a newbie and I am
trying to learn the concept of computing on the language. I have an example
that I think ought to work, but will not and I am not sure what I am doing
wrong.
I would like to sort a data frame by a list of columns.
bounds for each group. My question is, is there a function
in R that can do
the same thing for more complex and subtle groupings in
univariate data, and
** provide a statistical basis for the result? **
No. Others have suggested useful ways to **generate** reasonable hypotheses
Hi,
The R platform that I installed on my Windows XP crashes everytime that
I try to run some sophisticated graphics (e.g. Demo Graphics). Is that
to do with the configuration? Shall I reinstall it?
Thanks,
---
Cedric E. Ginestet
PsyPAG Chair
Psychology Department
Thames
Bela Bauer wrote:
I'm still looking for an efficient way to print the new summary. Is
there any easy way to tell the summary or print functions about the
corrected degrees of freedom?
The standard way is not to print the adjusted df, but rather just the
adjusted p-values (along with the
Hi R-gurus,
In my envirovment I have specified the LD_LIBRARY_PATH to
/data/opt/libmy/lib
After starting R the LD_LIBRARY_PATH variable is changed.
Sys.getenv(LD_LIBRARY_PATH)
/data/opt/R-devel//lib/R/lib:/usr/local/lib:/usr/X11R6/lib:/data/opt/libmy/lib
The problem is that in
It might be possible to translate 90% of this into R leaving
just a short shell portion in which case it would be easy
to translate just that part into Windows batch making it
available on that OS too. As someone already mentioned,
using texi2dvi could be helpful too.
Regards.
Date: Tue,
On Tue, 22 Feb 2005 16:50:51 -, [EMAIL PROTECTED] wrote
:
Hi,
The R platform that I installed on my Windows XP crashes everytime that
I try to run some sophisticated graphics (e.g. Demo Graphics). Is that
to do with the configuration? Shall I reinstall it?
You should say what version you
On Tue, 22 Feb 2005 [EMAIL PROTECTED] wrote:
The R platform that I installed on my Windows XP crashes everytime that
I try to run some sophisticated graphics (e.g. Demo Graphics). Is that
to do with the configuration? Shall I reinstall it?
Please consult the rw-FAQ.
It is likely to be a problem
Klaus Ladner wrote:
Does anyone know, how to insert a color bar as used with
filled.contour when using image?
See plot.im() (and im()) in the package spatstat.
cheers,
Rolf Turner
On Tue, 22 Feb 2005, W.E. Wolski wrote:
Hi R-gurus,
In my envirovment I have specified the LD_LIBRARY_PATH to
/data/opt/libmy/lib
After starting R the LD_LIBRARY_PATH variable is changed.
Sys.getenv(LD_LIBRARY_PATH)
/data/opt/R-devel//lib/R/lib:/usr/local/lib:/usr/X11R6/lib:/data/opt/libmy/lib
The
On Tue, 22 Feb 2005 08:41:07 -0800 Berton Gunter wrote:
bounds for each group. My question is, is there a function
in R that can do
the same thing for more complex and subtle groupings in
univariate data, and
** provide a statistical basis for the result? **
No. Others
Hello,
I found that round() does not behave as I expected.
Have you had similar experience as following?
x-seq(0.5,10.5,by=1)
x
[1] 0.5 1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5
round(x)
[1] 0 2 2 4 4 6 6 8 8 10 10
cbind(x,round(x))
x
[1,] 0.5 0
[2,]
This in NOT a bug. It is a convention (the IEEE standard).
Read the help on round(). I.e. RTFM.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
Dongseok Choi wrote:
Hello,
I found that round() does not behave as I expected.
Have you had similar experience as following?
You have found the DOCUMENTED round to even rule, hence not a bug...
Uwe Ligges
x-seq(0.5,10.5,by=1)
x
[1] 0.5 1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5
?round
'round' rounds the values in its first argument to the specified
number of decimal places (default 0). Note that for rounding off a
5, the IEEE standard is used, _go to the even digit_. Therefore
'round(0.5)' is '0' and 'round(-1.5)' is '-2'.
Dongseok Choi wrote:
Hello,
I found that
Why not just do:
do.call(order, A[, ind])
?? A data frame is already a list.
Andy
From: Brough, Tyler (FRS)
useRs,
I'm using version 2.0.1 on Windows XP. I am a bit of a
newbie and I am
trying to learn the concept of computing on the language. I
have an example
that I think
Hi
Hi,
About any graph gallery:
Philippe Grojean and me did have made some work. Our goal was to add a
clip
library to the SciViews project that would offer access to a graph
gallery.
I was workiong on the production of the gallery, where as Philippe is
still
working on his GUI API.
Hi:
I am writing a C program which need a gamma random
number generator. I download the source file of R and
compile, make it myself. There is a rgamma.c
function in the installing directory of
R(/home/zhliu/Backup/R-2.0.1/src/nmath/rgamma.c). My
question is how to call this function in my own
Paul Murrell wrote:
Along similar lines, it would be useful if gallery entries could be
submitted as a plain text or maybe an XML file. Something like:
plot title=blah
descyadda yadda/desc
codeplot(1)/code
/plot
plot title=blah blah
descyadda yadda yadda/desc
codeplot(2)/code
Gabor,
I definitely agree about use of R, however I am not so much in R
as I am in Bash. It took me less than two hours to write this script.
At home I use R under windows and I solve problem of portability
with Cygwin. I know that this is not optimal solution for everyone
but ...
Can you
Perhaps I read more into the question than others, but
I thought he was asking about a stable sort. If that is
the case, ?order says that the sort is stable except for method =
quick.
David L. Reiner
-Original Message-
From: Rau, Roland [mailto:[EMAIL PROTECTED]
Sent: Monday, February
On Tue, 22 Feb 2005, yyan liu wrote:
Hi:
I am writing a C program which need a gamma random
number generator. I download the source file of R and
compile, make it myself. There is a rgamma.c
function in the installing directory of
R(/home/zhliu/Backup/R-2.0.1/src/nmath/rgamma.c). My
question is
Hello:
I was looking for innovative a fun ways of teaching quality control.
Could you suggest me a resource where i could find material that might help
me develop something like this.
Thank you
avneet
(something like the paper helicopter or catapult etc.)
I have no data yet. It is a capital
On Tue, 22 Feb 2005 10:10:58 -0800, Dongseok Choi [EMAIL PROTECTED]
wrote :
Hello,
I found that round() does not behave as I expected.
Have you had similar experience as following?
x-seq(0.5,10.5,by=1)
x
[1] 0.5 1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5
round(x)
[1] 0 2 2
I have been using RODBC for a while with no complaints (R 2.0.1
patched under WinXP), then I saw a link
(http://grass.itc.it/statsgrass/r_and_dbms.html) showing that Rdbi
claims to run a query in 5 seconds that takes RODBC 4.3 minutes. That
was hard enough for me to believe that I wanted to try
Friedrich,
This texi2dvi is also very nice. Dou you know of any pros and cons
in comparison to rubber?
Is tex2dvi also shipped with R for windows? I am not able to find
a binnary in my windows installation, however I found help page.
Did I missed anything? I have 2.0.1
On Tue, 22 Feb 2005
There are two main TeX distributions for Windows, MiKTeX and
fptex. texi2dvi.exe ships with MiKTeX but not fptex.
I think fptex has a texi2dvi UNIX shell script that in principle
could be used with the Cygwin shell on Windows though I have
never tried it. Because of that MiKTeX is really the
You need to give the model formula that gave your output.
There are two sources of variation (at least), within and
between locations; though it looks as though your analysis
may have tried to account for this (but if so, the terms are
not laid out in a way that makes for ready interpretation.
The
Can someone tell me how to install ROracle package to PC? I couldn't find it
CRAN package list. And when I tried to install from zip file, it failed.
Thanks,
Jane
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
I'm sorry if the answer to my problem is buried in the archives. I
have limited experience with R and I couldn't find a solution to my
particular problem. I am running Mac OS X Aqua GUI v1.01 on a new G5
running os 10.3.8 with a 1.8Ghz processor and 1GB of sdram. I just
downloaded
Given a numeric vector of observations, does R have any generic way to estimate
the parameters of commonly used distributions (exponential, gamma, etc.)
without numerically optimizing the likelihood function?
Thanks,
David
___
David R. Bickel
Dear all,
I would like to know that R has the function for garch-t,gjr-
garch,qgarch and egarch.
Best Regards,
Luck
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hi, I have a set of twelve points and wonder how I can
get a function that can then be used to calculate the
area under the curve (most important). Thanks.
Eric
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Dear R users,
in some circumstances, try() shows a strange behaviour,
when the pixmap package is loaded.
The following piece of code works as expected, if it is
either sourced in an interactive session or invoked via
R CMD BATCH (the try-error is printed).
However, if i invoke R using ``R
I'm about to write my thesis in economics and will need to setup and
solve a system of non-linear equations. At our university we usually use
GAMS for this, and though GAMS is a fine program, it bugs me a that I
won't be able to run my code after I finish my thesis without buying a
license for
Thanx very much! I got to load the program into R with no problem!
Talita Perciano Costa Leite
Graduanda em Ciência da Computação
Universidade Federal de Alagoas - UFAL
Departamento de Tecnologia da Informação - TCI
Construção de Conhecimento por Agrupamento de Dados - CoCADa
From: Huntsinger,
Dear all,
I am trying to estimate the point density of a posterior distribution
(this to estimate a bayesian factor in aphylogenetic analysis ). The
problem that I would like to look at several variable at the same
time. From a quick look several of the paramaters that interest me
are
Dear R users,
I am using v2.0.1 on Windows 2000. I have read a .dat file with several
vectors, including 2 factors (2 levels x 3 levels), and a vector of
responses, into R. There are no unique row names. When I try plot(x,y) or
lm(y~x) the following error is returned:
Error in
On 22 February 2005 at 14:58, [EMAIL PROTECTED] wrote:
| side note 2: For make afficionados the follwing 2 rules in combination
| with the Sweave script from the FAQ do almost the same (that's what I
| use :-)
Sure, but it requires a (arguably small) Makefile in every sweave project
directory.
Following the suggestions in the Posting Guide would help us to help you
much better.
What command(s) did you use to get the data into R? What does the `.dat'
file look like? What are `x' and `y'?
Data frames in R, by definition, has to have (unique) rownames. If they are
not present in the
I believe
library(systemfit)
has nlsytemfit function.
Yh
T Petersen wrote:
I'm about to write my thesis in economics and will need to setup and
solve a system of non-linear equations. At our university we usually use
GAMS for this, and though GAMS is a fine program, it bugs me a that I
won't be
Le 22 Février 2005 18:01, Bickel, David a écrit :
Given a numeric vector of observations, does R have any generic way to
estimate the parameters of commonly used distributions (exponential, gamma,
etc.) without numerically optimizing the likelihood function?
Assuming you want to use maximum
Hi,
I use R to do some ARIMA forecasting and R runs out
of memory. The problem is that I have 20160
samples(which are quite alot) and when I try to fit
the model it runs out of memory. I tried with
memory.size() to change the limit, but it wouldn't
work. Is there anything you can suggest? Is
No, this doesn't seem right. What I look for is something that could
solve nonlinear systems with n unknowns and n equations. So there will
be zero degrees of freedom, and statistical methods can't be the right
way forward.
Specifically I can see that the litterature mentions's Scarf's
Greetings!
This is probably a simple problem, but it's really got me flummoxed. I
have some data, basically a vector with a factor, in the following
format, which I read in via 'auth - read.table(authorship.csv,
header=TRUE, sep=,)':
normalized.thomistic,category
1,116.0381739855475264,1
A system of n equations in n unknowns has a unique solution if the
n equations are linear and linearly independent. If the system is
nonlinear, then one must characterize the nonlinearity before saying
anything about whether a solution exists and if so how many solutions
are there?
Hello. Im hoping someone can help with a grouping
question related to the random= statement within the
nlme function. How do you specify that some grouping
levels are inner to others? I tried several things,
given below.
Lets say I have a data frame with five variables,
resp, cov1, ran1, ran2,
Looking at the posting guide will increase the chance to get a helpful
response from this list.
No one knows what kind of operating system are you running: Is it
Windows,MacOS or Linux (32 or 64 bit)?
Memory related problems are reported daily so it could be much
beneficial to browse the
Liu, Jane wrote:
Can someone tell me how to install ROracle package to PC? I couldn't find it
CRAN package list. And when I tried to install from zip file, it failed.
Thanks,
Jane
__
R-help@stat.math.ethz.ch mailing list
roger bos wrote:
I have been using RODBC for a while with no complaints (R 2.0.1
patched under WinXP), then I saw a link
(http://grass.itc.it/statsgrass/r_and_dbms.html) showing that Rdbi
claims to run a query in 5 seconds that takes RODBC 4.3 minutes. That
was hard enough for me to believe that
76 matches
Mail list logo