Bom Dia Pessoal, 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:
nde<-read.table("Tout_NDE.txt",header=T)# Todos os dados
#
dead2<-nde$Mortalite[nde$Temps==2]-nde$Mortalite[nde$Temps==1]# Mortalidade no segundo mês
mont2<-nde$Monticules[nde$Temps==2]-nde$Monticules[nde$Temps==1]# Número de insetos no segundo mês
parc<-1:108 # Número da parcela no campo
mois2<-cbind(parc,mont2,dead2)
#
nparc=30# Tamanho da amostra independente a ser sorteada com sample ()
Nsim=31# Numero de simulações da função sample
RES=NULL
for(i in 1:Nsim){
+ + res=NULL + spl=sample(parc,nparc, replace = TRUE) + resspl=mois2[spl[1:30],] + my.mont=mean(resspl[,2][i]) + my.dead=mean(resspl[,3][i]) + result=rbind(res,c(my.mont,my.dead,i)) + #print(result) + RES=rbind(RES,result) + colnames(RES)=c("monticule","mortalite","simulation") + }
print(RES)
monticule mortalite simulation [1,] 5 4.000000 1 [2,] 0 0.000000 2 [3,] 1 12.000000 3 [4,] 1 5.263158 4 [5,] 0 0.000000 5 [6,] 0 0.000000 6 [7,] 2 9.090909 7 [8,] 0 0.000000 8 [9,] 1 0.000000 9 [10,] 0 0.000000 10 [11,] 1 11.764706 11 [12,] 3 6.451613 12 [13,] 3 12.000000 13 [14,] 3 4.545455 14 [15,] 1 3.225806 15 [16,] 0 0.000000 16 [17,] 3 4.545455 17 [18,] 1 0.000000 18 [19,] 5 13.636364 19 [20,] 2 12.121212 20 [21,] 2 6.896552 21 [22,] 1 5.263158 22 [23,] 2 0.000000 23 [24,] 10 8.333333 24 [25,] 5 14.285714 25 [26,] 8 10.000000 26 [27,] 8 10.000000 27 [28,] 3 4.761905 28 [29,] 0 0.000000 29 [30,] 1 3.571429 30 [31,] NA NA 31
Obrigado, Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29