Re: [R] Survey imputation

2013-06-13 Thread Scott Raynaud
I paln on using R to do the imputation once I figure out how to do it.   - Original Message - From: Bert Gunter To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Thursday, June 13, 2013 12:45 PM Subject: Re: [R] Survey imputation Is this an R question? Seems like it be

[R] Survey imputation

2013-06-13 Thread Scott Raynaud
I'm working with NHIS survye data.  I'd like the to use muliple imputation to cover the missing data for the variables in which I'm interested.  My question concerns the use of certain variables in the imputation model.  For example, race would be an important predictor in the imputation model,

Re: [R] SPlus script

2013-06-06 Thread Scott Raynaud
for helping to diagnose and fix the problem. - Original Message - From: Ista Zahn To: Scott Raynaud Cc: William Dunlap ; "r-help@r-project.org" Sent: Thursday, June 6, 2013 9:15 AM Subject: Re: [R] SPlus script Presumably something like r <- sshc(50) print(r) But if you wer

Re: [R] SPlus script

2013-06-06 Thread Scott Raynaud
how to call the function appreciated. - Original Message - From: William Dunlap To: Scott Raynaud ; "r-help@r-project.org" Cc: Sent: Wednesday, June 5, 2013 2:17 PM Subject: RE: [R] SPlus script Both the R and S+ versions (which seem to differ only in the use of _ for assignm

Re: [R] SPlus script

2013-06-06 Thread Scott Raynaud
Ok. I tried copying the original program, changing all _ to <- and running in 3.0.1.  Still no error message or output.  It's a mystery to me. - Original Message - From: William Dunlap To: Scott Raynaud ; "r-help@r-project.org" Cc: Sent: Wednesday, June 5, 2013 2

Re: [R] SPlus script

2013-06-05 Thread Scott Raynaud
the historical group.  Maybe I'm forgetting something here...  It's been a while since I ran this.  Thanks for your help. - Original Message - From: William Dunlap To: Scott Raynaud ; "r-help@r-project.org" Cc: Sent: Wednesday, June 5, 2013 2:17 PM Subject: RE: [R] SPlus

Re: [R] SPlus script

2013-06-05 Thread Scott Raynaud
See my comments below. From: Pascal Oettli To: Scott Raynaud Sent: Wednesday, June 5, 2013 10:02 AM Subject: Re: [R] SPlus script Hello, 1) It is always nice to say something as "Hello", Tried, but could make it come out in anything other than fo

[R] SPlus script

2013-06-05 Thread Scott Raynaud
This originally was an SPlus script that I modifeid about a year-and-a-half ago.  It worked perfectly then.  Now I can't get any output despite not receiving an error message.  I'm providing the SPLUS script as a reference.  I'm running R15.2.2.  Any help appreciated.    ***

Re: [R] Error message

2012-08-29 Thread Scott Raynaud
rson so someone will have to tell me the command to obtain the Linux information requested.  Then maybe I can figure whether to go to the Debian board or not. - Original Message - From: Duncan Murdoch To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Tuesday, August 28, 20

Re: [R] Error message

2012-08-28 Thread Scott Raynaud
  - Original Message - From: Duncan Murdoch To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Tuesday, August 28, 2012 10:40 AM Subject: Re: [R] Error message On 28/08/2012 10:31 AM, Scott Raynaud wrote: > I suddenly started getting the error message below. > Not sure why.

[R] Error message

2012-08-28 Thread Scott Raynaud
I suddenly started getting the error message below. Not sure why.  If I type intalled.packages() it shows Matrix and lme4 installed.  Can someone tell what's going on and what I need to do to remedy the problem?  I'm running on a Linux box.             Loading required package: Matrix Loading requ

Re: [R] Apt-get

2012-03-16 Thread Scott Raynaud
Ok.  When I type in sudo add-apt-repository ppa:marutter/rrutter I get a prompt for a password.  I enter my domain password and get some print about adding the ppa and requeting that I press enter.  I press enter and this is followed by print stating that the key B04C661B is being requested from h

Re: [R] Apt-get

2012-03-15 Thread Scott Raynaud
iller To: Scott Raynaud ; "r-help@r-project.org" Cc: Sent: Wednesday, March 7, 2012 12:34 PM Subject: Re: [R] Apt-get Google is really useful for questions like this. http://cran.r-project.org/bin/linux/ubuntu/ -

[R] Apt-get

2012-03-07 Thread Scott Raynaud
I have a box set up with Kubuntu as the OS.  I didn't perform the R install but was told the version of R available via the apt-get command was 2.13.1.  Is there any way to get 2.14.0 in that same manner? __ R-help@r-project.org mailing list https://sta

[R] Access to OpenBLAS

2012-02-14 Thread Scott Raynaud
My IT people have set up a Kubuntu box with an RWkard front end.  I have OpenBLAS set up as a shared BLAS but I'm not sure how to get R to see it.  A.3.1 of the installation docs talks about it but I'm not clear if I need a option on my startup line or if I need to find a config file.  The BLAS is

[R] Simulation confidence interval

