setwd("~/Desktop/storico_2008_ascii")
c08 <- read.csv("cons2008.csv")
summary(c08)
f08 <- read.csv("comp2008.csv")
p08 <- read.csv("peso2008.csv")
View(`f08`)
f08f <- f08[f08$CFDIC==1,]
regress <- lm(c08$CN ~ c08$Y2 + f08f$NPERC)
summary(regress)
regress <- lm(c08$CN ~ c08$Y2 + f08f$NPERC + f08f$NCOMP)
summary(regress)
regress <- lm(c08$CN ~ c08$Y2 + f08f$NPERC + f08f$NCOMP + as.factor(f08f$QUALP10))
summary(regress)
tapply(c08$Y2,list(f08f$NPERC,f08f$QUALP10),mean)
regress <- lm(c08$CN ~ c08$Y2 + f08f$NPERC + f08f$NCOMP + as.factor(f08f$QUALP10), weights= p08$PESOFL2)
summary(regress)
tab <- table(f08f$NPERC,f08f$QUALP10)
summary(tab)
tab1 <- tapply(p08$PESOFL2*c08$Y2,list(f08f$NPERC,f08f$QUALP10),sum)
tab2 <- tapply(p08$PESOFL2,list(f08f$NPERC,f08f$QUALP10),sum)
tabfin <- tab1/tab2
pairs(cbind(c08$CN, c08$Y2, f08f$NPERC))