From: Jose Claudio Faria <joseclaudio.faria_at_terra.com.br>

Date: Sat 05 Aug 2006 - 06:09:23 EST

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. Received on Sat Aug 05 06:20:49 2006

Date: Sat 05 Aug 2006 - 06:09:23 EST

I would like to illustrate the origin of the Student t distribution using R.

#== begin options=====

# parameters

nchist = 150

#== end options=======

t = numeric()

pop = rnorm(10000, mean = mu, sd = sigma)

for (i in 1:nsim) {

amo.i = sample(pop, n, replace = TRUE) t[i] = (mean(amo.i) - mu) / (sigma / sqrt(n)) }

win.graph(w = 5, h = 7)

split.screen(c(2,1))

screen(1)

hist(t,

main = "histogram", breaks = nchist, col = "lightgray", xlab = '', ylab = "Fi", font.lab = 2, font = 2)

screen(2)

hist(t,

probability = T, main = 'f.d.p and histogram', breaks = nchist, col = 'lightgray', xlab = 't', ylab = 'f(t)', font.lab = 2, font = 2)

x = t

curve(dt(x, df = n-1), add = T, col = "red", lwd = 2)

Many thanks for any help,

___

Jose Claudio Faria

Brasil/Bahia/Ilheus/UESC/DCET

Estatística Experimental/Prof. Adjunto

mails: joseclaudio.faria@terra.com.br

jc_faria@uesc.br jc_faria@uol.com.br ______________________________________________R-help@stat.math.ethz.ch 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. Received on Sat Aug 05 06:20:49 2006

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.1.8, at Sat 05 Aug 2006 - 08:18:11 EST.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*