search for: medstat

This message was also sent to the MEDSTATS mailing list, so here is the reply I posted to that: Philippe, The machinery to use is to split follow-up time so finely that you can safely assume that rates are constant in each interval, and then just stuff it all into a Poisson model. This allows you to use any kind of time-dependent variabl...

Hi, Given that the R forum is mostly a site to discuss issues related to the language, are there any forums that discuss statistical issues? I am specifically interested in the statistical issues in experimental design. TIA, Anjan -- ============================= anjan purkayastha, phd bioinformatics analyst whitehead institute for biomedical research nine cambridge center cambridge, ma 02142

Please forgive the cross posting. I sent this to Medstats and have not received any response. I am looking for a good written explanation of components of variance (a.k.a variance components), in particular as they are used in genetic analyses. Does anyone have a suggestion? I am interested in books, articles and on-line material. I intent to give the m...

...) missed_meds <- c(TRUE, FALSE, FALSE, FALSE, FALSE, FALSE, TRUE) no_exercise <- c(FALSE, FALSE, TRUE, FALSE, TRUE, TRUE, TRUE) b2 <- data.frame(Date, Time, bg, missed_meds, no_exercise) b2$Date <- as.Date(b2$Date) # add "status" columns, could also be defined as factor. b2$medStat <- c("missed_meds",NA, NA, NA, NA, NA, "missed_meds") b2$exercise <- c(NA, NA, "missed_exercise",NA,"missed_exercise", "missed_exercise", "missed_exercise") Then your ggplot call would be like this: ggplot(data = b2, aes(x = Date,...

...FALSE, FALSE, FALSE, FALSE, TRUE) > no_exercise <- c(FALSE, FALSE, TRUE, FALSE, TRUE, TRUE, TRUE) > > b2 <- data.frame(Date, Time, bg, missed_meds, no_exercise) > > b2$Date <- as.Date(b2$Date) > # add "status" columns, could also be defined as factor. > b2$medStat <- c("missed_meds",NA, NA, NA, NA, NA, "missed_meds") > b2$exercise <- c(NA, NA, "missed_exercise",NA,"missed_exercise", > "missed_exercise", "missed_exercise") > > Then your ggplot call would be like this: > > gg...

Hello evyrone, I would be grateful if you could help me in (I hope) simple problem. I fit a gam model (from mgcv package) with several smooth functions . I don't know how to extract values of just one smooth function. Can you please help me in this? Kind regards, Daniel Rabczenko

Hello, I would be grateful if anybody can help me in finding an R function to compute weighted Spearman correlation coefficient? Kind regards, Daniel

Dear UseRs: Recently, I have read an article regarding the association between age and lymph node metastases. http://jco.ascopubs.org/content/27/18/2931.long In statistical analysis, the authors stated "Because a nonlinear relationship between age and lymph node involvement was expected based on existing literature, lymph node involvement was also regressed on age using nonparametric

Hi everyone, I'm a bit new here (and new to R), and I was trying to do an OLR, and testing the parallel slope assumption seems be very important. I browsed through past postings, and didn't find much to help me in this area. I was wondering if anyone knew how I could go about doing this. Thank you. -- View this message in context:

Hi, Suppose a binomial GLM with both continuous as well as categorical predictors (sometimes referred to as GLM-ANCOVA, if I remember correctly). For the categorical predictors = indicator variables, is then there a suggested minimum frequency of each level ? Would such a rule/ recommendation be dependent on the y-side too ? Example: N is quite large, a bit > 100. Observed however are

...5 5 [4,] 5 5 5 5 Any ideas would be much appreciated. Best regards, Hein Hein Putter Department of Medical Statistics Leiden University Medical Center P.O.Box 9604, 2300 RC Leiden, The Netherlands phone +31-71-5276827 fax +31-71-5276799 e-mail h.putter at lumc.nl http://www.medstat.medfac.leidenuniv.nl/MS/Hp -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !...

Hi All, I am trying to import a .sas7bdat file into R. I tried using Hmisc package's sasxport.get() function. temp <- sasxport.get("path\abcd.sas7bdat") I get an error that says "Error in lookup.xport(file) : file not in SAS transfer format" I am not familiar with SAS transfer format. Could somebody please clarify what is it that I am missing, Thanks for your