Progrmmazione su R caso di studio, Esercizi di Statistica Applicata

Scarica Progrmmazione su R caso di studio e più Esercizi in PDF di Statistica Applicata solo su Docsity! Codici utilizzati su R #Inserimento librerie library(gam) library(readxl) library(olsrr) library(fBasics) #Inserimento dati Data <- read.csv(file.choose()) #Trattamento valori mancanti #aggiunta media ai valori mancanti (N/A) EPData <- na.gam.replace(Data) #Analisi descrittiva basicStats(data.frame(EPData$Life_expectancy)) basicStats(data.frame(EPData$Adult_Mortality)) basicStats(data.frame(EPData$HIV)) basicStats(data.frame(EPData$Income)) basicStats(data.frame(EPData$T_expenditure)) basicStats(data.frame(EPData$Diphtheria)) #Status (Variabile fattoriale) Frequenza tab1(EPData$Status, short.group = "decreasing", cum.percent = TRUE) #Cambio in fattoriale per analisi correlazione e regressione EPData$Status=as.factor(EPData$Status) EPData$Status=as.integer(EPData$Status) #Trovare il miglior sottoinsieme model1 <- lm(Life_expectancy~Adult_Mortality+Infant_deaths+Alcohol+Percentage_expenditure+HepatitisB+ Measles+BMI+Under.5Death+Polio+T_expenditure+Diphtheria+HIV+GDP+Population+thinness1_19 +thinness5_9+Income+Schooling+Status,data =EPData) print(model1) FWDfit.p<-ols_step_forward_p(model1,penter=.05) #controllo delle variabili da accettare/selezionare FWDfit.p #Correlazione cor.test(EPData$Life_expectancy,EPData$Adult_Mortality) cor.test(EPData$Life_expectancy,EPData$HIV) cor.test(EPData$Life_expectancy,EPData$Income) cor.test(EPData$Life_expectancy,EPData$T_expenditure) cor.test(EPData$Life_expectancy,EPData$Diphtheria) cor.test(EPData$Life_expectancy,EPData$Status) #Regressione lineare multipla model2 <- lm(Life_expectancy~Adult_Mortality+HIV+Income+T_expenditure+Diphtheria+Status,data =EPData) print(model2) summary(model2) #Grafici Dataset1 <- as_tibble(EPData) mydataset <- Dataset1 %>% select(Life_expectancy,Adult_Mortality,HIV,Income,T_expenditure,Diphtheria,Status) #Boxplot boxplot(EPData$Life_expectancy, col= "lightblue", main= "Boxplot - Life Expectancy")