2012-02-01 Thread Scott Raynaud
The follwing is a code snippet from a power simulation program that I'm using:   estbeta<-fixef(fitmodel)  sdebeta<-sqrt(diag(vcov(fitmodel)))   for(l in 1:betasize)   {     cibeta<-estbeta[l]-sgnbeta[l]*z1score*sdebeta[l]     if(beta[l]*cibeta>0)  powaprox[[l]]<-powaprox[[l]]+1  

[R] Recompile

2012-01-23 Thread Scott Raynaud
Looked at several posts and the installtions docs and still not clear.  If I compile source codes and then somewhere down the line add a new package, then I have to recompile my entire installation, correct?  Seems like this is the sentiment of the emails I read. ___

[R] GUI

2012-01-18 Thread Scott Raynaud
I'm setting up a Linux box to run R.  I ususally run in a Windows envrionment but after reading the docs I'm not sure what to expect in terms of the front end appearance in Linux.  Does it resemble Windows or will I need Rkward or R Commander? __ R-h

[R] BLAS

2012-01-17 Thread Scott Raynaud
I'm setting up an Ubuntu virtual machine that will use 4-Intel Xeon CPU x5650.  I'd like to compile R with a BLAS but the question is whcih one.  Seems like the only free ones are GotoBLAS which I'm not sure is being maintained for newer CPUs and OpenBLAS for Loongson CPUs.  I saw a favorable re

Re: [R] grplasso

2012-01-13 Thread Scott Raynaud
So does anyone use this package?   - Original Message - From: Scott Raynaud To: "r-help@r-project.org" Cc: Sent: Tuesday, January 10, 2012 1:40 PM Subject: grplasso I want to use the grplasso package on a data set where I want to fit a linear model.  My interest is in i

[R] grplasso

2012-01-10 Thread Scott Raynaud
I want to use the grplasso package on a data set where I want to fit a linear model.  My interest is in identifying significant beta coefficients.  The documentation is a bit cryptic so I'd appreciate some help.   I know this is a strategy for large numbers of variables but consider a simple cas

Re: [R] rbinom

2011-12-30 Thread Scott Raynaud
but I'd need to square my x[,3] values before multiplying them by beta.  Can I say: x[,3]<-(rnorm(length,meanpred[3],sqrt(varpred[3])))^2 in lieu of x[,3]<-rnorm(length,meanpred[3],sqrt(varpred[3]))? - Original Message - From: peter dalgaard To: Scott Raynaud Cc: "r-hel

[R] rbinom

2011-12-27 Thread Scott Raynaud
I have the following code (which I did not write) that generates data based on a logistic model.  I'm only getting a single record with y=1.  It seems implausible that in 50k cases that have a single y=1.  Does that ring alarm bells for anyone else?   beta<-c(-1.585600,-0.246900) betasize<-lengt

[R] simulation

2011-12-16 Thread Scott Raynaud
I'm using an R program (which I did not write) to simulate multilevel data (subjects in locations) used in power calculations. It uses lmer to fit a mixed logistic model to the simulated data based on inputs of means, variances, slopes and proportions:   (fitmodel <- lmer(modelformula,data,famil

Re: [R] R/parallel

2011-12-08 Thread Scott Raynaud
recompiling with BLAS and pnmath would be a better option.   My main question is how to handle the random number generation when the child processes are spawned.  That's a problem no matter what method I choose to create the threads.   From: Tal Galili To: Scott Raynaud Cc: &quo

[R] R/parallel

2011-12-08 Thread Scott Raynaud
 I want to take advantage of my multicore CPU to speed up a loop in a simulation program.  I didn’t write the code, but the iterations appear independent to me, at least in the sense that the results of one loop do not depend on previous ones.  Right now I’m relegated to a Windows box that runs

Re: [R] modelling and R misconceptions; was: package installtion

2011-11-17 Thread Scott Raynaud
My responses are in brackets below, plus a final note after the main text. - Original Message - From: Uwe Ligges To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Thursday, November 17, 2011 9:16 AM Subject: Re: [R] modelling and R misconceptions; was: package installtio

Re: [R] package installtion

2011-11-17 Thread Scott Raynaud
See my responses in brackets below.   - Original Message - From: Rolf Turner To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Wednesday, November 16, 2011 6:04 PM Subject: Re: [R] package installtion On 17/11/11 05:37, Scott Raynaud wrote: > That might be an option if

[R] package installtion

2011-11-17 Thread Scott Raynaud
(fitmodel <- lmer(modelformula,data,family=binomial(link=logit),nAGQ=1)) mtrace(fitmodel) I added the mtrace to catch the error, but get the following: Error in mtrace(fitmodel) : Can't find fitmodel How can I debug this?   - Original Message - From: Rolf Turner To: Scott Rayn

[R] package installtion

2011-11-16 Thread Scott Raynaud
MLWin sets the associated fixed effects to 0.  When R choked, I increased from 20 to 60 as my minimum as suggested in the MLPowSim documentation.  Still no luck.   - Original Message - From: Uwe Ligges To: r-help@r-project.org Cc: Scott Raynaud Sent: Wednesday, November 16, 2011 1:01

[R] package installtion

