Displaying 1 result from an estimated 1 matches for "fueldata".
2008 Mar 02
3
Need help to locate my mistake
Dear readers
I would like to make General Linear Model (GLM) for the following data set
http://louise.hoffman.googlepages.com/fuel.csv
The code I have written is
fuelData<-read.table('fuel.csv',header=TRUE, sep=',')
n<-dim(fuelData)[1]
xOnes<- matrix(1,nrow=n,ncol=1)
x<-cbind(xOnes,fuelData[,3])
y<-fuelData[,4]
theta<-((t(x)%*%x)^(-1))%*%t(x)%*%y
which gives
> theta
[,1]
[1,] 215.8374077
[2,] 0.1083491
When I do...