similar to: RV: problems checking a package

Dear All, We are happy to announce, after a long gestation, the release of the new version of FRAILTYPACK (version 2.2-9.5) which is now available from CRAN. The package fit general frailty models using penalized likelihood estimation, for clustered or recurrent events. For instance : -- ADDITIVE FRAILTY MODELS for proportional hazards models with two correlated random effects (intercept

Dear All, We are happy to announce, after a long gestation, the release of the new version of FRAILTYPACK (version 2.2-9.5) which is now available from CRAN. The package fit general frailty models using penalized likelihood estimation, for clustered or recurrent events. For instance : -- ADDITIVE FRAILTY MODELS for proportional hazards models with two correlated random effects (intercept

Hi, when i run the following code i get this massege: "The instruction at 0x######## reference memory at 0x#######, the memory cannot be "read". and then i have to close R. what is the problem and how can i solve it? thanks in advance Avi my code # frailtypack library(frailtypack) cgd.ag <- read.csv("C:/rfiles/RE/cgd.csv") cgd.nfm <-frailtyPenal(Surv(TStart,

Dear, I want to analyze two-level survival data using a shared frailty model, for which I want to use the R package 'Frailtypack", proposed by Rondeau et al. The dataset was built using SAS software. I also tried to change the format using SPSS and Excell. My (reduced) dataset has following column names: ID entry time status family var1 I used following command: >

Hi, I am trying to convert the following NLMIXED code to NLME, but am running into problems concerning 'Singularity in backsolve'. As I am new to R/S-Plus, I thought I may be missing something in the NLME code. NLMIXED *********** proc nlmixed data=kidney.kidney; parms delta=0.03 gamma=1.1 b1=-0.003 b2=-1.2 b3=0.09 b4=0.35 b5=-1.43 varu=0.5; eta=b1*age+b2*sex+b3*gn+b4*an+b5*pkn+u;

Hi everybody, I'm trying to create a numerical data frame on which to perform PRCC. So far I have created a data frame that consists of function/vector output that displays in numerical form, but when I try and run PRCC (from epiR package) I get the following error message: "Error in solve.default(C) : Lapack routine dgesv: system is exactly singular" It appears this is because

Dear all, Does the frailty.object$history[[1]]$theta returns the Variance of random effect? Why is the value different? Here is an example with kidney data: > library(survival) > data(kidney) > frailty.object<-coxph(Surv(time, status)~ age + sex + disease + frailty(id), kidney) > frailty.object Call: coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id), data

Hi all, I'm new to R and have the following problem: I have a 2 factor design (a has 2 levels, b has 3 levels). I have an object kidney.aov which is an aov(y ~ a*b), and when I ask for model.tables(kidney.avo, se=T) I get the following message along with the table of effects: Design is unbalanced - use se.contrast() for se's but the design is NOT unbalanced... each fator level

Hi everyone, I emailed yesterday thinking I had a problem with reading matrix information as numerical output into epiR. After working on the data frame today, I realise it's a problem with the data itself rather than the syntax. An original data frame of min max T1 1.500000e+01 3.999954e+01 SE1 0.000000e+00 1.000000e+00 PRE

Hi everyone I am attempting to estimate a model with a frailty effect distributed as a lognormal variable.I am using the following code: frailtyPenal(formula, data, ..., RandDist = "LogN") I get the following error message: Error in frailtyPenal(Surv(,) ~ + , : unused argument(s) (RandDist = "") What can I do? Thanks for the help [[alternative HTML version deleted]]

Dear R users - I have a list of patient identifiers and diagnoses from inpatient admissions. I would like to reorganize the list, presently in a long format to a wide format in reshape, but in the absence of a "time" element, I am uncertain how to do this - any help greatly appreciated. ID Dx A nausea A diabetes A kidney failure A heart attack A fever B fever B

Hi, I am doing an analysis to see if these is tissue specific effects on the gene expression data . Our data were collected from 6 different labs (batch effects). lab 1 has tissue type 1 and tissue type 2, lab 2 has tissue 3, 4,5,6. The other labs has one tissue type each. The 'sample' data is as below:

Dear R-users, I have noticed small discrepencies in the reported estimate of the variance of the frailty by the print method for survreg() and the 'theta' component included in the object fit: # Examples in R-2.6.2 for Windows library(survival) # version 2.34-1 (2008-03-31) # discrepancy fit1 <- survreg(Surv(time, status) ~ rx + frailty(litter), rats) fit1 fit1$history[[1]]$theta

Dear friends - We have 25 rats, 14 of these subjected to partial removal of kidney tissue, 11 to sham operation, and then followed for 6 weeks. So far we have data on 26 urine metabolites measured by NMR 7 times during the observation. I have smoothed the measurements by b.splines in fda including a roughness penalty, and inspecting the mean curves for nephrectomized and sham animals indicate

Dear all, in survival5, kidney data set, appears in help page: "survival5 does not reproduce the original analysis." What does it means? thanks in advance TCM -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20010927/e006b08d/attachment.html

Hello R users Would somebody know how to estimate survival from a frailty cox model, using the function survfit and the argument newdata ? (or from any other way that could provide individual expected survival with standard error); Is the problem related to how the random term is included in newdata ? kfitm1 <- coxph(Surv(time,status) ~ age + sex + disease + frailty(id,

  To   R group Help Desk     I am a user of R software.  I am facing a problem while using "frailtyPenal" command in R.2.11.1. When I use these command, R closes completely without any prior alert message. Can I know what would be the reason? My data size is 7050 records with atleast 25 variables.   Looking forward to hear from you.   Truely Kalaivani M.Kalaivani, M.Sc. Scientist-I

I have data I'd like to plot using the formula interface to boxplot. I call boxplot like so: with(mydata, boxplot(count ~ geno * tissue)) I get a boxplot with x axis labels like "wt.kidney". I would like to change the '.' to a newline. Where is this separator configured? Thanks, -Ed

Hello Everyone, I am having problem in defining specific axis for plotting a vactor. vecAVG <- c(0.2, 0.4, 0.6, 0.2, 0.4) names(vecAVG)<-c("brain","heart","kidney","lung","blood") par(mar=c(12,4.1,4.1, 2.1))

Hello listers, Does anybody know how to add text to an xyplot without whipping out the existing curve? That's all. For instance, Lets say you generate a graph like this A <- data.frame(x = rnorm(100), y = rnorm(100)) xyplot(y ~ x, data = A) How would you add 'Hello world'? I tried 6.02E23 different partial solutions found on the web and failed. I just need one EXAMPLE that