[R] command for stopping and then hit return and continue (fwd)

2006-11-09 Thread Mark . Leeds
-- 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 that

[R] as.zoo behavior (

2006-11-03 Thread Mark . Leeds
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

Re: [R] did my searching but still couldn't find anything for bayesian dlm

2006-11-01 Thread Mark . Leeds
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

[R] strange algorithm needed

2006-11-01 Thread Mark . Leeds
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

Re: [R] strange algorithm needed

2006-11-01 Thread Mark . Leeds
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 need a strange

[R] ncp in dt

2006-11-01 Thread Mark . Leeds
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. __

Re: [R] ncp in dt

2006-11-01 Thread Mark . Leeds
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 to mu

Re: [R] estimating state space with exogenous input in measurement eq.

2006-09-11 Thread MARK LEEDS
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

Re: [R] pairwise.t.test vs. t. test

2006-09-07 Thread MARK LEEDS
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.

Re: [R] pairwise.t.test vs. t. test

2006-09-07 Thread MARK LEEDS
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,Qinghong,ST.LOUIS,Molecular Biology [EMAIL PROTECTED]; r-help

Re: [R] alpha in Holt-Winters method

2006-08-31 Thread MARK LEEDS
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

Re: [R] cumulative growth rates indexed to a common starting point over n series of observations

2006-08-31 Thread MARK LEEDS
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: r-help@stat.math.ethz.ch Sent: Thursday, August 31, 2006 4:22 PM Subject: [R] cumulative growth rates indexed

Re: [R] for() loop question

2006-08-26 Thread MARK LEEDS
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: r-help@stat.math.ethz.ch Sent: Saturday, August 26, 2006 1:06 PM Subject: [R] for() loop question Dear Lister, If I

Re: [R] fitting a gaussian to some x,y data

2006-08-25 Thread MARK LEEDS
to this with the purpose of determining systematically the center of the normal-like shape of the spatial feature. I have successfully done so in R with a polynomial but I 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

Re: [R] extremely slow recursion in R?

2006-08-24 Thread MARK LEEDS
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

Re: [R] aggregate example : where is the state.region variable?

2006-08-21 Thread MARK LEEDS
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.

Re: [R] E(1/(X+c)) ? where X ~ Chi-square(n) and c is a constant.

2006-05-15 Thread MARK LEEDS
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

[R] simple .bat question but i can't find out where to post. my apologies for posting a non R question

2006-05-14 Thread MARK LEEDS
- 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 out where

[R] two relatively simple questions

2006-05-12 Thread MARK LEEDS
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

[R] .Rprofile file and .First function

2006-05-12 Thread MARK LEEDS
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

[R] I think simple R question

2006-01-12 Thread Mark Leeds
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

[R] I think not so hard question

2006-01-12 Thread Mark Leeds
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

[R] a series of 1's and -1's

2006-01-11 Thread Mark Leeds
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 )

[R] complex matrix manipulation question

2006-01-10 Thread Mark Leeds
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

[R] R newbie example code question

2006-01-09 Thread Mark Leeds
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

Re: [R] R newbie example code question

2006-01-09 Thread Mark Leeds
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

[R] brown, durbin , evans ( 1975 )

2006-01-09 Thread Mark Leeds
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

Re: [R] .Rprofile files (was R newbie configuration)

2006-01-05 Thread Mark Leeds
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), function(...) grDevices::ps.options(horizontal

[R] R newbie configuration

2006-01-04 Thread Mark Leeds
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

[R] newbie R question

2006-01-03 Thread Mark Leeds
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

Re: [R] newbie R question

2006-01-03 Thread Mark Leeds
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] newbie R question

[R] newbie where to look question

2006-01-03 Thread Mark Leeds
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

[R] ESS and Emacs

2005-12-29 Thread Mark Leeds
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 more

[R] FW: R and read.irts

2005-12-28 Thread Mark Leeds
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

Re: [R] R and read.irts

2005-12-28 Thread Mark Leeds
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' Subject: FW: R