É o que o Benilton disse, tu ta tentando pegar o 31º vetor de algo com tamanho 30. O que por sinal tava errado, porque tu tava pegando somente a i-ésima observação da i-ésima simulação, e não a média da simulação.<div><br>
</div><div>Creio que isso aqui faz o que tu quer (em negrito o que eu mudei). <div><br></div><div><div>for(i in 1:Nsim){</div><div> res=NULL</div><div> spl=sample(parc,nparc, replace = TRUE)</div> resspl=mois2[spl[1:30],]<div>
<b> my.mont=mean(resspl[,2])</b></div><div><b> my.dead=mean(resspl[,3])</b></div><div> result=rbind(res,c(my.mont,my.dead,i))</div><div> #print(result)</div><div> RES=rbind(RES,result)</div><div> colnames(RES)=c("monticule","mortalite","simulation")</div>
<div>}</div><div><br></div><div><br></div><br><div class="gmail_quote">2011/4/25 Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
De um exemplo reproduzivel e tenho certeza que muitos terao prazer em<br>
ajuda-lo. b<br>
<br>
2011/4/25 Alexandre dos Santos <<a href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a>>:<br>
<div><div></div><div class="h5">> Obrigado Benilton,<br>
><br>
><br>
><br>
>      Agora já sei onde esta o problema, mais é que eu preciso do valor das<br>
> duas variáveis  na segunda e terceira colunas associadas aos valores<br>
> sorteados para a variável parc, tentei fazer resspl=mois2[spl,], também sem<br>
> sucesso, pois meu objetivo é fazer um sorteio de n=30 e repetir 99 vezes,<br>
><br>
> Alexandre<br>
><br>
><br>
><br>
><br>
><br>
> From: <a href="mailto:r-br-bounces@listas.c3sl.ufpr.br">r-br-bounces@listas.c3sl.ufpr.br</a><br>
> [mailto:<a href="mailto:r-br-bounces@listas.c3sl.ufpr.br">r-br-bounces@listas.c3sl.ufpr.br</a>] On Behalf Of Benilton Carvalho<br>
> Sent: segunda-feira, 25 de abril de 2011 08:15<br>
> To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
> Subject: Re: [R-br] Problema com sample()<br>
><br>
><br>
><br>
> O problema não é do sample, mas da sua implementação. Vc cria um obj com 30<br>
> linhas (resspl) e tenta extrair a 31a linha  (veja q eventualmente i=31).<br>
><br>
> On 25 Apr 2011 10:18, "Alexandre dos Santos"<br>
> <<a href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a>> wrote:<br>
><br>
>     Bom Dia Pessoal,<br>
><br>
>          Tenho um banco de dados com três variáveis parc (número da<br>
> parcela), mont2(número de insetos), dead2 (mortalidade), com 108 dados,<br>
> estou tentando fazer uma amostragem com a função sample(), para um n da<br>
> amostra igual a 30, ate ai tudo bem, mas quando repito a função sample()<br>
> para um número de vezes superior ao n da amostra começam a aparecer NA nos<br>
> sorteio, porque isso poderia estar acontecendo, uma vez que estou utilizando<br>
> a função com reposição, fiz:<br>
><br>
><br>
><br>
>> nde<-read.table("Tout_NDE.txt",header=T)# Todos os dados<br>
><br>
>> #<br>
><br>
>> dead2<-nde$Mortalite[nde$Temps==2]-nde$Mortalite[nde$Temps==1]#<br>
>> Mortalidade no segundo mês<br>
><br>
>> mont2<-nde$Monticules[nde$Temps==2]-nde$Monticules[nde$Temps==1]# Número<br>
>> de insetos  no segundo mês<br>
><br>
>> parc<-1:108 # Número da parcela no campo<br>
><br>
>> mois2<-cbind(parc,mont2,dead2)<br>
><br>
>> #<br>
><br>
>> nparc=30# Tamanho da amostra independente a ser sorteada com sample ()<br>
><br>
>> Nsim=31# Numero de simulações da função sample<br>
><br>
>><br>
><br>
>> RES=NULL<br>
><br>
>><br>
><br>
>>  for(i in 1:Nsim){<br>
><br>
> +<br>
><br>
> + res=NULL<br>
><br>
> + spl=sample(parc,nparc, replace = TRUE)<br>
><br>
> + resspl=mois2[spl[1:30],]<br>
><br>
> + my.mont=mean(resspl[,2][i])<br>
><br>
> + my.dead=mean(resspl[,3][i])<br>
><br>
> + result=rbind(res,c(my.mont,my.dead,i))<br>
><br>
> + #print(result)<br>
><br>
> + RES=rbind(RES,result)<br>
><br>
> + colnames(RES)=c("monticule","mortalite","simulation")<br>
><br>
> +   }<br>
><br>
>> print(RES)<br>
><br>
>       monticule mortalite simulation<br>
><br>
>  [1,]         5  4.000000          1<br>
><br>
>  [2,]         0  0.000000          2<br>
><br>
>  [3,]         1 12.000000          3<br>
><br>
>  [4,]         1  5.263158          4<br>
><br>
>  [5,]         0  0.000000          5<br>
><br>
>  [6,]         0  0.000000          6<br>
><br>
>  [7,]         2  9.090909          7<br>
><br>
>  [8,]         0  0.000000          8<br>
><br>
>  [9,]         1  0.000000          9<br>
><br>
> [10,]         0  0.000000         10<br>
><br>
> [11,]         1 11.764706         11<br>
><br>
> [12,]         3  6.451613         12<br>
><br>
> [13,]         3 12.000000         13<br>
><br>
> [14,]         3  4.545455         14<br>
><br>
> [15,]         1  3.225806         15<br>
><br>
> [16,]         0  0.000000         16<br>
><br>
> [17,]         3  4.545455         17<br>
><br>
> [18,]         1  0.000000         18<br>
><br>
> [19,]         5 13.636364         19<br>
><br>
> [20,]         2 12.121212         20<br>
><br>
> [21,]         2  6.896552         21<br>
><br>
> [22,]         1  5.263158         22<br>
><br>
> [23,]         2  0.000000         23<br>
><br>
> [24,]        10  8.333333         24<br>
><br>
> [25,]         5 14.285714         25<br>
><br>
> [26,]         8 10.000000         26<br>
><br>
> [27,]         8 10.000000         27<br>
><br>
> [28,]         3  4.761905         28<br>
><br>
> [29,]         0  0.000000         29<br>
><br>
> [30,]         1  3.571429         30<br>
><br>
> [31,]        NA        NA         31<br>
><br>
>><br>
><br>
> Obrigado,<br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
> Alexandre dos Santos<br>
><br>
> Ingenieur forestier, Msc.<br>
><br>
> INRA- Biostatistique et Processus Spatiaux (BioSP)<br>
><br>
> Domaine Saint-Paul<br>
> Site Agroparc<br>
> 84914 -  Avignon - France<br>
> Tél. : +33 (0)6 87 95 16 29<br>
><br>
><br>
><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>
> _______________________________________________<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>
><br>
<br>
<br>
<br>
</div></div><font color="#888888">--<br>
Successful people ask better questions, and as a result, they get<br>
better answers. (Tony Robbins)<br>
</font><div><div></div><div class="h5">_______________________________________________<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>
</div></div></blockquote></div><br></div></div>