similar to: What library is coxreg?

Dear List, After including cluster() option the coxreg (from eha package) produces results slightly different than that of coxph (from survival) in the following time-dependent treatment effect calculation (example is used just to make the point). Will appreciate any explaination / comment. cheers, Ehsan ############################ require(survival) require(eha) data(heart) # create weights

Dear mailing list, Do you know a specific package in R where I can create an artificial survival data with controlled censoring condition? (For example, I want to simulate a data set with 5 variables and 20% of the observations are censored.) Many thanks. Jei

Dear R-users, I need to use the aftreg function in package 'eha' to estimate failure times for left truncated survival data. Apparently, survreg still cannot fit such models. Both functions should be fitting the accelerated failure time (Weibull) model. However, as G?ran Brostr?m points out in the help file for aftreg, the parameterisation is different giving rise to different

Hi out there Is there a way to get the estimated coefficients in a logistic / Cox regression without having to specify a 'formula' but by only giving the model matrix? Example for Cox regression: ## predictors n <- 50 q1 <- rnorm(n) q2 <- rgamma(n, 2, 2) Z <- cbind(q1, q2) ## response ttf <- rexp(n) tf <- round(runif(n)) ## compute estimates res <- coxph(Surv(ttf,

Hi. I am using Windows version of R 1.8.1. Being somewhat new to survival analysis, I am trying to compare cph (Design) with coxph (survival) for use with a survival data set. I was wondering why cph and coxph provide me with different confidence intervals for the hazard ratios for one of the variables. I was wondering if I am doing something wrong? Or if the two functions are calculating hazard

Dear all, How can I analyze a life table (e.g. for a cohort of insects) in R? I have 20 insects in 200 cages with two different treatments, whose survival is followed over time, such that, e.g., in one treatment, the number of animals surviving is c(20,18,16,12,10,8,4,0), while in the other treatment the survival is c(20,20,18,18,16,15,15,14) at 8 subsequent time intervals. I would very

For time scale that are truly discrete Cox proposed the "exact partial likelihood". I call that the "exact" method and SAS calls it the "discrete" method. What we compute is precisely the same, however they use a clever algorithm which is faster. To make things even more confusing, Prentice introduced an "exact marginal likelihood" which is not

Hi, I am beginner with R and need firm guidance with my problem. I have seen some other threads discussing the subject of right censored data, but I am not sure whether or not this problem can be regarded as such. Data: I have a vector with laboratory test data (strength of wood specimens, example attached as txt-file). This data is the full sample. It is a common view that this kind of data

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130502/8bc15613/attachment.pl>

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130501/294047df/attachment.pl>

Dear all, all data I am talking about can be found at <http://www.metaanalyse.de/tmp/test.dat>. The R-code is located at <http://www.metaanalyse.de/tmp/test.R> I'd like to plot the frequencies for status against year. Typing table(tmp$year) it can be simple seen that some years are missing, e.g. 1975-1978, 1981 etc. Using xyplot to display the data I do not like the fact

Dear all, it seems to be that 'recode' can't handle any line breaks in its code. The following command causes no problem: datameta$smpid.r <- recode(datameta$smpid,"c(101,25,167,45,75)=25;c(104,51)=51") But if I type ... datameta$smpid.r <-recode(datameta$smpid, "c(101,25,167,45,75)=25; c(104,51)=51") ... the result is a syntax error.

Dear all, I am not able to install the lattice package. After finishing the installation, I got the following error massage: -------------------snip--------------------------- Delete downloaded files (y/N)? y Warning message: argument `lib' is missing: using /usr/local/lib/R/site-library in: install.packages("lattice") > library(lattice) Error in loadNamespace(i, c(lib.loc,

Dear all, I have two factors 'country' and 'status' which I would like to plot via barchart (lattice). 'status' consist of three different levels and should be the grouping variable, i.e. there should be drawn three different panels and within each panel a barchart of 'country'. barchart(daten$COUNTRY|daten$STATUS),

Pagarbiai Darius Evseicikas ================== UAB "PROGMERA" Maironio g. 3, LT-60149, Raseiniai tel./faks. +370 428 70329 mob. +370 699 97390 ICQ 297832308 Informacija šiame pranešime gali būti konfidenciali. Ji skirta tik asmeniui, kuriam yra adresuota. Jei Jūs nesate tas asmuo (arba atsakingas už šio pranešimo pristatymą tam asmeniui), Jūs negalite skaityti, kopijuoti ar bet kaip

Hi All, I am using RGui 2.3.0 library(stats) and using K-means Clustering to cluster data.matrix.(DDset) DDset (matrix 575 x 2) Exams (character) Slot (numeric) AAW326 35 BWA345 45 ............... ..... ZXE 453 67 DDset[,2] [1] 29 37 24 7 22 22 21 24 26 31 39 21 17 7 24 37 26 33 7 31 39 29 24 20 2 [26] 40 8 16 1 31 6 32 2

Dear Prof Harrell and everyone, My PC: Window XP service pack 3 and service pack 2 R version 2.8.1 and 2.9 alpha For the last 3 days, after updating R, my two computers have been facing problems when running existing and runable R commands that involves with Design package I attempt to use 'survest', but I failed all the times with R (both 2.8.1 and 2.9 alpha) being shut down

For an historical paper I'm working on, I have some Spanish plaintext, presently in the form of a Word .doc file, http://euclid.psych.yorku.ca/SCS/Gallery/images/Private/Langren/Verdadera-spanish-stripped.doc and also some ciphered text from the same original source. The ultimate goal is to use some frequency analysis of letters and word lengths in the plaintext to help decode the

Dear R mailing list, Greetings!!! I'm doing a graphical test of suitability for the aircondit data. I actually just followed the exercise in Davison and Hinkley's book of Graphical Test for model checking. The program for testing if it comes from an exponential distribution seems to be running okay but the program for the gamma distribution distribution is running away :-). I hope

Dear mailing list, I tried to run the example for the conditional bootstap written in the help file of censboot. I got the following result: STRATIFIED CONDITIONAL BOOTSTRAP FOR CENSORED DATA Call: censboot(data = aml, statistic = aml.fun, R = 499, F.surv = aml.s1, G.surv = aml.s2, strata = aml$group, sim = "cond") Bootstrap Statistics : original bias std. error t1*