Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anterior Revisão anterior Próxima revisão | Revisão anterior Próxima revisão Ambos lados da revisão seguinte | ||
disciplinas:ce223:comandos2008 [2008/04/23 15:41] ehlers |
disciplinas:ce223:comandos2008 [2008/05/07 16:47] ehlers |
||
---|---|---|---|
Linha 519: | Linha 519: | ||
by(d4[,4],d4$sexo,function(x)as.character(x)) | by(d4[,4],d4$sexo,function(x)as.character(x)) | ||
</code> | </code> | ||
+ | |||
+ | |||
Listas | Listas | ||
Linha 938: | Linha 940: | ||
6 4 2 4 | 6 4 2 4 | ||
</code> | </code> | ||
+ | |||
+ | === 30/04/2008 === | ||
+ | |||
+ | Gerando 1000 amostras de tamanho n=20 de uma ormal padrão | ||
+ | <code R> | ||
+ | rnorm(20, m=70, sd=10) | ||
+ | ams <- matrix(rnorm(20*1000, m=70, sd=10), ncol=20) | ||
+ | dim(ams) | ||
+ | ams[1,] | ||
+ | ams[2,] | ||
+ | </code> | ||
+ | |||
+ | Calculando o valor da estatística de interesse para a primeira e segunda amostra | ||
+ | <code R> | ||
+ | max(ams[1,])/quantile(ams[1,], prob=0.75) | ||
+ | unname(max(ams[1,])/quantile(ams[1,], prob=0.75)) | ||
+ | unname(max(ams[2,])/quantile(ams[2,], prob=0.75)) | ||
+ | </code> | ||
+ | |||
+ | Escrevendo uma função que calcula o valor da estatística de interesse e calculando ovamente o valor para a primeira e segunda amostras. | ||
+ | </code R> | ||
+ | T.est <- function(x) unname(max(x)/quantile(x, prob=0.75)) | ||
+ | T.est(ams[1,]) | ||
+ | T.est(ams[2,]) | ||
+ | </code> | ||
+ | |||
+ | Calculando valor da estatística de interesse agora para todas as amostras de uma só vez | ||
+ | <code R> | ||
+ | ts <- apply(ams, 1, T.est) | ||
+ | length(ts) | ||
+ | ts | ||
+ | </code> | ||
+ | |||
+ | Explorando os resultados: medidas resumo, grafico de densidade estimada e IC (95%) | ||
+ | <code R> | ||
+ | summary(ts) | ||
+ | plot(density(ts)) | ||
+ | quantile(ts, prob=c(0.025, 0.975)) | ||
+ | </code> | ||
+ | |||
+ | Aumentando o número de amostras para 5000. | ||
+ | <code R> | ||
+ | ams <- matrix(rnorm(20*5000, m=70, sd=10), ncol=20) | ||
+ | ts <- apply(ams, 1, T.est) | ||
+ | plot(density(ts)) | ||
+ | </code> | ||
+ | |||
+ | Distribuição amostral da média: empírica (por simulação) //versus// teórica | ||
+ | <code R> | ||
+ | medias <- apply(ams, 1, mean) | ||
+ | plot(density(medias)) | ||
+ | curve(dnorm(x,mean=70, sd=10/sqrt(20)), 60, 80, add=TRUE, col=2) | ||
+ | </code> | ||
+ | |||
+ | === 07/05/2008 === | ||
+ | |||
+ | Exercicios sobre o uso do Latex. | ||
+ | |||
+ | Um preambulo basico: | ||
+ | <code> | ||
+ | \documentclass[A4paper,12pt,oneside]{article}% classes basicas: book, article, report e letter | ||
+ | \usepackage[brazil]{babel}% português do Brasil. | ||
+ | \usepackage[latin1]{inputenc}% usar o conjunto de caracteres Europeu Ocidental. | ||
+ | </code> | ||
+ | |||
+ | Escreva os comandos Latex para as seguintes formulas matematicas: | ||
+ | * <latex>$E(X^2)=\int_{-\infty}^{\infty}x^2 f(x)dx$</latex> | ||
+ | * <latex>$X\sim N(\mu,\sigma^2)\rightarrow f(x)=\frac{1}{\sqrt{2\pi\sigma^2}}\exp\left\{-\frac{1}{2\sigma^2}(x-\mu)^2\right\}$</latex> | ||
+ | * <latex>$(a+b)^n = \sum_{k=0}^n \frac{n!}{k!(n-k)!} a^k b^{n-k}</latex> | ||
+ | |||
+ | Escreva os comandos do Latex para construir a seguinte matriz: | ||
<latex> | <latex> | ||
+ | \left[ | ||
+ | \begin{array}{cccc} | ||
+ | a_{11} & a_{12} & \dots & a_{1n}\\ | ||
+ | a_{21} & a_{22} & \dots & a_{2n}\\ | ||
+ | \vdots & \vdots & \vdots & \vdots\\ | ||
+ | a_{m1} & a_{m2} & \dots & a_{mn}\\ | ||
+ | \end{array} | ||
+ | \right] | ||
+ | </latex> | ||
- | $a^2$ | + | Escreva os comandos do Latex para montar as seguintes tabelas. Tente numerar as tabelas, colocar um comentário e referenciar as tabelas no texto. |
+ | <latex> | ||
+ | \begin{tabular}{ccc} | ||
+ | \hline | ||
+ | &masculino & feminino \\ | ||
+ | \hline | ||
+ | Não fumac &45 & 16 \\ | ||
+ | Fuma pouco &28 & 22 \\ | ||
+ | \hline | ||
+ | \end{tabular} | ||
</latex> | </latex> | ||
+ | |||
+ | <latex> | ||
+ | \begin{tabular}{|l|cc|} | ||
+ | \hline | ||
+ | &masculino & feminino \\ | ||
+ | \hline\hline | ||
+ | Não fuma &45 & 16 \\ | ||
+ | Fuma pouco &28 & 22 \\ | ||
+ | \hline | ||
+ | \end{tabular} | ||
+ | </latex> | ||
+ | |||
+ | Os comandos abaixo criam um arquivo Postscript com um histograma. | ||
+ | <code R> | ||
+ | postscript('histograma.ps') | ||
+ | hist(rnorm(1000)) | ||
+ | dev.off() | ||
+ | </code> | ||
+ | |||
+ | Use os comandos abaixo para incluir a figura no seu documento. Será preciso incluir o comando \usepackage{graphicx} no preambulo. | ||
+ | <code> | ||
+ | \begin{figure}[h] | ||
+ | \center | ||
+ | \includegraphics[width=10cm]{histograma.ps} | ||
+ | \caption{Histograma de uma amostra de tamanho 1000 da distribuição normal padrão.} | ||
+ | \label{hist} | ||
+ | \end{figure} | ||
+ | </code> | ||
+ | |||
+ | Pode-se fazer referencia a figura no texto: figura \ref{hist}, na pagina \pageref{hist}. | ||