<p>O problema não é do sample, mas da sua implementação. Vc cria um obj com 30 linhas (resspl) e tenta extrair a 31a linha (veja q eventualmente i=31).</p>
<p><blockquote type="cite">On 25 Apr 2011 10:18, "Alexandre dos Santos" <<a href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a>> wrote:<br><br>
<div lang="PT-BR" link="blue" vlink="purple">
<div>
<p class="MsoNormal"> Bom Dia Pessoal,</p>
<p class="MsoNormal"> Tenho um banco de dados com três variáveis parc
(número da parcela), mont2(número de insetos), dead2 (mortalidade), com 108
dados, estou tentando fazer uma amostragem com a função sample(), para um n da
amostra igual a 30, ate ai tudo bem, mas quando repito a função sample() para
um número de vezes superior ao n da amostra começam a aparecer NA nos sorteio,
porque isso poderia estar acontecendo, uma vez que estou utilizando a função
com reposição, fiz:</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">> nde<-read.table("Tout_NDE.txt",header=T)#
Todos os dados</p>
<p class="MsoNormal">> #</p>
<p class="MsoNormal">>
dead2<-nde$Mortalite[nde$Temps==2]-nde$Mortalite[nde$Temps==1]# Mortalidade
no segundo mês</p>
<p class="MsoNormal">>
mont2<-nde$Monticules[nde$Temps==2]-nde$Monticules[nde$Temps==1]# Número de
insetos no segundo mês</p>
<p class="MsoNormal">> parc<-1:108 # Número da parcela no campo</p>
<p class="MsoNormal"><span lang="EN-US">> mois2<-cbind(parc,mont2,dead2)</span></p>
<p class="MsoNormal">> #</p>
<p class="MsoNormal">> nparc=30# Tamanho da amostra independente a ser
sorteada com sample ()</p>
<p class="MsoNormal">> Nsim=31# Numero de simulações da função sample</p>
<p class="MsoNormal"><span lang="EN-US">> </span></p>
<p class="MsoNormal"><span lang="EN-US">> RES=NULL</span></p>
<p class="MsoNormal"><span lang="EN-US">> </span></p>
<p class="MsoNormal"><span lang="EN-US">> for(i in 1:Nsim){</span></p>
<p class="MsoNormal"><span lang="EN-US">+ </span></p>
<p class="MsoNormal"><span lang="EN-US">+ res=NULL</span></p>
<p class="MsoNormal"><span lang="EN-US">+ spl=sample(parc,nparc, replace = TRUE)</span></p>
<p class="MsoNormal"><span lang="EN-US">+ resspl=mois2[spl[1:30],]</span></p>
<p class="MsoNormal"><span lang="EN-US">+ my.mont=mean(resspl[,2][i])</span></p>
<p class="MsoNormal"><span lang="EN-US">+ my.dead=mean(resspl[,3][i])</span></p>
<p class="MsoNormal"><span lang="EN-US">+ result=rbind(res,c(my.mont,my.dead,i))</span></p>
<p class="MsoNormal"><span lang="EN-US">+ #print(result)</span></p>
<p class="MsoNormal"><span lang="EN-US">+ RES=rbind(RES,result)</span></p>
<p class="MsoNormal"><span lang="EN-US">+
colnames(RES)=c("monticule","mortalite","simulation")</span></p>
<p class="MsoNormal"><span lang="EN-US">+ }</span></p>
<p class="MsoNormal"><span lang="EN-US">> print(RES)</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span>monticule mortalite simulation</p>
<p class="MsoNormal"> [1,] 5 4.000000 1</p>
<p class="MsoNormal"> [2,] 0 0.000000 2</p>
<p class="MsoNormal"> [3,] 1 12.000000 3</p>
<p class="MsoNormal"> [4,] 1 5.263158 4</p>
<p class="MsoNormal"> [5,] 0 0.000000 5</p>
<p class="MsoNormal"> [6,] 0 0.000000 6</p>
<p class="MsoNormal"> [7,] 2 9.090909 7</p>
<p class="MsoNormal"> [8,] 0 0.000000 8</p>
<p class="MsoNormal"> [9,] 1 0.000000 9</p>
<p class="MsoNormal">[10,] 0 0.000000 10</p>
<p class="MsoNormal">[11,] 1 11.764706 11</p>
<p class="MsoNormal">[12,] 3 6.451613 12</p>
<p class="MsoNormal">[13,] 3 12.000000 13</p>
<p class="MsoNormal">[14,] 3 4.545455 14</p>
<p class="MsoNormal">[15,] 1 3.225806 15</p>
<p class="MsoNormal">[16,] 0 0.000000 16</p>
<p class="MsoNormal">[17,] 3 4.545455 17</p>
<p class="MsoNormal">[18,] 1 0.000000 18</p>
<p class="MsoNormal">[19,] 5 13.636364 19</p>
<p class="MsoNormal">[20,] 2 12.121212 20</p>
<p class="MsoNormal">[21,] 2 6.896552 21</p>
<p class="MsoNormal">[22,] 1 5.263158 22</p>
<p class="MsoNormal">[23,] 2 0.000000 23</p>
<p class="MsoNormal">[24,] 10 8.333333 24</p>
<p class="MsoNormal">[25,] 5 14.285714 25</p>
<p class="MsoNormal">[26,] 8 10.000000 26</p>
<p class="MsoNormal">[27,] 8 10.000000 27</p>
<p class="MsoNormal">[28,] 3 4.761905 28</p>
<p class="MsoNormal">[29,] 0 0.000000 29</p>
<p class="MsoNormal">[30,] 1 3.571429 30</p>
<p class="MsoNormal">[31,] NA NA 31</p>
<p class="MsoNormal">> </p>
<p class="MsoNormal">Obrigado,</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"><b>Alexandre dos Santos</b></p>
<p class="MsoNormal"><b>Ingenieur forestier, Msc.</b></p>
<p class="MsoNormal"><b>INRA- Biostatistique et Processus Spatiaux (BioSP)</b></p>
<p class="MsoNormal"><b><span lang="EN-US">Domaine Saint-Paul<br>
Site Agroparc <br>
84914 - Avignon - France<br>
Tél. : +33 (0)6 87 95 16 29</span></b></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
</div>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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>
<br></blockquote></p>