similar to: CCDF plotting - axis labels+overlay

Dear all, I would like to calculate the complementary cumulative distribution function. As it is known, the ccdf is the 1-ecdf(X)==1-F(x). (From ?ecdf help is shown that ecdf returns a function which is the F(x) I would like to use that function inside my function and after call curve (which accepts function as input). My code looks like that onVector<-seq(1:100)

Dear all, I'm trying to fit a curve to some 1 year failure-time data, so that I can extrapolate and predict failure rates up to 3 years. The data is in the general form: Treatment Time Status Treatment A 28 0 Treatment B 28 0 Treatment B 28 0 Treatment A 28

Hello, I have a survivor curve that shows account cancellations during the past 3.5 months.  Fortunately for our business, but unfortunately for my analysis, the survivor curve doesn't yet pass through 50%.  Is there a safe way to extend the survivor curve and estimate at what time we'll hit 50%? We started a new program 3.5 months ago, and I believe that this set of accounts behaves

I have a question about constructing the likelihood function where there is censoring at level 1 in a two-level random effects sum. In a conventional solution, the likelihood function is constructed using the density for failures and the survivor function for (in this case, right) censored results. Within (for example) an R environment, this is easy to do and gives the same solution as survreg

My apologies for this question but I?m stuck and I?m sure that there must be an easy answer out there (and hope that someone will have mercy and point me in the right direction). I have a data file that looks like: 1 77 3 1 8 1 1 7 2 1 1 5 1 42 7 1 0 2 1 23 1 2 83 9 2 8 2 2 6 5 2 23 3 3 11 3 3 8 1 . .... etc. . N 3 2 (FWIW, these are document, word reference, and word frequency counts.) I

I am experimenting with parallel processing using foreach and seem to be missing something fundamental. Cool stuff. I've gone through the list and seen a couple of closely related issues, but nothing I've tried seems to work. I know that the results from foreach are combined, but what if there is more than one variable within the loop? Below is a snippet (non-functioning) of code that I

Version 2.4.0 of Rattle has been released to CRAN. The rattle package (http://rattle.togaware.com) is a multi platform (GNU/Linux, Mac/OSX, MS/Windows) GTK based GUI for data mining (for exploring data and building descriptive and predictive models). It has undergone a lot of development over the past year. A companion book introducing data mining using Rattle is under development, with a draft

Version 2.4.0 of Rattle has been released to CRAN. The rattle package (http://rattle.togaware.com) is a multi platform (GNU/Linux, Mac/OSX, MS/Windows) GTK based GUI for data mining (for exploring data and building descriptive and predictive models). It has undergone a lot of development over the past year. A companion book introducing data mining using Rattle is under development, with a draft

The Building Painting (http://www.micaroo.com/home-and-garden/wall-art/building-paintings.html) Farmhouse in Provence also known as Entrance Gate to a Farm with Haystacks was made in 1888 by Vincent van Gogh in Arles in Provence at the height of his career. Oil Paintings for sale (http://www.micaroo.com/home-and-garden/wall-art.html) Partially due to having been inspired by painter Adolphe

I've been trying various ways of getting wine running on my 64 bit ubuntu machine. So far I haven't been very successful. I am trying the source for 0.9.20 and I'm using the latest k8 kernel in ubuntu. I have an nvidia card and am using the latest glx drivers available in ubuntu (1.0.8762). I am using the guide here: http://wiki.winehq.org/WineOn64bit with help from here:

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;

Please find release 3.5.0 of the CRAN package "randomSurvivalForest" now posted on CRAN. Thank you. ubk2101 at columbia.edu Udaya B. Kogalur, Ph.D. Kogalur Shear Corporation 5425 Nestleway Drive, Suite L1 Clemmons, NC 27012 --------------------------------------------------------------------------------- CHANGES TO RELEASE 3.5.0 RELEASE 3.5.0 represents a significant upgrade in the

Please find release 3.5.0 of the CRAN package "randomSurvivalForest" now posted on CRAN. Thank you. ubk2101 at columbia.edu Udaya B. Kogalur, Ph.D. Kogalur Shear Corporation 5425 Nestleway Drive, Suite L1 Clemmons, NC 27012 --------------------------------------------------------------------------------- CHANGES TO RELEASE 3.5.0 RELEASE 3.5.0 represents a significant upgrade in the

Dear R users, I am using the vars package to calculate the impulse response functions and the forecast error variance decomposition of a VAR model. Unfortunately I do not know whether these functions assume unit or one standard deviation shocks. I tried to look into the code of these functions, but in vain: neither irf, nor vars::irf, nor vars:::irf output the code of the functions. Does someone

Hi all. I may be wrong, (and often am), but in trying to determine how to calculate the erf function, the documentation for 'pnorm' states: ## if you want the so-called 'error function' erf <- function(x) 2 * pnorm(x * sqrt(2)) - 1 ## and the so-called 'complementary error function' erfc <- function(x) 2 * pnorm(x * sqrt(2), lower=FALSE) Should, instead, it read:

Dear R users, I need to to the the following. Let a= 1 2 3 4 5 6 and b= -1 -2 -3 be (2x3) matrices. -4 -5 -6 I need to combine the two matrices into a new (2x6) matrix like this: ab = ( 1 -1 2 -2 3 -3 ) 4 -4 5 -5 6 -6 How can this be done in R? ----------------------------------------------------------------- ?????? ???

Hi, I'm using Ubuntu 8.04 with samba 3.028a (ubuntu package) +LDAP and I'm getting an odd problem. I have a user and this user had 4 groups (1 principal and 3 complementary) if I run id <user> in the console of samba server the command returns 4 groups (so LDAP is working fine). But when I'm trying to access a samba share (with log level 5), I noticed that samba returns that

Dear R users, Please can you help me with a relatively straightforward problem that I am struggling with? I am simply trying to plot a baseline survivor and hazard function for a simple data set of lung cancer survival where `futime' is follow up time in months and status is 1=dead and 0=alive. Using the survival package: lung.wbs <- survreg( Surv(futime, status)~ 1, data=lung,

Dear all, What is the cumulative distribution (with parameterization) used within survreg with respect to the log-logistic distribution? That is, how are the parameters linked to the survivor function? Best regards, Mario [[alternative HTML version deleted]]

Dear R users, I need perform structural analysis on a no intercept VAR model. Unfortunately the functions irf.VAR and dfev that come with the MSBVAR package only work with objects output by the reduced.form.var function, which seems to only evaluate VAR models with intercept. Is there a way to suppress the estimation of intercept term in reduced.form.var? Do I need to modify the code, and if I