<div>Olá Hélio, o sorteio das 100 vezes com 20 amostras:<br /><br />L <- list()<br /><br />for (i in 1:100)<br />{<br />  r <- do.call(rbind, by(dados, dados[,1], function(x, ...) x[sample(1:nrow(x), ...), ], size=10))<br />  L[i] <- list(cbind(r$colheita,r$modelo1,r$modelo2))<br />} </div>
<div> </div>
<div>Agora, a segunda parte da sua solicitação não entendi muito bem, seje mais claro.<br /><br /><em>Att.</em><br /><em>André BVS</em><br /><br /></div>
<hr style="border-top: 1px solid #ccc;" />
<div>Em 18/06/2013 20:51, <strong>Hélio Gallo Rocha < heliogallorocha@gmail.com ></strong> escreveu:</div>
<div dir="ltr">Bom dia a todos,
<div> </div>
<div>Peço ajuda na seguinte questão destes dados que comparam dois modelos de previsão de safra em café, onde :</div>
<div>colheita = dados da parcela colhida</div>
<div>modelo 1 de previsão</div>
<div>modelo 2 de previsão</div>
<div> </div>
<div>CRM</div>
<div> </div>
<div>
<div>
<div>
<div>coleta=c(rep(1,25),rep(2,25))</div>
<div> </div>
<div>colheita=c(10.2,9.22,8.6,15.57,11.93,12.28,12.43,16.47,18.72,12.77,16.03,13.4,16.95,</div>
<div>15.4,18.64,14.38,15.02,16.98,14.38,16.35,18.28,16.06,17.92,12.85,13.83,4.48,3.35,7.24,</div>
<div>7.73,8.45,4.83,3.8,2.2,5.99,4.24,4.67,8.31,4.25,7.32,5.45,4.24,5.72,5.09,4.84,5.62,5.5,</div>
<div>4.54,4.72,7.69,3.45)</div>
<div> </div>
<div>modelo1=c(17.08,12.24,13.44,18.52,11.82,11.37,16.17,16.92,16.66,13.32,18.95,14.78,</div>
<div>15.54,16.12,25.62,18.2,15.86,20.5,20.47,20.45,17.15,22.15,20.29,13.82,12.86,3.13,</div>
<div>1.84,6.9,6.35,8.9,5.64,4.79,2,4.64,5.06,4.85,5.61,4.82,6.59,4.15,4.46,5.1,4.78,</div>
<div>4.05,5.02,5.65,3.04,3.61,9.48,4.13)</div>
<div> </div>
<div>modelo2=c(12.43,10.99,10.32,12.46,9.85,10.15,12.68,12.46,12.16,9.92,12.58,10.63,</div>
<div>11.18,10.6,15.71,13.22,10.82,12.85,13.48,14.13,11.41,14.09,12.23,10.28,8.75,</div>
<div>5.2,3.56,9.89,9.03,12,8.16,7.15,3.2,6.58,6.84,7.37,8.5,7.87,8.38,5.58,6.62,</div>
<div>6.36,6.84,6.48,7.58,7.35,4.99,5.3,11.99,6.21)</div>
<div> </div>
<div>dados=cbind(coleta,colheita,modelo1,modelo2);dados</div>
<div>x=cbind(colheita,modelo1,modelo2)</div>
<div>cor(x)</div>
<div> </div>
<div> </div>
<div>#sorteio de 20 medidas,</div>
<div>r=do.call(rbind, by(dados, dados[,1], function(x, ...) x[sample(1:nrow(x), ...), ], size=10))</div>
<div>y=cbind(r$colheita,r$modelo1,r$modelo2);y</div>
<div>cor(y)</div>
</div>
</div>
<div> </div>
<div> </div>
<div>Gostaria de sortear 100 vezes com 20 amostras, sendo o resultado a média da correlação e a média dos dados</div>
<div> </div>
<div>Grato a todos</div>
<div> </div>
<div> </div>
<div> </div>
-- <br />Hélio Gallo Rocha<br />IFSULDEMINAS - Câmpus Muzambinho</div>
</div>
<div> </div>