Hello,
suppose one has the following values 
x1 <- rnorm(10,5,1)
x2 <- rgamma(10,5,1)
y <- rnorm(10,4,1)
mydat <- data.frame(y,x1,x2)
then one can use glm like
mod <- glm(y~x1+x2, data=mydat, family=gaussian)
  But how could I estimate y_hat?
  Thanks alot!
  Sam

       
---------------------------------

	[[alternative HTML version deleted]]

