1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
|
\name{blinregexpected}
\alias{blinregexpected}
\title{Simulates values of expected response for linear regression model}
\description{
Simulates draws of the posterior distribution of an expected response for
a linear regression model with a noninformative prior}
\usage{
blinregexpected(X1,theta.sample)
}
\arguments{
\item{X1}{matrix where each row corresponds to a covariate set}
\item{theta.sample}{list with components beta, matrix of simulated draws of regression vector,
and sigma, vector of simulated draws of sampling error standard deviation}
}
\value{
matrix where a column corresponds to the simulated draws of the expected response for a given
covariate set
}
\author{Jim Albert}
\examples{
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)
}
\keyword{models}
|