# t=(media.amostral-media.pop)/(dp.amostral/sqrt(n)) ~ t-Student n=25 xbar=167 s=15 media.pop=160 t=(xbar-media.pop)/(s/sqrt(n)) tt <- seq(-5, 5,l=200) fx<-curve(dt(tt, n-1), min(tt), max(tt), ylab="f(t)",xname = "tt",xlab='t') legend(2,.4,c('t(24)'),bty='n') x <- seq(t,5, l=100) fx <- dt(x, n-1) polygon(c(x, rev(x)),c(fx, rep(0, length(fx))),col="gray70") pp=1-pt(t,n-1) pp=round(pp,3) text(3,0.05,pp) n=25 xbar=153 s=15 media.pop=160 t=(xbar-media.pop)/(s/sqrt(n)) tt <- seq(-5, 5,l=200) fx<-curve(dt(tt, n-1), min(tt), max(tt), ylab="f(t)",xname = "tt",xlab='t') legend(2,.4,c('t(24)'),bty='n') x <- seq(-5,t, l=100) fx <- dt(x, n-1) polygon(c(x, rev(x)),c(fx, rep(0, length(fx))),col="gray70") pp=pt(t,n-1) pp=round(pp,3) text(-3,0.05,pp)