Untitled
unknown
plain_text
2 years ago
1.0 kB
4
Indexable
###OUTLIERS e ptos influentes#### boxplot(mrl.stepforward2$residuals) #potenciais outliers plot(mrl.stepforward2) #ver leverage plot outlierTest(mrl.stepforward2) #teste de bonferroni influencePlot(mrl.stepforward2) #hat value > 0.5 -> 99,113 cooksD <- cooks.distance(mrl.stepforward2) plot(cooksD, main = "Cook's Distance") abline(h = 4/nrows, col = "red") #9? influence <- influence.measures(mrl.stepforward2) obsinfluentes<- which(apply(influence$is.inf, 1, any)) dados_treino_sem <- dados_treino[-c(obsinfluentes),] y_treino2_sem <- y_treino2[-c(obsinfluentes)] mrl.stepforward2_sem <- lm(y_treino2_sem ~ Shock.Type + Sex + Diastolic.Pressure + Systolic.Pressure + Age + Body.Surface.Area + Sex:Body.Surface.Area + Sex:Systolic.Pressure + Sex:Diastolic.Pressure + Diastolic.Pressure:Body.Surface.Area,data=dados_treino_sem) summary(mrl.stepforward2_sem) extractAIC(mrl.stepforward2_sem) anova(mrl.stepforward2_sem, mrl.int2) formula(mrl.stepforward2_sem)
Editor is loading...