-- Forwarded message --
Date: Thu, 9 Nov 2006 16:38:46 -0500
From: "Leeds, Mark (IED)" <[EMAIL PROTECTED]>
To: "Leeds, Mark (IED)" <[EMAIL PROTECTED]>
Subject: command for stopping and then hit return and continue
when i do sequential plotting, i use the park(ask=TRUE) command so
hi all : the code pasted below runs but then, a dput on
rollmeandifflogbidask gives me what is below the code. the structure of
rollmeandifflogbidask is a zoo object but with a "frequency"
so it's not the same structure as the original actual diff and this
really causes things to blow up in later
i meant dt function in my email below. my aopolgies. On Wed, 1 Nov 2006
[EMAIL PROTECTED] wrote:
>
> i just wanted to know if the "non centraility parameter" as it is referred
> to in the dt test is
> what i would think of as mu zero if i was doing a test of a
> hypothesis that the mean was equal
i just wanted to know if the "non centraility parameter" as it is referred
to in the dt test is
what i would think of as mu zero if i was doing a test of a
hypothesis that the mean was equal to mu zero and the variance was
unknown. thanks.
__
R-help@sta
thanks to all for your replies. very neat.
On Wed, 1 Nov 2006, Gabor
Grothendieck wrote:
> Try:
>
> rep(1:3, length = 17)
>
> or
>
> as.numeric(gl(3, 1, 17))
>
> or
>
> y <- 1:17
> replace(y, TRUE, 1:3) # ignore warning
>
>
> On 11/1/06, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> >
> >
> > i
i need a strange algorithm that i can easily do in a loop but need
o do without looping.
suppose i have a vector of length y filled with zeros and a number x.
then, i want a new vector which is (1,..x,1..x,1.x,1.y mod
x )
so if y was of length 17 and x was 3, the resultant vect
thanks gavin : i'll check those out.
On Wed, 1 Nov 2006, Gavin Simpson
wrote:
> On Tue, 2006-10-31 at 19:12 -0500, Leeds, Mark (IED) wrote:
> > I familarized myelf with kalmanlike and structts which are approaches
> > for building and estimating ( and forecasting ) state space models ( or
> > the
below is very close to a standard kalman filter setup except for the
exogenous u[k] so i would check on www.r-project.org
for any packages that do kalman filtering. ( also, good reference for state
space is a book by durbin and koopman but i forget the title exactly ).
- Original Message
thanks. i assumed we we were talking about the standard textbook difference
between the t test and pairwise t test.
my bad.
- Original Message -
From: "Chuck Cleland" <[EMAIL PROTECTED]>
To: "MARK LEEDS" <[EMAIL PROTECTED]>
Cc: "Li,Qinghon
no, because the formula for the test statistics ( even assuming that
variances are equal ) of the two different tests are different. in the
pairwise t test, the pairwise differences are
viewed as one sample so it turns into a one sample test. any intro stat book
will have the formulas.
i think you can add 1 to whatever series you are considering and then use
cumprod() but check that that works.
also do ?cumprod
- Original Message -
From: <[EMAIL PROTECTED]>
To:
Sent: Thursday, August 31, 2006 4:22 PM
Subject: [R] cumulative growth rates indexed to a common starting p
you would have to take your series, y_1, ... y_t and run it over various
values of alpha ( say zero to 1 insteps of .01 ) and see which one gives say
the least MSE.
but that will be the optimal alpha for the window you are looking at. it
doesn't mean it will be the optimal apha for the data you
let us know what you want to do because the beauty of R is that, in many
cases, you may not have to loop.
- Original Message -
From: "Wensui Liu" <[EMAIL PROTECTED]>
To:
Sent: Saturday, August 26, 2006 1:06 PM
Subject: [R] for() loop question
> Dear Lister,
>
> If I have a list of nu
> can't figure out how to do it with a gaussian.
>
> Does that make sense?
>
> Thanks!
> Michael
>
>
> On Aug 25, 2006, at 2:04 PM, MARK LEEDS wrote:
>
>> hi : i'm not clear on what you mean but someone else might be ? if you
>> say ( x,y), then it
i'm sure someone else will explain the recursion issue but , as far as your
program running a few days, you don't have to wait this long. if you are in
windows and do
a ctrl alt delete and then click on processes, if the memory usage being
used by that R process is staying EXACTLY the same and
these people/experts provide all these packages and documentation as a FAVOR
and for the fact that they enjoy spreading knowledge/statistical computing
abilities etc. It's not their job so I think criticism of the docs and the
fact that they use a variable from another place is kind of harsh.
i don't know it off the top of my head but there's a relation between powers
of the chi squared and the gamma dsitribution.
check it out in casella and berger or arnold.
mark
- Original Mes
- Original Message -
From: MARK LEEDS
To: [EMAIL PROTECTED]
Sent: Sunday, May 14, 2006 9:49 PM
Subject: simple .bat question but i can't find out where to post. my
apologiesfor posting a non R question
I have what I think is a simple .bat question but i can't figure ou
I am new to R ( and windows XP ) and am in the process of setting it up I put
an options (repos=http://.fhcrc.org";) file in my .Rprofile file
and it works but that's all I put in and i don'tknow what to put in my .First
function. I understand that .Rprofile is used for more general options
and
I am using R on windows XP and I really tried to find the answers to below but
I haven't been succesful.
the two questions are :
1) is there way of seeing what packages you have loaded ( search() gives you
the automatically included ones ).
2) is there a way of uninstalling packages that you h
I'm sorry to bother this list so much
But I haven't programmed in
A while and I'm struggling.
I have a vector in R of 1's and -1's
And I want to use a streak of size Y
To predict that the same value will
Be next.
So, suppose Y = 3. Then, if there is a streak of three
ones in a row, then I wil
I have a vector x with #'s ( 1 or -1 in them ) in it and I want to
"mark" a new vector with the sign of the value of the a streak
of H where H = some number ( at the next spot in the vector )
So, say H was equal to 3 and
I had a vector of
[1] [2] [3] [4] [5] [6] [7] [8] [9] [10]
1
Does anyone know of a simple test
in any R package that given
a series of negative ones and positive
ones ( no other values are possible in the series )
returns a test of whether the series is random or not.
( a test at each point would be good but
I can use the apply function to implement
that ) ?
Does anyone know of a newsgroup where
I can ask a question about formulating
a problem statistically ? I don't
think this is where this is done.
It has to do ( I think ) with logistic
regression.
Thanks
I've done stuff like this before but
it's been a while and I'm stuck.
Suppose I have a matrix with one
column x and another column y
and both are numeric and let the
row index of the matrix be i
Starting at index i ( i would equal on the first iteration )
when the cumulative sum of x_i+1 - x_i
Does anyone know where
I can get R code for plotting
the Brown , Durbin
and Evans cumsum
procedure ( 1975 ) ?
I wrote my own code but
I am a little worried
that my confiodence bands
may not be correct ( I find the formula
in the original paper confusing and S+Finmetrics
has a formula but that fo
wow. example() was exactly what
I wanted. you all have made all these
incredible facilities
for a new person. thanks.
mark
-Original Message-
From: Uwe Ligges [mailto:[EMAIL PROTECTED]
Sent: Monday, January 09, 2006 1:30 PM
To: Mark Leeds
Cc: R-Stat Help
Subject: Re
Sometimes I print out a package
and read about it and there
are sometimes nice examples
that I would like to run myself.
Is there a way to bring them
into R from the package or
are they only meant to be typed
in manually ? If manual is the
only way, that's fine. I was
just checking whether there
I just had a question for anyone who has done this.
I am currently using R in a research environment
but I was wondering if anyone has had experience using
it in a production/real time environment where data
is coming in quite quickly ( every second or so ) and
things have to be calculated quickly
y, January 05, 2006 4:04 AM
To: Petr Pikal
Cc: Mark Leeds; R-Stat Help
Subject: Re: [R] .Rprofile files (was R newbie configuration)
And here is one with a working setHook call.
options(show.signif.stars=FALSE)
setHook(packageEvent("grDevices", "onLoad"),
funct
I think I did enough reading on my
Own about startup ( part of the morning
And most of this afternoon )
to not feel uncomfortable asking
for confirmation of my understanding of this startup stuff.
Obviously, the startup process is more complicated
Than below but, for my R newbie purposes,
It seem
I don't want to bother anyone
with specific questions
because I am a R newbie and
I see that there is TON ( emphasis on
Ton ) of documentation
out there but could
someone just tell me the
best placed to look/read
for learning about ( for R-2-2.1 in Windows )
.Rprofile
.REnviron.
.Rdata
.First f
t has been done
in R by everyone. Truly
amazing. You should
all be quite proud
about what you have created.
Mark
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Tuesday, January 03, 2006 4:16 PM
To: Mark Leeds
Cc: R-Stat Help
Subject: Re: [R] ne
I'm sorry to bother everyone with a stupid
question but, when I am at an R prompt in Windows,
is there a way to see what packages
you already have installed from the R site so
that you can just do library(name_of_package)
and it will work.
I've looked at help etc but I can't find
a command like t
I have been using the document written by John Fox titled
Sn Introduction to ESS + XEmacs for Windows
Users of R.
It's a very nice document and
I went through it carefully but I got
an error when I finished it and launched XEmacs.
The error is "cannot open load file : ess-site".
So, I did mor
.10,3
09:52:19.8536,67.05,-1
09:52:48.1685,67.09,1
09:53:45.2387,67.08,-1
09:53:51.6611,67.08,-2
09:53:54.0639,67.09,3
09:55:07.0671,66.99,5
09:55:12.5074,66.90,-7
-Original Message-
From: Mark Leeds
Sent: Wednesday, December 28, 2005 6:42 PM
To: 'r-help@stat.math.ethz.ch'
Sub
I have never worked with R before so I am
sorry if this is a bad question but I've
tried and tried ( all day ) and I can't figure this
problem out. I have the code below and I included the
data file as an attachment.
The code works in term of reading in the data
correctly but when the graph gets
37 matches
Mail list logo