Displaying 2 results from an estimated 2 matches for "windprot".
2004 May 24
2
Manova and specifying the model
...a() funciton and the summary.manova() function to get the appropriate summary of test statistics.
I just don't manage to specify my model in the manova() call. How to specify a model with multiple responses and one explanatory factor?
If I type:
pcor.manova<-manova(isol+hcom+habarea+inclin+windprot+shrubcov+herbh+baregr+flowcov~pcor, data=pcor.df)
I always get error messages like:
Error in manova(isol + hcom + habarea + inclin + windprot + shrubcov + :
need multiple response
Can someone help me?
Thanks
Steffi
----------------------------------------
Stefanie von Felten
Käppeli...
2004 May 24
1
discriminant analysis
...data. With the CV=True option I was not able to perform the predict() call. What do I have to do? Or is there no possibility at all? You also need the predicted values to produce a plot of the analysis, as far as I know.
Here my code:
pcor.lda2<-lda(pcor~habarea+hcom+isol+flowcov+herbh+inclin+windprot+shrubcov+baregr, data=pcor.df, CV=T)
table2<-table(pcor.df$pcor, pcor.lda2$class)
table2
#doesn't work, becoause CV=True?
pcor.ld2<-predict(pcor.lda2, dimen=1)$x
plot(pcor.ld2)
plot(pcor.lda2, type="density", dimen=1) #kernel density estimates
I am happy if I get an answer fr...