similar to: ayuda con análisis de supervivencia

Queridos R-users: Tengo una duda que hace mucho tiempo que estoy intentando resolver, os explico a modo de ejemplo: Tengo estos efectos: Año(5 niveles),Localidad (10 niveles) y genotipo (3 niveles), año y localidad son aleatorios y genotipo es fijo (los he escogido yo). Me gustaría hacer obtener una tabla parecida a la Tabla Anova donde aparezca cada factor y sus interacciones y

Dear R users, Given this data: x <- seq(1,100,1) dx <- as.POSIXct(x*900, origin="2007-06-01 00:00:00") dfx <- data.frame(dx) Now to play around for example: subset(dfx, dx > as.POSIXct("2007-06-01 16:00:00")) Ok. Now for some reason I want to extract the datapoints between hours 10:00:00 and 14:00:00, so I thought well: subset(dfx, dx >

I am having trouble formatting some survival data to use in a time dependent cox model. My time dep. variable is habitat and I have it recorded for every day (with some NAs). I think it is working properly except for calculating the death.time. This column should be 1s or 0s and as I have it only produces 0s. Any help will be greatly appreciated.

Hi all, I apologize for this probably stupid question, but I really can't figure it out. I have a dataframe like this: group <- c(rep('A', 8), rep('B', 15), rep('C', 6)) time <- c(rep(seq(1:4), 2), rep(seq(1:5), 3), rep(seq(1:3), 2)) value <- runif (29, 1, 10) dfx <- data.frame (group, time, value) I want to calculate mean and standard deviation for all

Hi, I am trying in R to indicate in which quintile a value of a variable is for every month of my data frame in this case based on volatility. For each month I want to know for each stock if it is in the most volatile quintile of if it is in one of the others. So far I have come up with the following function (see below). Unfortunately, the function only works in some cases and often gives the

Dear R users, My question is more methodology related rather than specific to R usage. Using time on study as time in a cox model, eg: library(Design) stanf.cph1=cph(Surv(time, status) ~ t5+id+age, data=stanford2, surv=T) #In this case the 1000-day survival probability would be: stanf.surv1=survest(stanf.cph1, times=1000) #Age in this case is a covariate. #I now want to compare the above

Dear All, How can I sort a data frame with respect to more than one variable? I know that for one variable X one may use: df[order(df$X), ] where df is the data frame containing X. Many thanks, Bernard --------------------------------- [[alternative HTML version deleted]]

Dear all, I was wondering if anyone can help? I am an R user but recently I have resorted to SAS to calculate the probability of the event (and the associated confidence interval) for the Cox model with combinations of risk factors. For example, suppose I have a Cox model with two binary variables, one for gender and one for treatment, I wish to calculate the probability of survival for the

I tried: n <- 9 x <- sample(n) names(x) <- LETTERS[1:n] dotplot( sort(x) ) The lines of the dotplot are not ordered according to the values of x (but according to the names of x). So I did: dfx <- data.frame( x=x, label=factor( names(x), names(x)[order(x)] ) ) dotplot( label ~ x, data=dfx ) So I got what I wanted. My question: Is there an easier solution

This is just a curiosity question. Why do the two different syntaxes for df1 and df2 give such different results in the names(dfx)? Thanks df1 <- data.frame(nas = c("A", "B" , "B" ,"C" ,"B", "A", "D"),nums = c(3, 2, 1, 1, 2, 3, 7)) df2 <- data.frame(nas <- c("A", "B" ,

Buenas tardes, Estamos intentando hacer un estudio sobre diversidad en gatos callejeros ('stray cats') de la base de datos 'nancycats' (librería 'adegenet') de R. Sin embargo, no entendemos muy bien los datos. ¿Alguien puede explicarnos a qué corresponde la variable respuesta? Entendemos que las filas son los 237 gatos (observaciones), y las columnas se corresponden con

Estimado Gemma Ruiz Olalla De curioso miré una búsqueda en internet, y se la comparto, de esta observe la página 8 http://adegenet.r-forge.r-project.org/files/tutorial-basics.pdf Estoy de acuerdo con Jorge y Carlos, pero por las dudas, es solo una codificación, yo sin mirar más no interpretaría que el heterocigota tiene un "peso" de 0.5, porque si hay sobredominancia, dominancia, etc

Hi, I have the binned data (observed and generated from model)  that I would like to test using the anderson-darling goodness of fit test.  But I'm not sure which package in R to use. I tried ad.test(...) but it does not recognise the test by Vito Ricci in FITTING DISTRIBUTIONS WITH R   > ad.test(hist_hume_beec[,1],hist_hume_beec[,2]) Error: could not find function "ad.test"

Hi I'm having a problem with model.frame, encapsulated in this example: y1 <- matrix(c(3,1,0,1,0,1,1,0,0,0,1,0,0,0,1,1,0,1,1,1), nrow = 5, byrow = TRUE) y1 <- as.data.frame(y1) rownames(y1) <- paste("site", 1:5, sep = "") colnames(y1) <- paste("spp", 1:4, sep = "") y1 model.frame(~ y1) Error in model.frame(formula, rownames,

Hola, ¿qué tal? Sabréis que en unos días comienza un MOOC de Hastie y Tibshirani, https://class.stanford.edu/courses/HumanitiesScience/StatLearning/Winter2014/about basado en el libro "An Introduction to Statistical Learning, with Applications in R". Se me ocurre lo siguiente: 1) Montar un grupo de estudio en Madrid de manera que podamos reunirnos periodicamente (¿semanalmente?

I haven't done any R memory debugging lately, but https://www.mail-archive.com/rcpp-devel at lists.r-forge.r-project.org/msg10289.html shows how I used to have gdb break where valgrind finds a problem so you could examine the details. Also, running your code after running gctorture(TRUE) can help track down memory problems. -Bill On Wed, Feb 7, 2024 at 12:03?PM Therneau, Terry M., Ph.D.

?I've hit a roadblock debugging a new update to the survival package.?? I do debugging in a developement envinment, i.e. I don't create and load a package but rather? source all the .R files and dyn.load an .so file, which makes things a bit easier. ? Running with R -d "valgrind --tool=memcheck --leak-check=full" one of my test files crashes in simple R code a dozen lines

https://bugs.freedesktop.org/show_bug.cgi?id=46810 Bug #: 46810 Summary: Mesa Failing To Build Classification: Unclassified Product: Mesa Version: git Platform: x86-64 (AMD64) OS/Version: Linux (All) Status: NEW Severity: normal Priority: medium Component: Drivers/DRI/nouveau

Hola, Dado que hay varios matemáticos en el foro, ?alguien puede recomendar algún buen libro en español o ingles que ayude a entender el álgebra lineal de primero de matemáticas ?. Gracias y un saludo JM Arbones

Dear useRs, There is a new version 0.1-2 of the package surv2sample available on CRAN. Users of the previous versions should update because a bug in the function cif2.ks has been fixed. General information about the package: surv2sample provides various two-sample tests for right-censored survival data. Three main areas and corresponding methods are: * comparison of two survival