[R-br] Problema com sample()

[Alexandre dos Santos]

Obrigado Benilton,

 

     Agora já sei onde esta o problema, mais é que eu preciso do valor das
duas variáveis  na segunda e terceira colunas associadas aos valores
sorteados para a variável parc, tentei fazer resspl=mois2[spl,], também sem
sucesso, pois meu objetivo é fazer um sorteio de n=30 e repetir 99 vezes, 

Alexandre 

 

 

From: r-br-bounces em listas.c3sl.ufpr.br
[mailto:r-br-bounces em listas.c3sl.ufpr.br] On Behalf Of Benilton Carvalho
Sent: segunda-feira, 25 de abril de 2011 08:15
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] Problema com sample()

 

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).

On 25 Apr 2011 10:18, "Alexandre dos Santos"
<alexandresantosbr em yahoo.com.br> wrote:

    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

 


_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110425/252b4594/attachment.html>