<div><span style="font-size: medium; color: #0000ff;">Tentei rodar, mas o meu resultado dá muito diferente, qual o motivo?</span></div>
<div> </div>
<div>> dados <- read.table("dados_boot.txt", h=T)</div>
<div>> dados</div>
<div> X4 X5 X15 X6 X11 X12 X13 X16 X17 X18 X19 X20 X25 X26 X27 X29</div>
<div>1 3 4 3 2 3 4 4 2 3 3 3 5 3 4 3 3</div>
<div>2 2 5 5 5 3 5 5 4 5 4 5 4 4 3 5 3</div>
<div>3 4 5 4 3 5 5 5 4 5 6 6 5 4 4 5 4</div>
<div>4 2 6 4 4 2 6 6 4 4 3 6 6 3 4 4 2</div>
<div>5 4 3 1 3 1 4 5 3 3 1 3 1 2 1 2 1</div>
<div>6 5 5 6 5 5 6 6 6 5 5 6 4 5 5 5 4</div>
<div>7 1 7 4 5 7 6 7 5 4 5 5 3 4 4 7 6</div>
<div>8 5 7 6 5 4 6 6 3 7 5 4 3 5 7 5 6</div>
<div>9 5 6 4 3 5 5 5 4 4 4 3 2 2 5 2 5</div>
<div>10 3 6 2 5 5 6 6 6 4 4 5 5 2 2 6 3</div>
<div>11 2 2 3 2 2 4 4 4 3 4 5 4 5 3 4 3</div>
<div>12 2 3 2 2 3 2 5 3 2 4 5 2 4 5 3 2</div>
<div>13 7 7 4 3 3 6 7 6 6 6 6 6 4 4 6 5</div>
<div>14 5 7 6 6 7 6 6 7 7 6 5 5 3 4 6 6</div>
<div>15 3 6 5 6 5 6 6 5 6 6 5 5 4 4 6 6</div>
<div>16 1 2 1 4 1 5 2 2 3 6 5 4 1 2 1 1</div>
<div>17 4 5 2 4 3 4 3 5 5 4 3 3 2 4 3 4</div>
<div>18 5 6 5 5 3 5 5 6 7 6 6 5 4 3 6 6</div>
<div>19 4 3 5 3 1 5 4 2 3 4 4 5 1 4 2 5</div>
<div>20 3 3 3 4 4 3 4 2 2 3 3 3 2 3 1 2</div>
<div>> mean(dados[,1])</div>
<div>[1] 3.5</div>
<div>> library(sem)</div>
<div>> library(boot)</div>
<div>> </div>
<div>> #Bootstrap (médias das replicações)</div>
<div>> </div>
<div>> </div>
<div>> obs <- function(x){</div>
<div>+ </div>
<div>+ results <- c()</div>
<div>+ for(i in 1:500){</div>
<div>+ </div>
<div>+ results[i] <- mean(sample(dados,replace=T))</div>
<div>+ </div>
<div>+ saida <- paste(results,i,sep="")</div>
<div>+ </div>
<div>+ }</div>
<div>+ </div>
<div>+ return(saida)</div>
<div>+ }</div>
<div>> </div>
<div>> resultado<-apply(dados,2,obs)</div>
<div>There were 50 or more warnings (use warnings() to see the first 50)</div>
<div>> </div>
<div>> resultado[,1]</div>
<div> [1] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [10] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [19] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [28] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [37] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [46] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [55] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [64] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [73] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [82] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div> [91] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[100] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[109] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[118] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[127] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[136] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[145] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[154] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[163] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[172] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[181] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[190] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[199] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[208] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[217] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[226] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[235] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[244] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[253] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[262] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[271] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[280] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[289] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[298] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[307] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[316] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[325] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[334] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[343] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[352] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[361] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[370] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[379] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[388] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[397] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[406] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[415] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[424] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[433] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[442] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[451] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[460] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[469] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[478] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[487] "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>[496] "NA500" "NA500" "NA500" "NA500" "NA500"</div>
<div>> </div>
<div> </div>
<div><br />Como eu falei,</div>
<div>
<div> </div>
<div>Basta retirar a variável "saída" da função "obs" que a saída será do tipo numeric e aí funcionará seu comando</div>
<div> </div>
<div>segue função obs , corrigida</div>
<div>
<div> </div>
<div>obs <- function(x){<br /><br /> results <- c()<br /> for(i in 1:500){<br /><br /> results[i] <- mean( sample(x,replace=TRUE))<br /> <br /><br /> }<br /><br /> return(results)<br />}<br /><br /><br />resultado <- apply(dados,2,obs)<br /><br />str(resultado)</div>
<br />
<div>=========================================<br /> Fernando Souza<br /> Zootecnista, DSc. Produção e Alimentação Animal<br /> Celular: (31)99796-8781 (Vivo)<br /> <a href="https://n1.nylas.com/link/a5a0f79f704e8c5d398c43ed1aa06c5eac19980ecdd6befd76ba7e1c8ac49b0b/0?redirect=mailto%3Ae-mail%253Anandodesouza%40gmail.com&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">E-mail:nandodesouza@gmail.com</a><br /> Lattes: <a href="https://n1.nylas.com/link/a5a0f79f704e8c5d398c43ed1aa06c5eac19980ecdd6befd76ba7e1c8ac49b0b/1?redirect=http%3A%2F%2Flattes.cnpq.br%2F6519538815038307&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">http://lattes.cnpq.br/6519538815038307</a><br /> Blog: <a href="https://n1.nylas.com/link/a5a0f79f704e8c5d398c43ed1aa06c5eac19980ecdd6befd76ba7e1c8ac49b0b/2?redirect=https%3A%2F%2Fproducaoanimalcomr.wordpress.com%2F&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">https://producaoanimalcomr.wordpress.com/</a><br /> ==========================================</div>
<div> </div>
Sent from <a href="https://n1.nylas.com/link/a5a0f79f704e8c5d398c43ed1aa06c5eac19980ecdd6befd76ba7e1c8ac49b0b/3?redirect=https%3A%2F%2Fnylas.com%3Fref%3Dn1&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">Nylas Mail</a>, the best free email app for work
<div> </div>
</div>
</div>
<div><img class="n1-open" style="border: 0; width: 0; height: 0;" src="https://n1.nylas.com/open/a5a0f79f704e8c5d398c43ed1aa06c5eac19980ecdd6befd76ba7e1c8ac49b0b?r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" alt="" width="0" height="0" /></div>
<div class="gmail_quote nylas-quote nylas-quote-id-ff6e63ab34b009b9707ed2590a0d6bed3dedf7f11bd492bfe12e69b7a529cc60"><br /> On Mai 18 2017, at 9:44 am, Fernando Antonio de souza <nandodesouza@gmail.com> wrote: <br />
<blockquote class="gmail_quote" style="margin: 0 0 0 .8ex; border-left: 1px #ccc solid; padding-left: 1ex;">
<div dir="auto">
<div>Clodoaldo
<div dir="auto"> </div>
<div dir="auto">Estou sem acesso a meu computador aqui é por as não posso testar. Tente rodar a função obs se a variável saida. Acredito que ela está transformando a saída em carácter e por isso a saída.</div>
<div dir="auto"> </div>
<div dir="auto">Assim q acessar meu PC testo e lhe confirmo. </div>
<div dir="auto"> </div>
<div style="font-family: sans-serif;" dir="auto">obs <- function(x){<br /> <br /> results <- c()<br /> for(i in 1:500){<br /> </div>
<span style="font-family: sans-serif;"> results[i] <- mean(sample(dados,replace=T))</span>
<div style="font-family: sans-serif;" dir="auto"> <br /> }<br /> <br />return(results)<br />}</div>
<div><br />
<div>Em 18/05/2017 9:24 AM, "Clodoaldo José Figueredo" <<a href="../../../undefined//compose?to=clodoaldo.figueredo@ifc-araquari.edu.br" target="_blank">clodoaldo.figueredo@ifc-araquari.edu.br</a>> escreveu:<br />
<blockquote style="margin: 0 0 0 .8ex; border-left: 1px #ccc solid; padding-left: 1ex;">
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div>
<div>Caro Fernando (e demais colegas)<br /><br /></div>
Gostaria de agradecer sua grande ajuda, mas preciso tirar mais outras dúvidas.</div>
Abaixo está o código completo e a matriz dados que é base para os cálculos.</div>
A saída que inseri é a coluna resultado [,1]<br /><br /></div>
Porque não consigo calcular a média dos valores?</div>
As saídas serão usadas como vetores (individualmente) em um outro procedimento para uso do pacote "sem" (Modelagem de equações estruturais).</div>
Como tira as aspas, ou não tem nada a ver?</div>
<div>Porque aparecem valores 4500 e 3500 no meio dos calculados? (Não tem lógica)</div>
<br />setwd("C:/Dados")<br />dados <- read.table("Dados_Boot_txt.txt", h=T)<br /><br />dados<br />
<pre id="m_795982940808604400gmail-rstudio_console_output" style="font-family: 'lucida console'; font-size: 10pt; outline: medium none; border-width: medium; border-style: none; border-color: currentcolor; word-break: break-all; margin: 0px; white-space: pre-wrap; line-height: 15px; color: #000000; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; text-indent: 0px; text-transform: none; word-spacing: 0px; background-color: #ffffff;"><span style="color: blue; white-space: pre-wrap;">> </span><span style="color: blue;">dados
</span> X4 X5 X15 X6 X11 X12 X13 X16 X17 X18 X19 X20 X25 X26 X27 X29
1 3 4 3 2 3 4 4 2 3 3 3 5 3 4 3 3
2 2 5 5 5 3 5 5 4 5 4 5 4 4 3 5 3
3 4 5 4 3 5 5 5 4 5 6 6 5 4 4 5 4
4 2 6 4 4 2 6 6 4 4 3 6 6 3 4 4 2
5 4 3 1 3 1 4 5 3 3 1 3 1 2 1 2 1
6 5 5 6 5 5 6 6 6 5 5 6 4 5 5 5 4
7 1 7 4 5 7 6 7 5 4 5 5 3 4 4 7 6
8 5 7 6 5 4 6 6 3 7 5 4 3 5 7 5 6
9 5 6 4 3 5 5 5 4 4 4 3 2 2 5 2 5
10 3 6 2 5 5 6 6 6 4 4 5 5 2 2 6 3
11 2 2 3 2 2 4 4 4 3 4 5 4 5 3 4 3
12 2 3 2 2 3 2 5 3 2 4 5 2 4 5 3 2
13 7 7 4 3 3 6 7 6 6 6 6 6 4 4 6 5
14 5 7 6 6 7 6 6 7 7 6 5 5 3 4 6 6
15 3 6 5 6 5 6 6 5 6 6 5 5 4 4 6 6
16 1 2 1 4 1 5 2 2 3 6 5 4 1 2 1 1
17 4 5 2 4 3 4 3 5 5 4 3 3 2 4 3 4
18 5 6 5 5 3 5 5 6 7 6 6 5 4 3 6 6
19 4 3 5 3 1 5 4 2 3 4 4 5 1 4 2 5
20 3 3 3 4 4 3 4 2 2 3 3 3 2 3 1 2</pre>
<br />mean(dados[,1])<br />
<pre id="m_795982940808604400gmail-rstudio_console_output" style="font-family: 'lucida console'; font-size: 10pt; outline: medium none; border-width: medium; border-style: none; border-color: currentcolor; word-break: break-all; margin: 0px; white-space: pre-wrap; line-height: 15px; color: #000000; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; text-indent: 0px; text-transform: none; word-spacing: 0px; background-color: #ffffff;"><span style="color: blue; white-space: pre-wrap;">> </span><span style="color: blue;">mean(dados[,1])
</span>[1] 3.5</pre>
<br />library(sem)<br />library(boot)<br /><br />#Bootstrap (médias das replicações)
<div><br /><br />obs <- function(x){<br /> <br /> results <- c()<br /> for(i in 1:500){<br /> </div>
results[i] <- mean(sample(dados,replace=T))
<div><br /> saida <- paste(results,i,sep="")<br /> <br /> }<br /> <br />return(saida)<br />}<br /><br />resultado<-apply(dados,2,obs)<br /><br /></div>
resultado[,1]<br />
<pre id="m_795982940808604400gmail-rstudio_console_output" style="font-family: 'lucida console'; font-size: 10pt; outline: medium none; border-width: medium; border-style: none; border-color: currentcolor; word-break: break-all; margin: 0px; white-space: pre-wrap; line-height: 15px; color: #000000; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; text-indent: 0px; text-transform: none; word-spacing: 0px; background-color: #ffffff;"><span style="color: blue; white-space: pre-wrap;">> </span><span style="color: blue;">resultado[,1]
</span> [1] "3.75500" "4.1500" "4.05500" "3.5500" "3.8500" "4.1500" "3.45500" "4.75500"
[9] "4.25500" "3.45500" "3.6500" "3.85500" "3.9500" "4.05500" "3.6500" "3.75500"
[17] "4.2500" "4.4500" "3.65500" "4.2500" "4500" "4.15500" "5.05500" "3.75500"
[25] "3.8500" "4.45500" "4.25500" "3.65500" "3.65500" "4.05500" "4.05500" "3.3500"
[33] "3.4500" "4.15500" "3.15500" "3.8500" "3.7500" "4.3500" "4500" "3.15500"
[41] "4.55500" "3.95500" "4.1500" "4.05500" "4.2500" "3.6500" "3.95500" "4.35500"
[49] "3.65500" "4.75500" "3.7500" "3.65500" "3.6500" "4.85500" "3.8500" "3.05500"
[57] "3.85500" "3.8500" "3.8500" "3.35500" "3.7500" "4.25500" "4.25500" "3.8500"
[65] "3.95500" "3.25500" "3.9500" "3.75500" "3.7500" "4.45500" "3.65500" "3.4500"
[73] "4500" "4500" "3.55500" "3.9500" "3.75500" "4.3500" "3.3500" "4.25500"
[81] "4.1500" "4.2500" "4.2500" "4.2500" "3.95500" "4.15500" "4.35500" "4.1500"
[89] "4.1500" "4.05500" "3.9500" "3.65500" "3.4500" "3.65500" "4.15500" "3.35500"
[97] "3.55500" "3.75500" "3.8500" "3.8500" "4.25500" "4.15500" "4500" "3.95500"
[105] "3.95500" "3.65500" "4.4500" "3.6500" "3.75500" "3.55500" "3.9500" "4.3500"
[113] "3.9500" "3.55500" "3.35500" "4.1500" "3.65500" "3.75500" "3.6500" "3.05500"
[121] "4500" "3.9500" "4.05500" "3.65500" "4.15500" "4.25500" "4500" "4.05500"
[129] "4.7500" "3.1500" "3.2500" "3.85500" "3.5500" "3.5500" "4.4500" "3.45500"
[137] "3.7500" "3.65500" "4.1500" "4.05500" "3.35500" "4.1500" "3.9500" "3.1500"
[145] "4.2500" "4.15500" "4.45500" "3.95500" "3.25500" "3.9500" "3.9500" "3.75500"
[153] "3.65500" "3.6500" "3.5500" "3.65500" "3.9500" "4.05500" "4.05500" "3.85500"
[161] "4.05500" "3.45500" "3.05500" "3.45500" "4.3500" "4.6500" "3.8500" "4.35500"
[169] "4.05500" "4.1500" "4.15500" "4.8500" "3.6500" "3.05500" "4500" "3.95500"
[177] "3.6500" "4.4500" "4.25500" "4.35500" "3.6500" "4.2500" "3.9500" "4.55500"
[185] "3.9500" "4.3500" "4.15500" "3.35500" "3.6500" "3.9500" "3.85500" "3.45500"
[193] "3.65500" "3.45500" "3.7500" "3.45500" "4.3500" "3.7500" "3.4500" "3.55500"
[201] "3.75500" "3.85500" "3.8500" "4.65500" "3.8500" "3.6500" "3.7500" "3.45500"
[209] "3.65500" "4.2500" "3.5500" "3.75500" "4.4500" "3.8500" "3.45500" "3.5500"
[217] "4.5500" "4.1500" "3.85500" "3.6500" "3.15500" "4.55500" "3.5500" "3.65500"
[225] "3.8500" "4500" "4.15500" "3.65500" "3.35500" "3.5500" "3.45500" "3.65500"
[233] "3.5500" "4.1500" "3.5500" "4.15500" "4.15500" "4.15500" "4.2500" "4.05500"
[241] "4.15500" "3.8500" "3.85500" "4.05500" "4.1500" "3.75500" "3.75500" "3.9500"
[249] "3.25500" "4.15500" "3.85500" "4.25500" "4500" "3.9500" "4.15500" "4500"
[257] "3.9500" "3.95500" "4500" "4.1500" "3.8500" "4.5500" "4.4500" "3.9500"
[265] "3.95500" "4.4500" "3.65500" "3.5500" "4500" "3.8500" "3.4500" "4.2500"
[273] "3.05500" "3.6500" "3.6500" "4.2500" "3.4500" "3.55500" "4.25500" "3.95500"
[281] "3.4500" "3.9500" "4500" "4.05500" "3.75500" "3.85500" "3.95500" "4.3500"
[289] "3.1500" "4.35500" "3.75500" "3.25500" "3.95500" "3.65500" "3.75500" "4.45500"
[297] "4.2500" "4.2500" "3.8500" "4500" "3.3500" "4.3500" "3.6500" "4.25500"
[305] "3.45500" "3.65500" "5.05500" "4.1500" "3.85500" "4.1500" "3.8500" "3.65500"
[313] "3.55500" "4.05500" "3.45500" "3.8500" "4.15500" "3.3500" "3.85500" "4.45500"
[321] "4.05500" "3.8500" "3.6500" "4.65500" "3.95500" "3.55500" "3.7500" "3.95500"
[329] "3.8500" "3.7500" "3.55500" "3.95500" "3.55500" "3.15500" "4500" "3.35500"
[337] "4500" "4.2500" "3.85500" "3.6500" "3.8500" "3500" "4.2500" "3.55500"
[345] "3.8500" "4.25500" "3.5500" "3.85500" "3.8500" "3.25500" "3.65500" "4.4500"
[353] "4500" "3.6500" "4.4500" "3.3500" "4.45500" "4.25500" "4.05500" "3.45500"
[361] "4.3500" "3.85500" "4.15500" "3.65500" "3.55500" "4.3500" "3.9500" "4.4500"
[369] "4.15500" "4.15500" "3.85500" "4.25500" "3.7500" "3.8500" "4.2500" "3.35500"
[377] "3.2500" "3.8500" "4.2500" "3.95500" "3.85500" "3.25500" "4.35500" "3.1500"
[385] "3.7500" "3.35500" "4.1500" "4.15500" "3.85500" "4.55500" "3.65500" "3.75500"
[393] "3.8500" "3.25500" "3.2500" "3.5500" "3.75500" "3.45500" "4.65500" "4.15500"
[401] "3.25500" "4.1500" "4.3500" "4.2500" "3.4500" "3.9500" "4.6500" "4.05500"
[409] "4.35500" "3.55500" "3.8500" "3.8500" "3.55500" "3.8500" "3.7500" "4.05500"
[417] "3.65500" "3.55500" "3.5500" "4.15500" "3.75500" "3.6500" "3.85500" "3.3500"
[425] "3.85500" "2.95500" "3.85500" "4.45500" "3.85500" "3.45500" "3.65500" "4.1500"
[433] "3.6500" "3.9500" "3.8500" "3.15500" "3.65500" "3.4500" "3.5500" "4.25500"
[441] "4.2500" "4500" "3.6500" "3.6500" "3.25500" "4.25500" "3.95500" "3.1500"
[449] "4.2500" "3.15500" "3.65500" "4500" "3.4500" "4.55500" "4.1500" "4.5500"
[457] "3.75500" "3.85500" "4500" "3.85500" "4.55500" "3.6500" "4.1500" "4.3500"
[465] "4.05500" "3.7500" "3.95500" "3.55500" "3.1500" "3.95500" "3.85500" "3.3500"
[473] "3.95500" "4.25500" "3.95500" "4.1500" "4.05500" "4.25500" "3.75500" "4.3500"
[481] "4.3500" "4.3500" "4.7500" "4.45500" "2.9500" "3.5500" "3.6500" "3.55500"
[489] "4.4500" "4500" "3.4500" "3.5500" "3.45500" "3.35500" "4.65500" "4.3500"
[497] "3.8500" "4.3500" "3.95500" "3.65500"</pre>
<br />mean(resultado[,1])<br />
<pre id="m_795982940808604400gmail-rstudio_console_output" style="font-family: 'lucida console'; font-size: 10pt; outline: medium none; border-width: medium; border-style: none; border-color: currentcolor; word-break: break-all; margin: 0px; white-space: pre-wrap; line-height: 15px; color: #000000; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; text-indent: 0px; text-transform: none; word-spacing: 0px; background-color: #ffffff;"><span style="color: blue; white-space: pre-wrap;">> </span><span style="color: blue;">mean(resultado[,1])
</span>[1] NA
<span style="color: #c5060b;">Warning message:
</span><span style="color: #c5060b;">In mean.default(resultado[, 1]) :</span><span style="color: #c5060b;"> argumento não é numérico nem lógico: retornando NA<br /></span><br />Muito obrigado </pre>
<div>
<div><br />
<div>
<div>
<div dir="ltr">
<div><strong>"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." </strong></div>
<div><a href="http://www.pensador.info/frase/NTQwOTE1/&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">Alexander Graham Bell</a></div>
<div> </div>
<div>Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348</div>
<div>Métodos Numéricos para Engenharia - Matemática Aplicada</div>
<div>Instituto Federal Catarinense - Campus Araquari<br />Rodovia BR 280 - km 27 - Cx. Postal 21<br />CEP 89245-000 - Araquari/SC<br />Fone: (47) 3803-7240</div>
<div> </div>
</div>
</div>
</div>
</div>
<div>
<div>Em 17 de maio de 2017 17:42, Fernando Antonio de souza <span dir="ltr"><<a href="../../../undefined//compose?to=nandodesouza@gmail.com" target="_blank">nandodesouza@gmail.com</a>></span> escreveu:</div>
<blockquote style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid #cccccc; padding-left: 1ex;">
<div>
<div>Clodoaldo</div>
<div> </div>
<div>O código abaixo , ler seus dados e para cada coluna, aplica a função "obs" a saída será uma tabela contendo o número de colunas do seu dado fornecido e cada coluna conterá 500 linhas correspondente as quinhentas amostragens realizadas pela função "sample".</div>
<div> </div>
<div>Assim você pode indentificar cada vetor contendo os 500 dados amostrados através do index de coluna</div>
<div> por exemplo : resultado[,1]</div>
<div> </div>
<div>Espero ter ajudado</div>
<div> </div>
<div>A disposição</div>
<div> </div>
dados <- matrix(nrow=10,ncol=16,data=rnorm(160))<br /><br />obs <- function(x){<span><br /> <br /> results <- c()<br /> for(i in 1:500){<br /> <br /></span> results[i] <- mean( sample(x,replace=TRUE))<br /> saida <- paste(results,i,sep="")<br /> <br /> }<br /> <br /> return(saida)<br />}<br /><br />resultado<-apply(dados,2,obs)<br /><br />
<div>=========================================<br /> Fernando Souza<br /> Zootecnista, DSc. Produção e Alimentação Animal<br /> Celular: (31)99796-8781 (Vivo)<br /> <a href="https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/0?redirect=mailto%3Ae-mail%253Anandodesouza%40gmail.com&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">E-mail:nandodesouza@gmail.com</a><br /> Lattes: <a href="https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/1?redirect=http%3A%2F%2Flattes.cnpq.br%2F6519538815038307&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">http://lattes.cnpq.br/6519538815038307</a><br /> Blog: <a href="https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/2?redirect=https%3A%2F%2Fproducaoanimalcomr.wordpress.com%2F&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">https://producaoanimalcomr.wordpress.com/</a><br /> ==========================================</div>
<div> </div>
Sent from <a href="https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/3?redirect=https%3A%2F%2Fnylas.com%3Fref%3Dn1&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">Nylas Mail</a>, the best free email app for work
<div> </div>
</div>
<img style="border-width: 0px; border-style: none; border-color: currentcolor; width: 0px; height: 0px;" src="https://n1.nylas.com/open/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e?r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy" alt="" width="0" height="0" />
<div>
<div>
<div>
<div><br /> On Mai 17 2017, at 4:19 pm, Clodoaldo José Figueredo via R-br <<a href="../../../undefined//compose?to=r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>> wrote: </div>
</div>
<blockquote style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid #cccccc; padding-left: 1ex;">
<div>
<div>
<div dir="ltr">
<div>Na rotina abaixo eu preciso gravar cada vetor results para os 16 valores de j. Preciso colocar um índice porque esses 16 resultados formarão uma nova matriz com 16 colunas e 500 linhas.</div>
Como indexar os vetores results no comando "for" e montar a nova matriz?<br />
<div>
<div><br />for(j in 1:16){<br />y <- dados[,j]<br />qqnorm(y);qqline(y)<br />mean(y)<br />results <- c()<br />for(i in 1:500){ results[i] <- mean( sample(y,replace=T))}<br />results<br />qqnorm(results);qqline(results)<br />mean (results) <br />}<br /><br /></div>
<div>Obrigado<br /><br /></div>
<div>Clodoaldo<br /><br /></div>
<div>
<div>
<div>
<div dir="ltr">
<div><strong>"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." </strong></div>
<div><a href="http://www.pensador.info/frase/NTQwOTE1/&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">Alexander Graham Bell</a></div>
<div> </div>
<div>Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348</div>
<div>Métodos Numéricos para Engenharia - Matemática Aplicada</div>
<div>Instituto Federal Catarinense - Campus Araquari<br />Rodovia BR 280 - km 27 - Cx. Postal 21<br />CEP 89245-000 - Araquari/SC<br />Fone: (47) 3803-7240</div>
<div> </div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<pre>_______________________________________________
R-br mailing list
<a href="../../../undefined//compose?to=R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia&r=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</pre>
</blockquote>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<div>_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<br /><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br />Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div>