library(foreign) base=read.spss("2poblaciones-Mismotrabajo-DiferenteNutricion.sav",F,T) color=c("#FF000044","#0000FF44") colesterol=base$Colesterol grupo=base$Grupo+1 plot(colesterol,rep(1,length(colesterol)),col=color[grupo+1]) # library(foreign) df=read.spss("supervivencia.sav",F,T) names(df) df$t=df$tiempo df$status=df$estado df$grupo=df$tumorcat superv=survfit(Surv(t, status) ~ grupo, data = df,type = "kaplan-meier") qq=survdiff(Surv(t, status) ~ grupo, data = df) names(qq) pchisq(qq$chisq,df = 2) qq$n qq$obs qq$exp qq$var qq$chisq library(tidyverse) library(openxlsx) library(haven) df=read.xlsx("coolKids.xlsx") %>% as_factor() %>% as_tibble() %>% mutate(Sexo=factor(Sexo), Curso=integer(Curso), Edad=integer(Edad), Raza=factor(Raza), Entorno=factor(Entorno) )