search for: muhaz

Dear R users, I am trying to use "muhaz" and "plot.muhaz" functions in "muhaz" package to estimate and plot hazard funciton. However function "muhaz" always gives error message "Error in Surv(times, delta) : object "times" not found". I could not even run their sample codes in t...

...QC problems, but have CRAN and BioC packages depending on them which would be broken by the archival (and hence need archiving alongside). Therefore we are looking for new maintainers taking over maintainership for one or more of the following packages: R2HTML SemiPar cghseg hexbin lgtdl monreg muhaz operators pamr Package maintainers whose packages depend on one of these may be natural candidates to become new maintainers. Hence this messages is addressed to all these maintainers via BCC and to R-devel. See <http://CRAN.R-project.org/package=R2HTML> <http://CRAN.R-project.o...

Hello, Is it possible to create confidence intervals for hazard rates? ?I'm creating two muhaz objects:? haz1 <- muhaz(NumDaysCustomer[cRV=="true"],status[cRV=="true"]) haz2 <- muhaz(NumDaysCustomer[cRV=="false"],status[cRV=="false"]) and plotting them. ? There are many, many more observations in the cohort cRV=="false" than ==&...

...w that if the hazard function is essentially constant, the mean lifetime is 1/hazard rate. ?But if I'm confident that the empirical hazard function is not constant, I'm not sure how to go about calculating an estimate of mean lifetime. For instance, suppose someone goes about generating a muhaz object by doing the following: >? haz <- muhaz(runif(1000,30,40)) and then gives me only the haz object. ?I can plot it, and see pretty quickly that it isn't constant. ?The person who generated the muhaz object would know that the average lifetime is approx. 35. ?Is there anyway that...

Hello. ?I am hoping someone will be willing to help me understand something about hazard plots created with muhaz(...). ?I have some background in statistics (minor in grad school), but I haven't been able to figure one thing about hazard plots. ?I am using hazard plots to track customer cancellations. ?I figure I can treat a cancellation as a "death", and if someone is still a customer today, th...

Dear All, I would like to ask if it is possible to estimate a hazard function using the muhaz command when all the data is right-censored. My data has information of the number of weeks people has been unemployed but all of them are unemployed at the date of the survey, that is, I cannot observed when the individuals leave the unemployment state. I appreciate your help, Best, Monica L.

...bqtl car cclust cfa chron cluster cmprsk coda conf.design cramer date diamonds dse e1071 ellipse event.chart exactRankTests fastICA fdim fields foreign fracdiff gafit gee geoR gld gregmisc gss ineq leaps lgtdl lmtest locfit logspline lokern lpridge maptree maxstat mclust mda meanscore mgcv mlbench muhaz multiv mvnmle mvtnorm netCDF nlme nlrq norm odesolve oz panel pcurve permax pinktoe pixmap polymars polynom princurve pspline quadprog quantreg rmeta rpart rpvm scatterplot3d sem sgeostat sm sma sn sna splancs sptests spweights strucchange survival tensor tree tripack tseries twostage vegan wavethr...

Dear colleagues, I have the following dataset. It is modelled on the data included in Box-Seteffenheiser and Jones "Event History Modelling" Using the following code, I try to find the baseline hazard function haz_1<-muhaz(bpa$time, bpa$censored, subset=(bpa$year=="2010" | bpa$ban=="1"), min.time=1, max.time=3) I think I'm doing everything right, but what I don't understand is how to derive a duration dependency coefficient rom the values contained in the muhaz object as per Box-Steffenhe...

...roduce a plot of log(relative risk) or relative risk (also known as hazard ratio) against the continuous variable. I have spent a long time looking for advice on how to do this but my search has proved fruitless - sorry if I've missed something obvious. It seems that there are options such as muhaz, survfit, coxph and cph that may enable some plots to be produced but none that specifically look at the relative risk one. In addition to the survival analysis, I have incorporated the mfp function (from package mfp). I currently use code such as, library(mfp) library(Design) coxfit1 <- cox...

Dear List, Is there an automated way to use the survival package to generate survival rate estimates and their standard errors? To be clear, *not *the survivorship estimates (which are cumulative), but the survival *rate * estimates... Thank you in advance for any help. Best, Brian [[alternative HTML version deleted]]

Hi, I want to plot non parametric estimates of the empirical hazard function for right censored data. I've tried many functions from different packages (muhaz, Design, survival, eha, event), but none of them gave me what I wanted. Am I missing something? Here's what I want. The data below is the same used by Kiefer (J. Economic Literature, 1988), which in turns use a subset of the data used by Kiefer (J. Econometrics, 1985) in his study of strikes d...

Hi all, Please help on this. We will be teaching epidemiology using opensource software. What are R built-in functions or functions in available packages that are capable of doing these: a) Logistic regression (glm?) b) Conditional logistic regression c) Logistic regression with random effects d) Beta-binomial regression e) Poisson regression f) Weibull regression (eha?) g) Exponential

...ger evaluate highr modelr progress rprojroot vctrs cli fansi hms modeltools promises rstudioapi viridisLite clipr farver htmltools mstate ps rvest webshot coda.base flexmix htmlwidgets muhaz purrr scales whisker colorspace flexsurv httr munsell qpdf selectr withr crayon flextable isoband mvtnorm quadprog spatial xfun crosstalk forcats jsonlite nlme R6 stringi...

...<- c(TEMP,"mlbench","mlica","mlmRev","mmlcr","mnormt","MNP","moc","modeltools","moments","monoProc","monreg","mprobit","MPV","MSBVAR","msm","muhaz","multcomp","multilevel","multinomRob","multtest","mvbutils","mvnmle","mvnormtest","mvoutlier","mvpart","mvtnorm","NADA","ncomplete","ncvar","negenes",&qu...

Hi everyone. There's been several threads on baseline hazard in Cox model but I think they were all on cumulative baseline hazard, for instance http://tolstoy.newcastle.edu.au/R/help/01a/0464.html http://tolstoy.newcastle.edu.au/R/help/01a/0436.html "basehaz" in package survival seems to do a cumulative hazard. extract from the basehaz function: sfit <- survfit(fit) H

...nAids (0.6-4), minet (2.4.0), minqa (1.1.8), miscTools (0.6-2), mitools (2.0.1), mixAK (1.0), mixdist (0.5-3), mixlow (1.0.0), mixOmics (2.7-1), mixtools (0.4.4), mlegp (3.1.1), mlmmm (0.3-1.2), mlogit (0.1-7), mokken (2.3), monomvn (1.8-3), mrdrc (1.0-5), msm (0.9.7), mstate (0.2.4), mtsc (0.0.2), muhaz (1.2.5), multcomp (1.1-7), multinomRob (1.8-4), multtest (2.4.0), mutatr (0.1.2), mvabund (0.1-10), mvbutils (2.5.1), mvtnorm (0.9-92), nanop (1.0-1), ncdf (1.6.1), ncvreg (2.1), NetCluster (0.2), NetData (0.2), networksis (1.4), nFactors (2.3.2), nFDR (0.0-1), nleqslv (1.6.1), nltm (1.4.1), NMF (0...