search for: babyak

.... The ticks show up fine, but I'd like to alter the color. I know the ticks are produced by scat1d, but after spending a fair bit of time going through documentation, it still isn't clear to me how to do this in the context of lattice. Guidance would be greatly appreciated. Thanks, Mike Babyak Duke University Medical Center My code using R 2.10.1/windows XP myx<-c(1,2,3,4) myy<-c(1,2,3,5) library(rms) d<-datadist(myx) options(datadist="d") myfit<-ols(myy~myx,x=T,y=T) p1<-Predict(myfit,myx =.) library(lattice) #change line to black line <- trellis.par.g...

...I run the code below using something like plot(x~y) and plot(a~b) instead of the two rms summary objects, it works fine, so I am assuming there is something about rms I am missing. I'd be happy to provide some data if that would help. Any guidance would be greatly appreciated. Thanks, Mike Babyak Duke University Medical Center ######################################################## #summary of dietary variables broken out by Groups A, B, and C sumfc<-summary(Group~fcgrain+fcveg+fcfruit+fcmeat+fcdairy+fcnsl+ fcfat+fcsatfat+fcsweet+fcsod, method='reverse', overall=F, te...

...e larger one? Of secondary importance, anything you can tell me about improving my clumsy little simulation code would be welcome--I have little to no idea what I'm doing there. for example, would there be a way to produce the group variable with the disparate Ns more directly? Thanks, Mike Babyak Behavioral Medicine Research Center Duke University Medical Center #question about jitter/llist in rms #R v 2.14.1 under windows 7 #################################################################### #question about jitter/llist in rms require(rms) #simulate some data n = 5000 age = runif...

...some SAS code for power estimation, but I haven't had any luck contacting the authors. Does anyone have any experience with this issue or pointers to where I might do some reading on this topic? Of course, any thoughts on the approriate model would also be more than welcome. Thanks, Mike Babyak

...pdf. I'm sure I'm missing some fundamental concept here, but I'm afraid I'm not seeing it. I'd appreciate if anyone could point me in the right direction. ( I have no trouble writing my own simple latex code and converting it to pdf using pdftex in miktex). Thanks, Mike Babyak Department of Psychiatry and Behavioral Science Duke University Medical Center ************************** Here's the code: library(Hmisc) library(survival) getHdata(pbc) pbc<-upData(pbc, moveUnits=TRUE, labels=c(stage='Histologic Stage\nLudwig Criteria')) kmsurv <- function(S...

Hi, 1. I want 95% not 99% confidence intervals in my summary.Design plot using the Design package. Putting conf.int=.95 as an argument in plot does not work. The default appears to be .99 not .95 as stated in the package Design manual (p. 164). 2. My sweave chuck is below and my output is attached as well as linked here: http://www.sonoma.edu/users/s/stanny/330A/project/ciplot.pdf 3.