similar to: Parametric survival models with left truncated, right censored data

Hi all, Sorry if this has been answered already, but I couldn't find it in the archives or general internet. Is it possible to implement the gompertz distribution as survreg.distribution to use with survreg of the survival library? I haven't found anything and recent attempts from my side weren't succefull so far. I know that other packages like 'eha' and

Dear list: I am attempting to use what I thought would be a pretty straightforward practical application of Cox regression. I figure users of the survival package must have come across this problem before, so I would like to ask you how you dealt with it. I have set up an illustrative example and included it at the end of this post. I took a sample of 100 data points from each of two populations

Hello, all. I have co-opted a number of functions that can be used to plot the hazard/survival functions and associated density distribution for a Gompertz mortality model, given known parameters. The Gompertz hazard model has been shown to fit relatively well to the human adult lifespan. For example, if I wanted to plot the hazard (i.e., mortality) functions: pop1 <- function (t) {

Hello. I am trying to write an R function to plot the survival function (and associated hazard and density) for a Siler competing hazards model. This model is similar to the Gompertz-Makeham, with the addition of a juvenile component that includes two parameters---one that describes the initial infant mortality rate, and a negative exponential that describes typical mortality decline over the

Hi all, I am having a problem trying to declare a slot representation to hold a file connection. With V2.8.0 I had been using: setClass("Element", representation( [other slots removed], datafileConn = setOldClass(c("file","connection")) ) ) This resulted in a warning but still ran okay. No however with V2.9.0 I am getting

My stepAIC function works for one set of data but not another....one set of data shows the steps of eliminating variables, versus another set of data doesn't throw away any variables. Can anyone please explain why? Thanks ______________________________________________________________________ The information transmitted, including any attachments, is intended only for the person or entity

See error message below: can someone please help with this? Thanks! Residuals: ALL 9 residuals are 0: no residual degrees of freedom! Residual standard error: NaN on 0 degrees of freedom Multiple R-squared: 1, Adjusted R-squared: NaN F-statistic: NaN on 8 and 0 DF, p-value: NA ______________________________________________________________________ The information transmitted,

Can someone please help with the below - thanks! Warning messages: 1: Not plotting observations with leverage one: 7 2: Not plotting observations with leverage one: 7 > print(qqPlot(fit),envelop=.95); Error in model.frame.default(formula = Y ~ X - 1, drop.unused.levels = TRUE) : variable lengths differ (found for 'X') In addition: Warning message: In matrix(yhat, n, reps) : data

Can someone please help with the error message below? thanks! Start: AIC=-Inf value ~ 1 + Core_CPI__ + GDP_change + Unemployment + housing + interest + S_P + d1 + d2 + d3 Error in if (any(ch)) { : missing value where TRUE/FALSE needed In addition: Warning message: attempting model selection on an essentially perfect fit is nonsense

After installing a new version of R, how do I make the packages that I have already installed and use on the old version available to the new version? Thank you, Roy Robertson

Can someone please help with the error message below? Warning messages: 1: Not plotting observations with leverage one: 7 2: Not plotting observations with leverage one: 7 > print(qqPlot(fit),envelop=.95); Error in model.frame.default(formula = Y ~ X - 1, drop.unused.levels = TRUE) : variable lengths differ (found for 'X') In addition: Warning message: In matrix(yhat, n, reps) :

Michael Robertson, founder of both MP3.com and Lindows, has launched a new company to supply inexpensive SIP phones ($129 for two) and related services. See today's press release at http://www.sipphone.com/tiki-index.php?page=SIPphone%20Inc My question for the list is who will be the first to report on the compatibility and usability of the SIPphone with Asterisk? The functionality

Hi all, We have been fighting with telco for a entire week. Today they came here with a LITE3000 to analyze what is going on. When I configure zaptel with no external clock, E1 gets aligned/synchronized with bit rate in 2048000 bps, both me and telco. span=4,0,0,ccs,hdb3,crc4 But when I configure span4 to get clock source from telco they become unsynchronized. TElco bit rate stays in

Dear Tomas, thank you for looking into this. Here's the output: # number of logical processors - what detectCores() should return out <- system("wmic cpu get numberoflogicalprocessors", intern=TRUE) [1] "NumberOfLogicalProcessors \r" "22 \r" "22 \r" [4] "20 \r"

Hi all, Prefix: I am a frustrated Java coder in R. I am coding a medium sized ecosystem modelling program in R. I have changed to using S4 objects and it has cost me an order of magnitude in execution speed over the functional model. I cannot afford this penalty and have found that it is the result of all the passing-by-value of objects. I see that you can now safely inherit from

fit lwr upr 1 218.4332 90.51019 346.3561 2 218.3906 90.46133 346.3198 3 218.3906 90.46133 346.3198 4 161.3982 44.85702 277.9394 5 192.4450 68.39903 316.4909 6 179.8056 56.49540 303.1158 7 219.5406 91.52707 347.5542 8 162.6761 46.65760 278.6945 9 193.8506 70.59838 317.1029 10 181.3816 58.11305 304.6502 11 221.2871 92.14366 350.4305 12 164.2947 47.91081 280.6785 13

Hello, I have a problem with some variables in 1.8.6.0. I set on extension the following lines: exten => h, n, Set (CDR (LLP) = $ {CHANNEL (rtpqos, audio, local_lostpackets)}) ; lost packets by local end ** exten => h, n, Set (CDR (PCR) = $ {CHANNEL (rtpqos, audio, remote_lostpackets)}) ; lost packets by remote end exten => h, n, Set (CDR (ljitt) = $ {CHANNEL (rtpqos,

Dear R users, I am using lmer combined with AIC model selection and averaging (in the MuMIn package) to try and assess how isotope values (which indicate diet) vary within a population of animals. I have multiple measures from individuals (variable 'Tattoo') and multiple individuals within social groups within 4 locations (A, B, C ,D) crucially I am interested if there are

Hi, I'm getting the following error: "The mapped network drive could not be created because the following error has occured - there are currently no logon servers available to service the logon request" Any ideas? Kind Regards, Mark Parsons. Ernst & Young is proud to bring you Entrepreneur Of The Year - the prestigious global business awards for entrepreneurs.

Hi All - Quick question. Is there anyway that when a user first invokes Samba to map to a network share that they are asked to change their (samba) password. I've got a whole bunch of (samba) users who all have been given the same password and I want them to have to change it to one of their own choosing upon first login. Kind Regards, Mark Parsons Ernst & Young is proud to bring