Che materia stai cercando?

Anteprima

ESTRATTO DOCUMENTO

#################################

# regressione lineare semplice #

#################################

data(longley)

dati1<-longley

fix(dati1)

summary(dati1)

names(dati1)

attach(dati1)

plot(Employed,GNP)

mod1<-lm(GNP~Employed)

coef(mod1)

abline(mod1)

summary(mod1)

n<-nrow(dati1)

X<-cbind(rep(1,n),Employed)

solve(t(X)%*%X)%*%t(X)%*%GNP

sigma.hat<-sqrt(sum(residuals(mod1)^2/(n-2)))

sigma.betas<-sigma.hat^2*solve(t(X)%*%X)

sqrt(sigma.betas[1,1])

sqrt(sigma.betas[2,2])

2*pt(coef(mod1)[1]/sqrt(sigma.betas[1,1]),df=n-2)

2*(1-pt(coef(mod1)[2]/sqrt(sigma.betas[2,2]),df=n-2))

grid.x<-as.data.frame(seq(min(Employed),max(Employed),length=100))

names(grid.x)<-"Employed"

pred<-predict(mod1, newdata=grid.x, se.fit = TRUE)

pred

pred$fit[1]

grid.x$Employed[1]

coef(mod1)[1]+coef(mod1)[2]*grid.x$Employed[1]

pred$se.fit[1]

sigma.hat*sqrt(1/n+(grid.x$Employed[1]-mean(Employed))^2/((n-1)*var(Employed)))

lower<-pred$fit+qt(0.005,df=n-2)*pred$se.fit

upper<-pred$fit-qt(0.005,df=n-2)*pred$se.fit

lines(grid.x$Employed,lower,col="red")

lines(grid.x$Employed,upper,col="red")

detach(dati1)

data(LifeCycleSavings)

dati2<-LifeCycleSavings

summary(dati2)

names(dati2)

attach(dati2)

plot(dpi,sr)

mod2<-lm(sr~dpi)

abline(mod2)

summary(mod2)

grid.xx<-as.data.frame(seq(min(dpi),max(dpi),length=100))

names(grid.xx)<-"dpi"

pred<-predict(mod2, newdata=grid.xx, se.fit = TRUE)


PAGINE

2

PESO

33.18 KB

AUTORE

Atreyu

PUBBLICATO

+1 anno fa


DETTAGLI
Esame: Econometria
Corso di laurea: Corso di laurea magistrale in politiche pubbliche
SSD:
A.A.: 2011-2012

I contenuti di questa pagina costituiscono rielaborazioni personali del Publisher Atreyu di informazioni apprese con la frequenza delle lezioni di Econometria e studio autonomo di eventuali libri di riferimento in preparazione dell'esame finale o della tesi. Non devono intendersi come materiale ufficiale dell'università Roma Tre - Uniroma3 o del prof Lagona Francesco.

Acquista con carta o conto PayPal

Scarica il file tutte le volte che vuoi

Paga con un conto PayPal per usufruire della garanzia Soddisfatto o rimborsato

Recensioni
Ti è piaciuto questo appunto? Valutalo!

Altri appunti di Econometria

Esercitazioni con R
Dispensa
Intervalli di confidenza
Dispensa
Richiami sulla distribuzione normale
Dispensa
R - Esercitazioni II parte
Dispensa