On 07/04/04, Duncan Murdoch ([EMAIL PROTECTED]) wrote:
S?ren H?jsgaard wrote:
I'm on a PC and have VIM, a VI clone that I love. Is it possible and
relatively easy to replace the Emacs editor with one like that which
you already have?
Does VIM have a way to run a program that uses stdin,
That is because the Kalman filter is only an estimation method. You will have to
specify the
matrices A and C in:
y_t = A_t * x_t + noise_t
x_t = C *x_t-1 + noise_t
where you plug the explanatory variables in A_t. The Kalman filter will then give you
estimates
of the state vector
Dear R-helpers,
I'm a beginner using R 1.8.1 on Windows 2000. I'm trying to replicate some
examples in Franses' Paap's Quantitative Models in Marketing Research.
t - 1:1000
e1 - rnorm(1000)
e2 - rnorm(1000)
x - 0.0001*t+e1
y2 - -2+x+e2
y - ifelse(y20,1,0)
plot(x, y, pch = 16, col =
Jan P. Smit wrote:
Dear R-helpers,
I'm a beginner using R 1.8.1 on Windows 2000. I'm trying to replicate some
examples in Franses' Paap's Quantitative Models in Marketing Research.
t - 1:1000
e1 - rnorm(1000)
e2 - rnorm(1000)
x - 0.0001*t+e1
y2 - -2+x+e2
y - ifelse(y20,1,0)
plot(x, y, pch =
You probably mean something like:
ti - 1:1000
e1 - rnorm(1000)
e2 - rnorm(1000)
x - 0.0001*ti+e1
y2 - -2+x+e2
y - ifelse(y20,1,0)
plot(x, y, pch = 16, col = darkblue,
main = expression(paste(Scatter diagram of ,
italic(y[t]), against , italic(x[t]))),
xlab =
I have a problem with nlme on Windows 2000, and I'm having a devil of a
time determining whether the problem is with my computer or with
something in R. I'm running v1.8.1 on a Dell Pentium III with 512MB of
RAM and all of the recommended Windows 2000 updates applied.
If I use Rterm, I can run
Ed L Cashin [EMAIL PROTECTED] writes:
Julian Taylor [EMAIL PROTECTED] writes:
...
You are better off using other tools to give you the right subsets. Try
d - do.call(rbind, lapply(split(d, factor(paste(d$a, d$b, sep =
))),
function(el) el[el$c ==
Richard == Richard A O'Keefe [EMAIL PROTECTED]
on Thu, 8 Apr 2004 16:48:05 +1200 (NZST) writes:
Richard I just got a G4 running MacOS 10.3 today, so I
Richard immediately downloaded R from CRAN.
good move; one of the first thing a new computer needs,
(just after emacs) ;-)
Martin Maechler [EMAIL PROTECTED] writes:
Richard a designated sysadmin, which I'm not. Is there any
Richard alternative to downloading a source distribution
Richard and doing the classic ./configure;make;make install
Richard waltz?
yes, there's always been.
You have
I am an applied statistician new to R, and I was wondering if R can calculate expected
mean squares for unbalanced mixed models. Any help would be appreciated.
Thanks,
Bill Moore
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED]
Hie!
I do not manage to make several plots on a same graph ( for example two
densities on the same graph )
thanks for the answer
Pierre Clauss
xxx
Ce message et toutes les pièces jointes (ci - apres le \message\ )
sont établis a l'attention
Hello!
This works:
hist(rnorm(100, mean = 20, sd =12), xlim=range(0,100), ylim=range(0,50))
par(new = TRUE)
hist(rnorm(100, mean = 88, sd = 2), xlim=range(0,100), ylim=range(0,50))
Matthias Templ
-Ursprüngliche Nachricht-
Von: Pierre Clauss [mailto:[EMAIL PROTECTED]
Gesendet:
Douglas Bates [EMAIL PROTECTED] writes:
...
The point is that factor(paste(...)) returns a factor with length of
nrow(d) and with factor levels determined by the combination of levels
of a and b (provided that you don't get ambiguities as described
above). The split function does not require
Moore, William G [EMAIL PROTECTED] writes:
I am an applied statistician new to R, and I was wondering if R can
calculate expected mean squares for unbalanced mixed models. Any
help would be appreciated.
Not that I know of but you can obtain REML or ML estimates of the
parameters in an
Ed L Cashin [EMAIL PROTECTED] writes:
Douglas Bates [EMAIL PROTECTED] writes:
...
The point is that factor(paste(...)) returns a factor with length of
nrow(d) and with factor levels determined by the combination of levels
of a and b (provided that you don't get ambiguities as described
Douglas Bates [EMAIL PROTECTED] writes:
Ed L Cashin [EMAIL PROTECTED] writes:
...
I think I need to look at the source of split when I have more time.
?split is probably a better place to start.
I've read that, so I think I need more nitty gritty.
[1] If you really want to be cautious you
Dear all
I would like to represent the outliers in the plot. These few outliers
are much larger than the limit of 50 in the ylim-argument.
plot(daten$month~daten$no,ylim=c(0,50))
I know that it is possible to introduce the information about the
presence of outliers without changing the range
I have the following dataset:
96 plots
12 varieties
2 time points
The experiment is arranged as follows:
A single plot has two varieties tested on it.
With respect to time points, plots come in 3 kinds:
(1) varietyA, timepoint#1 vs. variety B, timepoint#1
(2) varietyA timepoint #2 vs.
On Thu, 8 Apr 2004, Ed L Cashin wrote:
One thing I notice is that I get funny row names using this method. I
can change them back easily enough, though.
I think you can use interaction() rather than the factor(paste()) approach
interaction(d$a, d$b, drop=TRUE) gives a factor with one level
On Thu, 8 Apr 2004, Peter Dalgaard wrote:
Martin Maechler [EMAIL PROTECTED] writes:
Richard a designated sysadmin, which I'm not. Is there any
Richard alternative to downloading a source distribution
Richard and doing the classic ./configure;make;make install
Richard
Hi, Bob.
As I said, I just begin to learn R, so I am not in a proper position to
compare the languages.
But, two things of R immediately captured my attention.
Pro. A lot of standard libraries including a powerful visualization, like
MATLAB.
As a scientific programmer, I really missed this kind
On Thu, 8 Apr 2004, Thomas Lumley wrote:
I believe, though I haven't yet tried, that the 1.9.0 setup will be
better. I think you can put the R Framework in ~/Library/Frameworks and
then put R.app anywhere.
Almost true. You need to edit R.app/Contents/MacOS/R
and change the first line of
hello,
I'm using R 1.8.1 with the lastest snow package on FreeBSD 4.9.
However, when I try to using socket clusters, it's very unstable.
Sometimes it dies half way when I run parSapply(), sometimes
it dies when cluster connection is idle.
I create a socket cluster by following cmd
Dear R-Help,
I am testing a web based erosion model and I would like to get R to submit
the http requests. Now I make the requests manually and use R to extract
the results by scan()ing a url() connection. It would be nice if I could
fully automate the process. Any ideas/examples? The model
If you absolutely have to use the socket approach you will probably
have to track down the particular reason for this failure in your
environment by looking at what is going on under the hood and checking
the pieces. The master sets up a server socket, the slave noted
connect to that socket, and
I don't think the socket stuff is as well tested as PVM and MPI.
We've had a few problems with PVM via rpvm in the past, but I've not
had problems recently (maybe this message will shake them out) and
LAM-MPI via Rmpi seems a bit more stable (though perhaps less
tested).
best,
-tony
Rong-En
From: Shin, Daehyok sdhyok at catchlab.org
Con. Poor design for object-oriented programming, compared with Python.
Have a look at oo.R package if you
want a more traditional oo approach.
__
[EMAIL PROTECTED] mailing list
Hi,
I have run rpart to construct a regression tree. Is there any simple
method to draw a nice picture of it, as it is usually done in books and
paper to visualize the tree?
Thank you,
Christian
***
Christian Hennig
Fachbereich
Shin, Daehyok [EMAIL PROTECTED] writes:
Con. Poor design for object-oriented programming, compared with
Python. I am quite sure that in the aspect of language design
Python is far better than MATLAB and R. Simple, elegant, and
powerful. S3 supports only a basic level of OOP, while S4
See ?post.rpart. That's probably as good as it'll get. You might also want
to look at the maptree package on CRAN.
HTH,
Andy
From: Christian Hennig
Hi,
I have run rpart to construct a regression tree. Is there any simple
method to draw a nice picture of it, as it is usually done in
Christian Hennig wrote:
Hi,
I have run rpart to construct a regression tree. Is there any simple
method to draw a nice picture of it, as it is usually done in books and
paper to visualize the tree?
Thank you,
Christian
***
Gabor Grothendieck ggrothendieck at myway.com writes:
From: Shin, Daehyok sdhyok at catchlab.org
Con. Poor design for object-oriented programming, compared with Python.
Have a look at oo.R package if you
want a more traditional oo approach.
As several people have pointed out, this
Hi,
I like to add class labels on dendograms (as plot(agnesobject)).
Is it possible to be done?
Thank You in advance,
francesca
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
To all,
Does anyone know whether there are any
Discrete Choice Modeling modules for R??
Thanks, Charlie
-
Charles H. Rosa, Ph.D.
Staff Research Engineer
General Motors Corporation
RD and Planning
Mail Code 480-106-359
30500 Mound Rd
Warren, MI 48090-9055
Tel/work/cell:
R Users,
I was just wondering if anyone has written a program (or if there
is a package) out there that calculates the different derivations of AIC
(e.g. AIC, AICc, QAIC, etc.) along with AIC differences (delta's), model
likelihoods, Akaike weights and evidence ratio's (from Burnham and
First, here's the problem I'm working on so you understand the context. I
have a data frame of travel activity characteristics with 70,000+ records.
These activities are identified by unique chain numbers. (Activities are
part of trip chains.) There are 17,500 chains.
I use the chain numbers as
[EMAIL PROTECTED] writes:
What I've found, however, is that it is not easy (or I have not found the
easy way) to split a named vector into a list that retains the vector names.
For example, splitting an unnamed vector (70,000+) based on the chain
numbers takes very little time:
[EMAIL PROTECTED] wrote:
Dear all
I would like to represent the outliers in the plot. These few outliers
are much larger than the limit of 50 in the ylim-argument.
plot(daten$month~daten$no,ylim=c(0,50))
I know that it is possible to introduce the information about the
presence of outliers
Bret Collier wrote:
R Users,
I was just wondering if anyone has written a program (or if
there is a package) out there that calculates the different derivations
of AIC
...
I did not
see anything on the help pages (or in packages, but I could have missed
it)
...
[EMAIL PROTECTED]
Corey Moffet wrote:
Dear R-Help,
I am testing a web based erosion model and I would like to get R to submit
the http requests. Now I make the requests manually and use R to extract
the results by scan()ing a url() connection. It would be nice if I could
fully automate the process. Any
Hi Bret!
There is a step function and an AIC function in the base package. Also
there is a logLik function in base as well.
Finally, there is a stepAIC function in MASS
Hope this helps!
Sincerely,
Erin Hodgess
mailto: [EMAIL PROTECTED]
__
[EMAIL
[I thought this has been covered before, but I can't find it in the
archive...]
PRESS = sum of squares of deletion residuals
Any good book on regression will show you that deletion residual for
multiple regression fitted by least squares is just the raw residual divided
by one minus the hat
Hi,
I have found a strange thing. If I copy and paste the R source file, then
the result is okay with. However, if I use R GUI (File icon and Source
File),
then the result is different. Please let me know possible reasons.
Many thanks in advance.
Best
Ja-Yong
[[alternative
On Fri, Apr 09, 2004 at 01:33:18PM +0900, Ja-Yong Koo wrote:
I have found a strange thing. If I copy and paste the R source file, then
the result is okay with. However, if I use R GUI (File icon and Source
File),
then the result is different. Please let me know possible reasons.
Next time
44 matches
Mail list logo