<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=PT-BR link=blue vlink=purple>
<div class=Section1>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> Obrigado Rodrigo e Beniltol,<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'> Funcionou era isso mesmo que estava fazendo errado.<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Abraço,<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Alexandre dos Santos<o:p></o:p></span></b></p>
<p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Ingenieur forestier, Msc.<o:p></o:p></span></b></p>
<p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>INRA- Biostatistique et Processus Spatiaux (BioSP)<o:p></o:p></span></b></p>
<p class=MsoNormal><b><span lang=EN-US style='font-size:11.0pt;font-family:
"Calibri","sans-serif";color:#1F497D'>Domaine Saint-Paul<br>
Site Agroparc <br>
84914 - Avignon - France<br>
Tél. : +33 (0)6 87 95 16 29<o:p></o:p></span></b></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'>
<p class=MsoNormal><b><span lang=EN-US style='font-size:10.0pt;font-family:
"Tahoma","sans-serif"'>From:</span></b><span lang=EN-US style='font-size:10.0pt;
font-family:"Tahoma","sans-serif"'> r-br-bounces@listas.c3sl.ufpr.br
[mailto:r-br-bounces@listas.c3sl.ufpr.br] <b>On Behalf Of </b>Rodrigo Coster<br>
<b>Sent:</b> segunda-feira, 25 de abril de 2011 10:26<br>
<b>To:</b> r-br@listas.c3sl.ufpr.br<br>
<b>Subject:</b> Re: [R-br] Problema com sample()<o:p></o:p></span></p>
</div>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>É 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.<o:p></o:p></p>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal>Creio que isso aqui faz o que tu quer (em negrito o que eu
mudei). <o:p></o:p></p>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<div>
<p class=MsoNormal>for(i in 1:Nsim){<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal> res=NULL<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal> spl=sample(parc,nparc, replace = TRUE)<o:p></o:p></p>
</div>
<p class=MsoNormal> resspl=mois2[spl[1:30],]<o:p></o:p></p>
<div>
<p class=MsoNormal><b> my.mont=mean(resspl[,2])</b><o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><b> my.dead=mean(resspl[,3])</b><o:p></o:p></p>
</div>
<div>
<p class=MsoNormal> result=rbind(res,c(my.mont,my.dead,i))<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal> #print(result)<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal> RES=rbind(RES,result)<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal> colnames(RES)=c("monticule","mortalite","simulation")<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>}<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<p class=MsoNormal><o:p> </o:p></p>
<div>
<p class=MsoNormal>2011/4/25 Benilton Carvalho <<a
href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>><o:p></o:p></p>
<p class=MsoNormal>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>>:<o:p></o:p></p>
<div>
<div>
<p class=MsoNormal style='margin-bottom:12.0pt'>> 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>
<o:p></o:p></p>
</div>
</div>
<p class=MsoNormal><span style='color:#888888'>--<br>
Successful people ask better questions, and as a result, they get<br>
better answers. (Tony Robbins)</span><o:p></o:p></p>
<div>
<div>
<p class=MsoNormal>_______________________________________________<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><o:p></o:p></p>
</div>
</div>
</div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
</div>
</div>
</body>
</html>