search for: lizzylaw

Hey, Can anyone answer this question. I am working with really large datasets and most of the programs I have been running take quite some time. I heard that R may be faster in Unix. I sthis true and if so can anyone reccomend which system and requirements may allow things to go faster for? Thanks!! Elizabeth Lawson --------------------------------- [[alternative

...Then you can replace the wavelet coefficients and the scaling coefficients using putC and putD from wavethresh. This will leave you with a wd object that you can then reconstruct using wr from wavethresh. I hope this works and is still useful to you. Elizabeth Lawson Elizabeth Lawson <lizzylaws@yahoo.com> wrote: Date: Wed, 19 Oct 2005 08:04:02 -0700 (PDT) From: Elizabeth Lawson <lizzylaws@yahoo.com> Subject: Re: [R] Wavelet reconstruction To: Amir Safari <amsa36060@yahoo.com> Using wavethresh you can recompose a decompose signal using wr. Here is an example...

Hey, Does anyone know if there is a way to get back from goodfit what it estimated the parameters to be? I used the code fit<-goodfit(round(data$PLX_NRX),type="nbinomial" and got a pretty good fit. I could not however duplicate this good fit with any parameter estimates that I had. Any ideas??? Thanks, Elizabeth Lawson ---------------------------------

Hey eveyone! I hope someone can help wiht this question. I have a matirux of all zeros and ones and I would like to indentify all unique patterns in the rows andthe number of times the pattern occurs. I changed all zeros to NA tried to use prelim.norm to identify all patterns of missing data in the rows. I got the message Warning message: NAs introduced by coercion Any ideas of how

hey, I am running R 2.0.1 on a the terminal of my MAC OS X and I have written some C code that I would like to use dynload to bring into R. When I try to compile the code i get the error R.defines.h no such file or directory. I know that the file is in resources/include but where do I have to move the file to so that the compiler can access it? I tried copying Rdefines.h and R_ext folders

Hey, I don;t know if anyone has come across this error before... I am running R on the terminal of my MAC OS X 10.3.4 and I have written C code and compiled it using R CMD SHLIB mycode.c There were no problems in compiling so I now have mycode.o and mycode.so. I used dyn.load("mycode.so") and again, no problems. But when I try to use the code .C("mycode",x)

Hey everyone, I have read the paper by Ritter and Tanner(1992) on Griddy-Gibbs sampler and I am trying to implement it in R without much luck. I was wondering if anyone had used this or could point me to any example code. Thanks, Liz --------------------------------- [[alternative HTML version deleted]]

Hey, I was using lme and proc mixed in SAS to run a empirical bayesian model. I used the same method for both lme and proc mixed (the default REML). I got very similar, but not identical results. I am just wondering if anyone knows what the differneces may be between proc mixed and lme. Any thoughts would be appreciated! Thanks, Liz

Hi, I am using gamm and I run out of memory and R crashes. I have tried to change the amount of memory in R using memory.limit as so far the largest it has allowed memory.limit(size=3072) still ahs this problem. If I try to use max I get the following error > memory.limit(memory.size(max=TRUE)) Error in memory.size(size) : don't be silly!: your machine has a 4Gb address limit Does

Hi, I am hoping some one can help with this. I am using nlme to fit a random coefficients model. It ran for hours before returning Error: Singularity in backsolve at level 0, block 1 The model is > plavix.nlme<-nlme(PLX_NRX~loglike(PLX_NRX,PD4_42D,GAT_34D,VIS_42D,MSL_42D,SPE_ROL,XM2_DUM,THX_DUM,b0,b1,b2,b3,b4,b5,b6,b7,alpha), + data=data, + fixed=list(b0 +

Hi, I am hoping some one can help me. I am learning to use C and would like to learn how to call c code in R. I have look at Writing R Extensions and I tried to copy the example on page 38 void convolve(double *a, int *na, double *b, int *nb, double *ab) { int i, j, nab = *na + *nb - 1; for(i = 0; i < nab; i++) ab[i] = 0.0; for(i = 0; i < *na; i++) for(j = 0; j <

Hey, I am estiamting a partially linear model y=X\beta+f(\theta) where the f(\theta) is estiamted using wavelets. Has anyone heard of methods to test if the betas are significant or to address model fit? Thanks for any thoughts or comments. Elizabeth Lawson __________________________________________________ [[alternative HTML version deleted]]