<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 http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<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;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-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-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.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> Bom Dia Pessoal,<o:p></o:p></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:<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>> nde<-read.table("Tout_NDE.txt",header=T)#
Todos os dados<o:p></o:p></p>
<p class=MsoNormal>> #<o:p></o:p></p>
<p class=MsoNormal>>
dead2<-nde$Mortalite[nde$Temps==2]-nde$Mortalite[nde$Temps==1]# Mortalidade
no segundo mês<o:p></o:p></p>
<p class=MsoNormal>>
mont2<-nde$Monticules[nde$Temps==2]-nde$Monticules[nde$Temps==1]# Número de
insetos no segundo mês<o:p></o:p></p>
<p class=MsoNormal>> parc<-1:108 # Número da parcela no campo<o:p></o:p></p>
<p class=MsoNormal><span lang=EN-US>> mois2<-cbind(parc,mont2,dead2)<o:p></o:p></span></p>
<p class=MsoNormal>> #<o:p></o:p></p>
<p class=MsoNormal>> nparc=30# Tamanho da amostra independente a ser
sorteada com sample ()<o:p></o:p></p>
<p class=MsoNormal>> Nsim=31# Numero de simulações da função sample<o:p></o:p></p>
<p class=MsoNormal><span lang=EN-US>> <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>> RES=NULL<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>> <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>> for(i in 1:Nsim){<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ res=NULL<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ spl=sample(parc,nparc, replace = TRUE)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ resspl=mois2[spl[1:30],]<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ my.mont=mean(resspl[,2][i])<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ my.dead=mean(resspl[,3][i])<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ result=rbind(res,c(my.mont,my.dead,i))<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ #print(result)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ RES=rbind(RES,result)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+
colnames(RES)=c("monticule","mortalite","simulation")<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>+ }<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>> print(RES)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> </span>monticule mortalite simulation<o:p></o:p></p>
<p class=MsoNormal> [1,] 5 4.000000 1<o:p></o:p></p>
<p class=MsoNormal> [2,] 0 0.000000 2<o:p></o:p></p>
<p class=MsoNormal> [3,] 1 12.000000 3<o:p></o:p></p>
<p class=MsoNormal> [4,] 1 5.263158 4<o:p></o:p></p>
<p class=MsoNormal> [5,] 0 0.000000 5<o:p></o:p></p>
<p class=MsoNormal> [6,] 0 0.000000 6<o:p></o:p></p>
<p class=MsoNormal> [7,] 2 9.090909 7<o:p></o:p></p>
<p class=MsoNormal> [8,] 0 0.000000 8<o:p></o:p></p>
<p class=MsoNormal> [9,] 1 0.000000 9<o:p></o:p></p>
<p class=MsoNormal>[10,] 0 0.000000 10<o:p></o:p></p>
<p class=MsoNormal>[11,] 1 11.764706 11<o:p></o:p></p>
<p class=MsoNormal>[12,] 3 6.451613 12<o:p></o:p></p>
<p class=MsoNormal>[13,] 3 12.000000 13<o:p></o:p></p>
<p class=MsoNormal>[14,] 3 4.545455 14<o:p></o:p></p>
<p class=MsoNormal>[15,] 1 3.225806 15<o:p></o:p></p>
<p class=MsoNormal>[16,] 0 0.000000 16<o:p></o:p></p>
<p class=MsoNormal>[17,] 3 4.545455 17<o:p></o:p></p>
<p class=MsoNormal>[18,] 1 0.000000 18<o:p></o:p></p>
<p class=MsoNormal>[19,] 5 13.636364 19<o:p></o:p></p>
<p class=MsoNormal>[20,] 2 12.121212 20<o:p></o:p></p>
<p class=MsoNormal>[21,] 2 6.896552 21<o:p></o:p></p>
<p class=MsoNormal>[22,] 1 5.263158 22<o:p></o:p></p>
<p class=MsoNormal>[23,] 2 0.000000 23<o:p></o:p></p>
<p class=MsoNormal>[24,] 10 8.333333 24<o:p></o:p></p>
<p class=MsoNormal>[25,] 5 14.285714 25<o:p></o:p></p>
<p class=MsoNormal>[26,] 8 10.000000 26<o:p></o:p></p>
<p class=MsoNormal>[27,] 8 10.000000 27<o:p></o:p></p>
<p class=MsoNormal>[28,] 3 4.761905 28<o:p></o:p></p>
<p class=MsoNormal>[29,] 0 0.000000 29<o:p></o:p></p>
<p class=MsoNormal>[30,] 1 3.571429 30<o:p></o:p></p>
<p class=MsoNormal>[31,] NA NA 31<o:p></o:p></p>
<p class=MsoNormal>><o:p> </o:p></p>
<p class=MsoNormal>Obrigado,<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><b>Alexandre dos Santos<o:p></o:p></b></p>
<p class=MsoNormal><b>Ingenieur forestier, Msc.<o:p></o:p></b></p>
<p class=MsoNormal><b>INRA- Biostatistique et Processus Spatiaux (BioSP)<o:p></o:p></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<o:p></o:p></span></b></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
</div>
</body>
</html>