Dear All!
I'm trying to determine if a plant genotype is resistant or susceptible to an insect. For that purpose, I made a test with five insects in a Petri dish. Inside Petri dish, was one leaf from one plant of genotype C1. We used five plants from each Genotype and three leaves from each plant. Here is the subset of my data https://www.dropbox.com/scl/fi/6k8gd17t9qa7xaj5t9z0h/SUZ.Sobrevivencia.l.dat?rlkey=g7mg5i345eb8zng8xqjz7wwq3&dl=0
I trying to modeling it in Survival package.
library(survival)
library(survminer)
SUZ.Sobrevivencia.l.dat <- read.table("SUZ.Sobrevivencia.l.dat", header=TRUE, dec=",", sep="\t")
head(SUZ.Sobrevivencia.l.dat)
SUZ.survival <- survfit(Surv(as.numeric(Day), Alive==5) ~ Genotype, data = SUZ.Sobrevivencia.l.dat)
surv_summary(SUZ.survival)
ggsurv <- ggsurvplot(SUZ.survival, data = SUZ.Sobrevivencia.l.dat)
ggsurv$plot + theme_bw()+facet_wrap(~strata)
My doubts are if this approach is suitable to:
1. determine if Genotype is resistant or susceptible;
2. how to to a threshold to divide genotypes in four groups: Resistant, moderately resistant, moderately susceptible, and susceptible.
If this approach is not suitable for that purpose, can you advise in a good way to do it?
Thank you co much!