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 I
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
https://stat.ethz.ch/mailman/listi
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 help
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
This does not really answer your question and frankly
I have not really spent the time to figure out what your
script does but assuming ...\R\rw\bin is in your PATH
then the following one line batch file will run Sweave
from the command line in the Windows console (assuming
its placed in a ba
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?
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
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
algor
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 it
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
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
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 da
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. W
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 model.frame(f
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 correlat
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,
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
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 --van
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
PLEASE
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! http://www.R-proje
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 http://davidbick
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 bioconducte
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
PLEA
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
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 on
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 20
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 so
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(
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 mi
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 i
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
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 p
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:
>
>
> yadda yadda
> plot(1)
>
>
> yadda yadda yadda
> plot(2)
>
>
> This would allow people with existing sets of plots to generat
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 ow
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 GU
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
?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 ro
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
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]
__
R-help@stat.math.ethz.ch
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
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? **
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"
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 [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 wit
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 versio
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, 22
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 /usr/loc
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 \epsilo
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 Vall
> > 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 hy
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.
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? Becau
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 te
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! http://www.R-project.org
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
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 wonder
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
http://www.biostat.jhs
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
https://stat.
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
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 =
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
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 Tu
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 + Error(fac
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 devi
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 <-
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
> 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 hesita
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
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
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,
Grego
Bela Bauer wrote:
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 also
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 la
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
sta
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 yo
"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-par
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,
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 informa
77 matches
Mail list logo