Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças

Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

Link para esta página de comparações

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 09:44]
ehlers
disciplinas:ce223:comandos2008 [2008/05/26 14:57]
ehlers
Linha 328: Linha 328:
 [2,]   ​15 ​  ​30 ​  29 [2,]   ​15 ​  ​30 ​  29
  
-# Cada matrix 2x3 contem as contagens por sexo (linhas) e estado (colunas). A ultima dimensao refere-se ao habito de fumar.+# Cada matrix 2x3 contem as contagens por sexo (linhas) e estado (colunas). ​ 
 +A ultima dimensao refere-se ao habito de fumar.
  
 nomes = list(c('​M','​F'​),​c('​PR','​SC','​RS'​),​c('​nao fuma','​fuma pouco','​fuma muito'​)) nomes = list(c('​M','​F'​),​c('​PR','​SC','​RS'​),​c('​nao fuma','​fuma pouco','​fuma muito'​))
Linha 518: 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 941:
 </​code>​ </​code>​
  
 +=== 30/04/2008 ===
 +
 +Gerando 1000 amostras de tamanho n=20 de uma distribuição normal 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 novamente 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[12pt]{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>​
 +Para usar Unicode, substitua a última linha por
 +<​code>​
 +\usepackage[utf-8]{inputenc}
 +</​code>​
 +
 +Apos o preambulo coloque o titulo, autoria e data do artigo.
 +<​code>​
 +
 +\Title{Título do Trabalho}
 +\author{Nome do Autor}
 +\date{\today}
 +</​code>​
 +
 +Agora sim começa o documento,
 +<​code>​
 +\begin{document}
 +\maketitle
 +% seu texto ...
 +% Inclua seções e subseções ​
 +\section{Uma seção}
 +\subsection{Uma subseção}
 +\end{document}
 +</​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>​
 +\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>​
 +
 +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>​
 +\begin{tabular}{|l|cc|}
 +\hline
 +           &​masculino & feminino \\
 +\hline\hline
 +Não fuma   &​45 ​       & 16       \\
 +Fuma pouco &​28 ​       & 22       \\
 +\hline
 +\end{tabular}
 +</​latex>​
 +
 +Incluindo figuras.
 +
 +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[dvips]{graphicx} no preambulo.
 +<​code>​
 +\begin{figure}[htbp]\centering
 +\includegraphics[width=10cm,​ height=10cm]{histograma.ps}
 +\caption{Histograma de uma amostra de tamanho 1000 da distribuição normal padrão.}
 +\label{fig:​hist}
 +\end{figure}
 +</​code>​
 +
 +Pode-se fazer referencia a figura no texto: figura \ref{hist}, na pagina \pageref{hist}.
  

QR Code
QR Code disciplinas:ce223:comandos2008 (generated for current page)