search for: km

Dear R users: I want to convert some character vectors into numeric vectors. > head(price) [1] "15450 EUR" "7900 EUR" "13800 EUR" "3990 EUR" "4500 EUR" [6] "4250 EUR" >head(mileage) [1] "21000?km" "119000?km" "36600?km" "92000?km" "140200?km" [6] "90000?km" in the first example I can use: price <- sub(" EUR", "", price) to get "15450" "7900" "13800" "3990" "4...

Dear R users: I encountered difficulties in michaelis-menten equation. I found that when I use right model definiens, I got wrong Km vlaue, and I got right Km value when i use wrong model definiens. The value of Vd and Vmax are correct in these two models. #-----right model definiens-------- PKindex<-data.frame(time=c(0,1,2,4,6,8,10,12,16,20,24), conc=c(8.57,8.30,8.01,7.44,6.88,6.32,5.76,5.20,4.08,2.98,1.89)) mm.mo...

I am new to the world of R/programming so this may be a really easy question. I thank you for your patience and help in advance I would like the characters km^2 to be displayed on the plot subtitle as km squared - two as a superscript. I would also like to have the numbers from the data set for longitude and latitude to be rounded to four decimal places. Thank you. plot ( decade[['date']], decade[['value']], type = 'l',...

...ta and code below. I am plotting the number of trials and the number of successes in the spatial binomial experiments and would like to draw the spatial cells were the trials and successes were counted, i.e. a partial grid in the plot only for those cells where there is a number. The cells are 2x2 km cells. The count of Trials and Success should ideally appear in the middle of the square cell. I know there is the 'grid' package but it seems the plots made using 'graphics' are not compatible with the plots made using 'grid' (as warned in the grid help pages). Thanks in a...

Dear R-help list, I have a problem with the tick marks of a Kaplan-Meier survival plot. Here is a sample: follow.up<-c(10,20,30,40,50,60,70,80,90,100) #months dead<-c(1,1,1,0,1,1,0,0,0,0) KM <-survfit(Surv(follow.up, dead)) plot(KM) The result is a nice plot. However, our research group thinks it may be a better idea to place the ticks to the years on the time scale, i.e. 0, 12, 24, 36 etc. months. Is this possible with R? I tried to look it up and the most relevant manual page wa...

I'm testing R version 2.7.0 on windows and there seems to be a compatibility issue with objects that were created by "approxfun" in older versions. As long as the objects were created in version 2.7.0 things work ok, but calling the interpolated functions from R version 2.0.1 causes this error: Error in .C("R_approx", as.double(x), as.double(y), as.integer(n), xout =

Hello again, I'm able to run KoolMoves 6.0 (KM) under Wine but it's missing some functionality. KM is a Flash authoring tool and the missing functionality is mostly the ability to use the internal preview functions. Internal preview requires a browser and a flash plug-in, both of which I have. I also tried KM under CrossOver Office and i...

...stor.Mine [i] = summary(data_Pb.boot)$coefficients[5,1] } hist(stor.r2, xlab="R-squared",main="Distribution of R-squared - Lead (log)") hist(stor.inter, xlab="Intercept",main="Distribution of Intercept - Lead (log)") hist(stor.Ind5, xlab="Industrial 5 km",main="Distribution of Industrial 5 km - Lead (log)") hist(stor.LNPRI25, xlab="NPRI 25 km (log)",main="Distribution of NPRI 25 km - Lead (log)") hist(stor.NPRI10, xlab="NPRI 10 km",main="Distribution of NPRI 10 km - Lead (log)") hist(stor.Mine,...

Dear All, I am trying to put 10^-8 st km^-2day^-1 on x-axis of my plot. I tried using : ylab = expression(paste("st / ", plain(km)^2, " / day")) to see if I can at least get the unit before thinking about the power of 10 (10^-8). However, ylab = expression(paste("st / ", plain(km)^2, " / day")) di...

...the scaled schoenfeld residuals and it also provides some built-in ones and the capability to provide your own. It also appears to me that different forms look at different departures from proportionality. So, my question is what are the advantages and disadvantages of the default transform="km" compared to say, identity or log? Thank you. Kevin -- Kevin E. Thorpe Biostatistician/Trialist, Knowledge Translation Program Assistant Professor, Department of Public Health Sciences Faculty of Medicine, University of Toronto email: kevin.thorpe at utoronto.ca Tel: 416.946.8081 Fax: 41...

...última revisión de la cuchilla > 3. TIEMPO: el tiempo transcurrido entre las dos fechas anteriores > 4. ESTADO: estado de la cuchilla cuando se revisó la última vez: > 0 = Buen estado > 1 = Para tirar > > Y utiliza el método kaplan-meier: > > library(survival) > KM <- survfit(Surv(DATO$TIEMPO, DATOS$ESTADO) ~ 1) > summary(KM) > plot(KM) > > Un ejemplo: > > > library(survival) > > AML.KM <- survfit(Surv(aml$time, aml$status) ~ 1) > > summary(AML.KM) > Call: survfit(formula = Surv(aml$time, aml$status) ~ 1) > >...

Hello, This is Elaine. I am drawing a plot with x-axis label with "km square" as the unit. Now I want to print km square in the form of km2 and output "2" as the uppercase. Please kindly help suggest command to show the uppercase. Thank you. Elaine [[alternative HTML version deleted]]

What is the best way to report the standard error when publishing Kaplan-Meier plots? In my field (Vascular Surgery), practitioners loosely refer to the "10% error" cutoff as the point at which to stop drawing the KM curve. I am interpreting this as the *standard error of the cumulative hazard*, although I'm having a difficult time finding some guidelines about this (perhaps I am not searching the correct terms or references). My KM figures contain typically two curves that I am comparing using the logran...

I am trying to get the following line in a plot margin using mtext: 100% Area = 120.000 km^2 Where I intend that "100% Area =" is text, 120.000 is a number that varies according to different data, and km^2 should be a neat km-with-superscript-2. The expression function fails me, since it apparently cannot coerce an expression when a variable is involved. In the R help file t...

Hi, I am wondering if anyone can explain to me if cumulative incidence (CI) is just "1 minus kaplan-Meier survival"? Under what circumstance, you should use cumulative incidence vs KM survival? If the relationship is just CI = 1-survival, then what difference it makes to use one vs. the other? And in R how I can draw a cumulative incidence plot. I know I can make a Kaplan-Meier survival plot using plot(survfit()), for example: fit<-survfit(Surv(time,status)~group,data=dat...

Buenas, Como pudeo calcular el tiempo de vida? Os cuento, tengo una serie de cuchillas y quiero ver el consumo de las mismas y he pensado en hacer un estudio por tiempo de vida. No se como hacerlo con R Gracias Jesús [[alternative HTML version deleted]]

Hi All, I'm trying to figure out the cumulative incidence curve in R in some limited time. I found in package "cmprsk", the command "plot.cuminc" can get this curve. But I noticed that there is no mark for the censored time there, comparing with the KM curve by "plot.survfit". Here are my codes (attached is the data): ---------------- dat<-read.table("F://wendy/BMT data analysis/final data.txt",header=TRUE,sep="\t") library(cmprsk) library(survival) attach(dat) par(mfrow=c(2,1)) curve<-...

...ka. > passwd-file). Application (or rather device) specific passwords are > implementing by using having an additional "username" with a specific > password for a particular application or device. E.g. some entries for > myself: > > bbmutt:*:10001:10001::/krot/mail/km::userdb_mail=maildir:~/Maildir userdb_quota_rule=*:bytes=10240M > kmozilla:*:10001:10001::/krot/mail/km::userdb_mail=maildir:~/Maildir userdb_quota_rule=*:bytes=10240M > sailpad:*:10001:10001::/krot/mail/km::userdb_mail=maildir:~/Maildir userdb_quota_rule=*:bytes=10240M > wo...

Hi all, I'm plotting to get the intersection value of three curves. Defining the x-axis as dsm, the following code works; dsm = c(800,600,NA,525,NA,450,400,NA,NA,NA,0) s3 = seq(0.05,1.05,0.1) plot(dsm,s3,col="blue",las=1,ylab="fraction",xlab="distance (km)") fc <- function(x,a,b){a*exp(-b*x)} fm <- nls(s3~fc(dsm,a,b),start=c(a=1,b=0)) co <- coef(fm) curve(fc(x,a=co[1],b=co[2]),add=TRUE,col="black",lwd=1) r <- range(dsm,na.rm=TRUE) n <- 1/2.71 val <- optimize(f=function(x) abs(fc(x,a=co[1],b=co[2])-n),c(r[1],r[2]))...

...quot;) 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: motorette <- read.table("http://stat.ethz.ch/Teaching/Datasets/NDK/motorette.dat",header=T) km <- survfit(Surv(time,status)~1,data=motorette,type="kaplan-meier") plot(km,fun="cloglog") # Error message and wrong scales plot(km,fun="cumhaz",log="xy") # Does the same. #If I do the plot "by hand", it does work: log.cum.H <- log(-log(km$s...