<div dir="ltr"><div>Olá, Alexandre</div><div><br></div><div>Tente criar um novo data.frame para o resultado e iterar pelos valores de Amb:</div><div><br></div><div>Repare que o parâmetro "na.rm = T" do comando sd irá desconsiderar os registros com rate = NA no cálculo do desvio padrão.<br></div><div><br></div><div>d$Amb = factor(d$Amb)<br>ambs = levels(d$Amb)<br><br>resultado = data.frame(<br> State = character()<br> , Amb = character()<br> , Year = integer()<br> , rate = double()<br> , desv_pad = double())<br><br># quantas vezes irá realizar o experimento:<br># "sortear 5 linhas dentro de cada variável Amb e calcular o desvio padrão da variável rate"<br>repeticoes = 2 <br><br>for(amb in ambs) {<br> subpopulacao = d[d$Amb == amb]<br> for (i in 1:repeticoes) {<br> amostra = subpopulacao[sample(NROW(subpopulacao), 5),]<br> amostra$desv_pad = sd(amostra$rate, na.rm = T)<br> resultado = rbind(resultado, amostra)<br> }<br>}</div><div><br></div><div>Atenciosamente,</div><div>Rodrigo Ângelo<br></div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Nov 5, 2018 at 2:38 PM ASANTOS por (R-br) <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Prezados Membros,<br>
<br>
Segue um exemplo do que estou tentando fazer:<br>
<br>
#Tenho um banco de dados d com 4 variáveis (State, Amb, Year e rate)<br>
<br>
library(data.table)<br>
d <- fread ("<br>
State Amb Year rate<br>
GA soil 2015 .4<br>
GA soil 2015 .4<br>
GA soil 2015 .5<br>
GA soil 2015 .2<br>
GA soil 2015 .3<br>
GA soil 2015 NA<br>
GA soil 2016 1.0<br>
GA soil 2016 2.0<br>
GA soil 2016 .5<br>
GA soil 2016 1.0<br>
GA veg 2015 .4<br>
GA veg 2015 .4<br>
GA veg 2015 .5<br>
GA veg 2015 .2<br>
GA veg 2015 .3<br>
GA veg 2015 NA<br>
GA veg 2016 1.0<br>
GA veg 2016 2.0<br>
GA veg 2016 .5<br>
GA veg 2016 1.0<br>
")<br>
#<br>
<br>
Gostaria de sortear 5 linhas dentro de cada variável Amb, calcular o <br>
desvio padrão da variável rate e repetir<br>
os valores obtidos nas linhas envolvidas em uma nova coluna (desv_pad) <br>
ao fim cada sorteio, sendo meu output desejado em um sorteio hipotético:<br>
<br>
State Age Year rate desv_pad<br>
1: GA soil 2015 0.4 0.7549834<br>
2: GA soil 2015 0.4 0.7549834<br>
3: GA soil 2015 0.5 0.3082207<br>
4: GA soil 2015 0.2 0.3082207<br>
5: GA soil 2015 0.3 0.3082207<br>
6: GA soil 2015 NA 0.7549834<br>
7: GA soil 2016 1.0 0.7549834<br>
8: GA soil 2016 2.0 0.7549834<br>
9: GA soil 2016 0.5 0.3082207<br>
10: GA soil 2016 1.0 0.3082207<br>
11: GA veg 2015 0.4 0.580517<br>
12: GA veg 2015 0.8 0.8736895<br>
13: GA veg 2015 0.5 0.580517<br>
14: GA veg 2015 1.5 0.580517<br>
15: GA veg 2015 0.3 0.8736895<br>
16: GA veg 2015 NA 0.8736895<br>
17: GA veg 2016 1.0 0.580517<br>
18: GA veg 2016 NA 0.8736895<br>
19: GA veg 2016 2.0 0.8736895<br>
20: GA veg 2016 1.7 0.580517<br>
<br>
Alguém teria uma dica de como resolver isso?<br>
<br>
Obrigado<br>
<br>
-- <br>
======================================================================<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
Fone: <a href="tel:+55%2065%2099686-6970" value="+5565996866970" target="_blank">(+55) 65 99686-6970</a> (VIVO) <a href="tel:+55%2065%203221-2674" value="+556532212674" target="_blank">(+55) 65 3221-2674</a> (FIXO)<br>
<a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a><br>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a><br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" rel="noreferrer" target="_blank">http://lattes.cnpq.br/1360403201088680</a><br>
OrcID: <a href="http://orcid.org/0000-0001-8232-6722" rel="noreferrer" target="_blank">orcid.org/0000-0001-8232-6722</a> - ResearcherID: A-5790-2016<br>
Researchgate: <a href="http://www.researchgate.net/profile/Alexandre_Santos10" rel="noreferrer" target="_blank">www.researchgate.net/profile/Alexandre_Santos10</a><br>
LinkedIn: <a href="http://br.linkedin.com/in/alexandre-dos-santos-87961635" rel="noreferrer" target="_blank">br.linkedin.com/in/alexandre-dos-santos-87961635</a><br>
Mendeley:<a href="http://www.mendeley.com/profiles/alexandre-dos-santos6/" rel="noreferrer" target="_blank">www.mendeley.com/profiles/alexandre-dos-santos6/</a><br>
======================================================================<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e fornea cdigo mnimo reproduzvel.</blockquote></div>