similar to: Having trouble with plot.survfit and fun="cloglog"

Hello everybody, does anybody know how the function plot.survfit exactly works? I'd like to plot the log of the cummulative hazard against the log time by using plot.survfit(...fun="cloglog") which does not work correctly. The scales are wrong and there is an error message about infinit numbers. It must have something to do with the censored data, doesn't it? #Example:

windows xp R 2.8.1 I am trying to plot the -log(log(survival)) to visually test the proportional hazards assumption of a Cox regression. The plot, which should give two lines (one for each treatment) gives only one line and a warning message. I would appreciate help getting two lines, and an explanation of the warning message. My problem may the that I have very few events in one of my strata,

I'm having trouble using binomial(link="cloglog") with GLMM in lme4, Version: 0.5-2, Date: 2004/03/11. The example in the Help file works fine, even simplified as follows: fm0 <- GLMM(immun~1, data=guImmun, family=binomial, random=~1|comm) However, for another application, I need binomial(link="cloglog"), and this generates an error for me: >

I would like to do an R glm() with family = binomial(link="loglog") Right now, the cloglog link exists, which is nice when the data have a heavy tail to the left. I have the opposite case and the loglog link is what I need. Can someone suggest how to add the loglog link onto glm()? It would be lovely to have it there by default, and it certainly makes sense to have the two opposite

I am having problems with the package dse. I just installed R 2.2.1 and reinstalled all packages. I am running Windows XP Pro with all updates. Below there are two examples of error messages generated when trying to execute some simple programs. The code was taken directly from the package documentation. Any help on this will be greatly appreciated. Merry Christmas Fernando

hie when i use plot.survfit to plot more than one graph why I only see the last graph how do i see the other graphs.for example n=20 n1=n/2 n2=n/4 a11=4;a12=4 ;a21=4 ;a22=4 t1<-array(1,c(n1)) t2<-array(2,c(n1)) treatgrp=matrix(c(t1,t2))

Dear R-help, I am trying to create groupedData objects using the nlme library. I'm missing something basic, I know: Here is the first example in ch.1 of Pinheiro & Bates (2000): library(nlme) x2=Rail$travel;x1=Rail$Rail;eg1=data.frame(x1,x2);eg1gd=Rail print(eg1gd) x11();print(plot(eg1gd)) femodel=lm(x2~x1-1,data=eg1gd) print(femodel$coefficients) Result: x12 x15 x11

Hi, If I try to do this: if(is.R()) data("eg1.DSE.data.diff", package="dse1") model <- est.VARX.ls(eg1.DSE.data.diff) (Page 14 - DSE Package Manual) I obtain a Segment Violation. I use R-1.3.0 and the last dse package version Maximino Ameneiro Gomez -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I am trying to figure out how to apply a loglog link to a binomial model (dichotomous response variable with far more zeros than ones). I am aware that there are several relevant posts on this list, but I am afraid I need a little more help. The two suggested approaches seem to be: 1) modify the make.link function in GLM, or 2) use the loglog or cloglog functions in the VGAM package.

I don't off-hand remember, but the folks on the user list probably do (;-)) You accidentally sent to the development list --dave Serhat G?RSOY wrote: > *********************** > This mail and its attachments have been scanned by InterScan MSS. > ***********-*********** > > > hi folks, > > I am a new samba user and I have a problem like this; > > When I copy

Hello Experts I have the following codes and data for 2 interpolation plots. http://www.nabble.com/file/7206/3d_plot_data.txt 3d_plot_data.txt data<-read.table("3d_plot_data.txt", header=T) attach(data) par(mfrow=c(1,2)) library(akima) interpolation<-interp(rr,veg_r,predict) persp(interpolation,theta = -45, phi = 30, ticktype = "detailed", nticks=4, cex=0.8,

I'm trying to use GLMM in library(lme4), R 1.9.0pat, updated just now. I get an error message I can't decipher: library(lme4) set.seed(1) n <- 10 N <- 1000 DF <- data.frame(yield=rbinom(n, N, .99)/N, nest=1:n) fit <- GLMM(yield~1, random=~1|nest, family=binomial, data=DF, weights=rep(N, n)) Error in eval(expr, envir, enclos) : ..1 used in an incorrect

hi, Just a very simple question: is there an R equivalent to the matlab command find(X) which returns the indices of vector X that store non-zero elements? e.g. > find( [1 0 0 1 0]) ans = 1 4 so, in R, how do I do: ans <- rfind( c(1,0,0,1,0)) so that ans is the vector c(1,4) thanks, stephen -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

Dear R users, I try to write an rcode for a paranetric weibull model with unobserved heterogeneity. The data I have are giuven below. ( The data comes from Keiding and Klein artricle published in the Statistics in Medicine 1996) . While the simple weibull model runs perfectly ( I get -2*loglikelihood 1093,905), the weibull frailty model does not run well. The code I write is :

Hi there! I set up a test enviroment and did a log level 256 of a domain user log on with a Samba 3.5.5 and a windows 7 machine to deliver more data for that _netr_ServerAuthenticate3: netlogon_creds_server_check failed. Rejecting auth request from client TESTPC1 machine account TESTPC1$ problem. Please find the logs enclosed (if a ~300kb rar log file is allowed here) and maybe somebody

I have a need to build a time series and there are a couple of aspects about the time series object that are confusing me. First it seems that ts.union is not doing what I would expect. For example: x0 <- rep(0,10) x1 <- rep(1,10) xt0 <- ts(x0, frequency=10) xt1 <- ts(x1, frequency=10) st2 <- ts.union(xt0, xt1) > xt2 Time Series: Start = c(1, 1) End = c(1, 10) Frequency = 10

Hello, The following code: > binomial(identity) Generates an error message: Error in binomial(identity) : link "identity" not available for binomial family; available links are ?logit?, ?probit?, ?cloglog?, ?cauchit?, ?log? While : > binomial("identity") Yields an identity-binomial object that works as expected with stats::glm The error in the first example mislead

Hello, I have different results from these two softwares for a simple binomial GLM problem. >From Genmod in SAS: LogLikelihood=-4.75, coeff(intercept)=-3.59, coeff(x)=0.95 >From glm in R: LogLikelihood=-0.94, coeff(intercept)=-3.99, coeff(x)=1.36 Is there anyone tell me what I did wrong? Here are the code and results, 1) SAS Genmod: % r: # of failure % k: size of a risk set data

Hello, I'm not a Stata user so I'm trying to reproduce Stata results that are given to me in R. I would like to use a GLM with a complementary log-log function. The stata code I have is: glm c IndA fia, family(binomial s) link(cloglog) offset(offset) The R code is: glmt <- glm(data=dataset, c ~ IndA + fia, offset = offset, family = binomial(link = cloglog)) Which yields

Hi Hadley and everyone, here's a patch for ggplot2 that fixes the behavior of opts(legend.position={left,top,bottom}). If you try the following code in an unmodified ggplot2 options(warn = -1) suppressPackageStartupMessages(library("ggplot2")) data <- data.frame( x = c(1, 2, 3, 4, 5, 6), y = c(2, 3, 4, 3, 4, 5), colour = c(TRUE, TRUE, TRUE, FALSE, FALSE, FALSE))