Not sure why that is.  You can always use manova() to get the multivariate 
summary statistics and then use individual models to get predictions for 
each response.  For example,

mydata <- data.frame(y1=rnorm(50), y2=rnorm(50), x1=rnorm(50), 
x2=rnorm(50), x3=rnorm(50))

myfit <- manova(cbind(y1, y2) ~ x1 + x2 + x3, data=mydata)
summary(myfit)

myfit1 <- lm(y1 ~ x1 + x2 + x3, data=mydata)
myfit2 <- lm(y1 ~ x1 + x2 + x3, data=mydata)
predict(myfit1, se.fit=TRUE)
predict(myfit2, se.fit=TRUE)

Jean


BrutishFruit <brutishfruit@hotmail.com> wrote on 08/25/2012 04:09:04 PM:
> 
> Hi,
> I have problem getting the standard deviation from the manova output.
> 
> I have used the manova function:          myfit <- manova(cbind(y1, y2) 
~ x1
> + x2 + x3, data=mydata) .
> I tried to get the predicted values and their standard deviation by 
> using: 
> predict(myfit, type="response", se.fit=TRUE)
> 
> But the problem is that I don't get the standard deviation values, I 
only
> get the predicted values for y1 and y2.
> 
> But if I type:     predict*.lm*(myfit, type="response", se.fit=TRUE)
> I get the predicted values and standard deviation, but only for y1 (and
> nothing from y2...).
> 
> //BF

	[[alternative HTML version deleted]]


