[R] LME & data with complicated random & correlational structures

2005-12-01 Thread Keith Chamberlain
Dear List, This is my first post, and I'm a relatively new R user trying to work out a mixed effects model using lme() with random effects, and a correlation structure, and have looked over the archives, & R help on lme, corClasses, & etc extensively for clues. My programming experience is minimal

[R] extend.series not zero padding

2005-12-05 Thread Keith Chamberlain
Dear List, I was trying to verify that I could use extend.series in the wavelets package and kept getting an error when trying to use method="zero". I'm not seeing where my syntax has gone awry. According to the documentation, [see ?extend.series] " method: A character string indicating which ex

[R] Obtaining the adjusted r-square given the regression

2006-01-12 Thread Keith Chamberlain
Hi people, I want to obtain the adjusted r-square given a set of coefficients (without the intercept), and I don't know if there is a function that does it. Exist Dear Alexandra, Without knowing what routine you were using that returned an R-Square value too you, it is a little d

[R] Relating Spectral Density to Chi-Square distribution

2006-01-24 Thread Keith Chamberlain
Dear list, I had some confusion regarding what function too use in order too relate results from spec.pgram() too a chi-square distribution. The documentation indicates that the PSD estimate can be approximated by a chi-square distribution with 2 degrees of freedom, but I am having trouble figu

[R] spec.pgram() normalized too what?

2006-01-24 Thread Keith Chamberlain
Dear list, What on earth is spec.pgram() normalized too? If you would like to skip my proof as to why it's not normed too the mean squared or sum squared amplitude of the discrete function a[], feel free too skip the rest of the message. If it is, but you know why it's not exact in spec.pgram() wh

Re: [R] R-help Digest, Vol 35, Issue 24

2006-01-24 Thread Keith . Chamberlain
wer your question about the reference distribution. The help page does not say it is `normalized' at all: it says it computes the peridogram, and you seem unaware of the definitions of the latter (and beware, there are more than one). On Tue, 24 Jan 2006, Keith Chamberlain wrote: __

Re: [R] How do I normalize a PSD?

2006-02-01 Thread Keith Chamberlain
Dear Tom, Short answer, if your using spec.pgram(), use the smoothing kernel to get a better estimate at the frequency centered in the bandwidth. If your frequency bin of interest is wider than the bandwidth of the kernel, average across frequencies (I think). The estimate appears to be normalized

[R] LME Correlation Component using spectrum()?

2006-03-06 Thread Keith Chamberlain
Dear List-mates, Is there a corStruct constructor already written to calculate the shape of the spectral density for linear models using Fourier estimates (e.g. terms for a linear model derived from the frequency domain)? I have data with a long memory process but do not want to destroy it by taki

[R] Trellis - setting xlim or ylim by data range in whole column or row

2006-03-09 Thread Keith Chamberlain
Dear List-mates, I have been trying to set up a 4x8 trellis plot (that is, 4 columns, 8 rows), and would like to have the axis limits set based on the data range of rows (for ylim) and columns (for xlim). I've been using the call: foo<-xyplot(y~x|Epoch+Subject, type=c("l","r"), pa

Re: [R] Trellis - setting xlim or ylim by data range in whole column or row

2006-03-09 Thread Keith Chamberlain
To: Keith Chamberlain Cc: r-help@stat.math.ethz.ch Subject: Re: Trellis - setting xlim or ylim by data range in whole column or row On 3/9/06, Keith Chamberlain <[EMAIL PROTECTED]> wrote: > Dear List-mates, > > I have been trying to set up a 4x8 trellis plot (that is, 4 columns, 8 > ro

Re: [R] GAM using R tutorials?

2006-03-16 Thread Keith Chamberlain
Dear Templ (& Others), Thank you for pointing out the contributed papers. It is one of those things that I have been introduced to more than once, but I personally never remembered to check the section out while in the acute phase of "problem-can't do it-what the heck is going on." Dear Michael,

[R] Platform independent dialogs & menus?

2006-03-20 Thread Keith Chamberlain
Dear list mates, Are {utils} dialog box functions, and winMenuAdd... functions used to change (e.g. Console) menus, platform dependent? I'm writing a script loaded with .First that provides first time users in a lab course with the ability to select, load, and change between what I called a 'ses

Re: [R] Platform independent dialogs & menus?

2006-03-20 Thread Keith Chamberlain
ines? Respecftully, KeithC. -Original Message- From: Philippe Grosjean [mailto:[EMAIL PROTECTED] Sent: Monday, March 20, 2006 6:18 AM To: Keith Chamberlain Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Platform independent dialogs & menus? Keith Chamberlain wrote: > Dear list ma

[R] MenuRead() Question

2006-03-25 Thread Keith Chamberlain
Dear List-mates, I'm trying to read a tk window menu from a file using {svWidgets} and 'menus.txt' but am receiving Warnings without seeing the desired consequences of the call. library(svWidgets) tkWinAdd("KSesnMain",title="kLab Session Manager for R", pos="+0+0") MenuRead(file="menus.txt") War

Re: [R] MenuRead() Question

2006-03-26 Thread Keith Chamberlain
string vectors going on that I don't understand well enough yet. Do I need to install PCRE-6.6 for this to start working? Rgds, KeithC. -Original Message- From: Keith Chamberlain [mailto:[EMAIL PROTECTED] Sent: Saturday, March 25, 2006 4:26 PM To: 'r-help@stat.math.ethz.ch&

Re: [R] MenuRead() Question

2006-03-26 Thread Keith Chamberlain
that it is a Tk menu. So, rewrite your menu definition file as: > > $Tk.KSesnMain > |$MenuTest > ||Objects~~ ls() > ||- > ||Path ~~ search() Note also that KSesnMain must point to a valid Tk window previously constructed. Best, Philippe Grosjean Kei

[R] {nuSpectral} irreg sampled ts, phase based periodogram

2006-03-28 Thread Keith Chamberlain
Dear List-mates, I just read Mathias, et al (2004) about their {nuSpectral} package in R for a weighted leased squares method for unevenly sampled time-series. What kind of experiences have users had with the package? R-Site Search came up blank. I've never asked about a package that was downlo

[R] Extracting tolerance in R?

2006-12-13 Thread Keith . Chamberlain
Dear list, How is the tolerance for a model parameter in an lm() call extracted? I did not see a solution in the documentation for lm(), or predict(), nor in the archives using 'tolerance' as the search string. I also checked into the nlme package, though nothing popped out at me. Sincerely, Kei