2011-11-16 Thread Scott Raynaud
I like R since it's free, but I can't work around the problem I'm currently having.   - Original Message - From: Uwe Ligges To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Wednesday, November 16, 2011 9:48 AM Subject: Re: [R] package installtion On 16.11

[R] package installtion

2011-11-16 Thread Scott Raynaud
nitial par (65) The first 25 are like 26.  So, it seems I'm having the same problem as before.  Again, if this is due to a column of zeroes in my x matrix, the best solution would be to assign zeroes to the fixed effects, but I'm not sure if there's a way to do this.   - Forwarde

Re: [R] package installtion

2011-11-16 Thread Scott Raynaud
an be done.  I've looked at the documentation but it's still not clear.   - Original Message - From: Uwe Ligges To: Scott Raynaud Cc: "r-help@r-project.org" Sent: Wednesday, November 16, 2011 2:44 AM Subject: Re: [R] package installtion On 15.11.2011 21:34, Scott Rayn

Re: [R] package installtion

2011-11-15 Thread Scott Raynaud
gn matrix, but I'm not sure.  Any other possibilities here and how can I diagnose? ----- Original Message - From: Scott Raynaud To: "r-help@r-project.org" Cc: Sent: Tuesday, November 15, 2011 2:11 PM Subject: Re: package installtion Never mind-I fixed it. My script is throwing the fo

Re: [R] package installtion

2011-11-15 Thread Scott Raynaud
no longer available on lme4 but I have AGQ specified.    Here's the line that fits my model: (fitmodel <- lmer(modelformula,data,family=binomial(link=logit),method="AGQ")) If I change it to nAGQ I still get an error. Any ideas as to what's going on?   ----- Original Mess

[R] package installtion

2011-11-15 Thread Scott Raynaud
I'm getting the following error in a script: "Error: could not find function "lmer."    I'm wondering of my lme4 package is installed incorrectly.  Can someone tell me the installation procedure?  I looked at the support docs but couldn't translate that into anything that would work. __

[R] mlpowsim

2011-11-15 Thread Scott Raynaud
I'm using Bill Browne's MLPowSim to do some sample size estimation for a multilevel model.  It creates an R program to carry out the estimation using lmer in the lme4 library.  When there are predictors with more than two categories one has to modify the code generated to account for the multino

[R] Random multinomial variable

2011-11-04 Thread Scott Raynaud
I need some help interpreting the following code which is part of a mutlilevel model simulation with 2 levels.  I've put in comments with my understanding of the code, but I'm not sure how [i2id] is functioning.  It's defined in another part of the program as l2id<-rep(c(1:n2),each=n1) which loo

[R] Power mixed model ordinal logistic regression

2011-10-27 Thread Scott Raynaud
Is there a package that will perform power calculations for mixed model ordinal logistic regression?  I searched and came up with nothing. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Power mixed model ordinal logistic regression

2011-10-26 Thread Scott Raynaud
Is there a package that will perform power calculations for mixed model ordinal logistic regression?  I searched an came up with nothing. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

Re: [R] SPlus to R

2011-10-06 Thread Scott Raynaud
istics From: William Dunlap @r-project.org> Sent: Wednesday, October 5, 2011 11:25 AM Subject: RE: [R] SPlus to R I think you only have to change the multi-argument returns to call list.  You can remove the name from the single argument return, as it will be ign

Re: [R] SPlus to R

2011-10-06 Thread Scott Raynaud
I did some testing and I believe the program is operating properly.  It takes some time to finish, especially as sample sizes get larger, but I seem to be able to reproduce the results from the original paper.  Right now I'm most interested in method 3.  I set nc=40 and d=.2 as in the paper.  Th

Re: [R] SPlus to R

2011-10-05 Thread Scott Raynaud
It seems I have things set up correctly.  I suspect that the arguments sshc(100,10) are the isuue.  It seems that the 100,10 is not necessary since the code itself specifies the arguments.  It runs and produces a power curve if I simply type sshc() but it also seems to try to keep running someth

Re: [R] SPlus to R

2011-10-05 Thread Scott Raynaud
igges Cc: Barry Rowlingson ; "r-help@r-project.org" Sent: Wednesday, October 5, 2011 8:03 AM Subject: Re: [R] SPlus to R On 05.10.2011 13:44, Scott Raynaud wrote: > Hope I did this right.  I repeated what I'd done before: > > 1) Opened script > 2) Selected run all (

Re: [R] SPlus to R

2011-10-05 Thread Scott Raynaud
Hope I did this right.  I repeated what I'd done before:   1) Opened script 2) Selected run all (this produced my inital post   Then as suggested I:   3) Typed ls() 4) Saw that the function was present and issued sshc(100,10)   Here's what I got:   > ls() [1] "c.searchd" "convex"    "Epower"    "ne

[R] SPlus to R

2011-10-04 Thread Scott Raynaud
I'm trying to convert an S-Plus program to R.  Since I'm a SAS programmer I'm not facile is either S-Plus or R, so I need some help.  All I did was convert the underscores in S-Plus to the assignment operator <-.  Here are the first few lines of the S-Plus file:   sshc _ function(rc, nc, d, meth