Dear All,
#I have the following data:
primate<-log(c(100,95,90,85,80,75,60,40,20,10))
rank<-log(1:10)
#from which I can compute a linear regression, and obtain the confidence
interval for the slope coefficient as follows:
res<-lm(primate~rank)
confint(res,level=0.9)[2,]
#Now, what I would like to obtain is, for a given sequence of "possible" slope
coefficients, their probability of occurrence (not sure #how to call this),
based on the t-distribution.
#Thus for instance if I'm looking at the following sequence
seq(-1,0,0.01)
#I would like to have a range of probability values with the same length, which
sum() should be approximately equal to 1.
#Ultimately what I would like to obtain, is, for any given value q, the
probability that this is the correct estimate for the slope coefficient.
Many Thanks, and apologies if the question sounds stupid!
Enrico
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.