Simulates values of expected response for linear regression model
Simulates draws of the posterior distribution of an expected response for a linear regression model with a noninformative prior
blinregexpected(X1,theta.sample)
X1 |
matrix where each row corresponds to a covariate set |
theta.sample |
list with components beta, matrix of simulated draws of regression vector, and sigma, vector of simulated draws of sampling error standard deviation |
matrix where a column corresponds to the simulated draws of the expected response for a given covariate set
Jim Albert
chirps=c(20,16.0,19.8,18.4,17.1,15.5,14.7,17.1,15.4,16.2,15,17.2,16,17,14.1) temp=c(88.6,71.6,93.3,84.3,80.6,75.2,69.7,82,69.4,83.3,78.6,82.6,80.6,83.5,76.3) X=cbind(1,chirps) m=1000 theta.sample=blinreg(temp,X,m) covset1=c(1,15) covset2=c(1,20) X1=rbind(covset1,covset2) blinregexpected(X1,theta.sample)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.