EAE_ayuda_trabajo_2018, R

Language: Layout:

media<-86
desv<-6.2
id<-c(1:60)

actitud<-round(rnorm(60, media, desv),0)

sexo<-rep(sample(c("H","M"),size=30, prob = c(0.35,0.65), replace = T),2)
nota<-c(round(rnorm(30, 6.7, 1.2),2), round(rnorm(30, 5.8, 1.1),2))

viaAcceso<-rep(sample(c("GSFP","M25","PAU","Otros"),size=30, prob = c(5,7,10,8)/30, replace = T),2)

grupo<-c(rep("E",30),rep("C",30))

df<-data.frame(id, actitud, sexo, nota, viaAcceso, grupo)
df$sexo<-factor(sexo, levels =c("H","M"))
df$viaAcceso<-factor(viaAcceso,levels=c("GSFP","M25","PAU","Otros"))

cat("\n","listar los 10 primeros registros", "\n")
head(df, 10)

cat("\n","listar los 10 ultimos registros", "\n")
tail(df, 10)

# estudio descriptivo
cat("\n","*** estudio descriptivo", "\n")
summary(df[, 2:6])
cat("varianza Actitud ", round(var(df$actitud),3), "\n")
cat("varianza Nota ", round(var(df$nota),3), "\n")

# variable Actitud
cat("\n","*** variable Actitud", "\n")
stem(df$actitud/10)
hist(df$actitud, ylim=c(0,0.1), probability = T,main="Distribución de la actitud hacia la estadística", xlab="Actitud hacia la Estadística", ylab="Densidad")
curve(dnorm(x, mean=mean(df$actitud), sd=sd(df$actitud)), add=TRUE, col="red")

# boxplot(df$actitud~df$grupo, main="Distribución de la actitud hacia la estadística")

# variable sexo
cat("\n","*** variable sexo", "\n")
table(df$sexo)
slices<-table(df$sexo)
lbls<-names(slices)
pct<-round(slices/sum(slices)*100,1)
lbls <- paste(lbls, pct) # add percents to labels
lbls <- paste(lbls,"%",sep="") # ad % to labels
pie(slices, col=c("blue", "red"), main="Distribución del Sexo", labels = lbls)

# variable nota
cat("\n","*** variable nota", "\n")
stem(df$nota/10)
boxplot(df$nota, ylim=c(0,10), main="Distribución de las Notas Finales")
text(1,mean(df$nota), "*", col="red")

# boxplot(df$nota~df$grupo, ylim=c(0,10), main="Distribución de las Notas Finales")

# variable ViaAcceso 
cat("\n","*** variable ViaAcceso ", "\n")
table(df$viaAcceso)
barplot(table(df$viaAcceso), main="Distribución por forma de acceso")

# variable grupo
cat("\n","*** variable grupo", "\n")
table(df$grupo)
slices<-table(df$grupo)
lbls<-names(slices)
pct<-round(slices/sum(slices)*100,1)
lbls <- paste(lbls, pct) # add percents to labels
lbls <- paste(lbls,"%",sep="") # ad % to labels
pie(slices, col=c("orange", "green"), main="Distribución del grupo", labels = lbls)

[ + ] Show input
Absolute running time: 1.56 sec, cpu time: 0.95 sec, memory peak: 34 Mb, absolute service time: 1,59 sec