<div dir="ltr">André,<div>Será bem vindo no trabalho</div><div><br></div><div style>Será uma dissertação e artigo.</div><div style>Estamos em fase de colheita de café, vou pro campo jájá..., por isso vou olhar sua rotina com calma de noite.</div>
<div style>Mas como disse vou olhar com calma, mas não tenho intenção de combinar os sorteios</div><div style>Seria "apenas" fazer o sorteio com cada vez menos pontos até o momento em que a correlação e a regressão não seja satisfatória.</div>
<div style>Para isso pensei em fazer a média da correlação dos sorteios.</div><div style>Acho que não fui claro... de noite sento com mais calma e te escrevo.</div><div style><br></div><div style>Obrigado</div><div style>
<br></div><div style>Hélio</div><div><div><br></div><div style><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 20 de junho de 2013 09:28, andrebvs [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4659719h57@n4.nabble.com" target="_blank">ml-node+s2285057n4659719h57@n4.nabble.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>Olá Hélio, estou organizando sua rotina, e pintou um dúvida, ou seja, você deseja saber a média de todas as combinações dos sorteios? tipo:<br><br>Sorteio1[1],Sorteio1[2]<br>Sorteio1[1],Sorteio1[3]</div>
<div>Sorteio1[1],Sorteio1[4]</div>
<div>.<br>.<br>.<br>SorteioSorteio1[99],SorteioSorteio1[100]</div>
<div> </div>
<div>e da mesma forma para o sergundo sorteio, tipo:</div>
<div> </div>
<div>Sorteio2[1],Sorteio2[2]<br>Sorteio2[1],Sorteio2[3]
<div>Sorteio2[1],Sorteio2[4]</div>
.<br>.<br>.<br>SorteioSorteio2[99],SorteioSorteio2[100] <br><br>É isso? Se for, serão 4950 combinações possiveis para cada sorteio.<br><br>Esse trabalho que estais fazendo é para publicar em artigo?<br><br><em>Att.</em><br>
<em>André BVS</em><br><br></div>
<hr style="border-top:1px solid #ccc">
<div>Em 19/06/2013 14:50, <strong>Hélio Gallo Rocha < <a href="http://user/SendEmail.jtp?type=node&node=4659719&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a> ></strong> escreveu:</div>
<div dir="ltr"><div><div class="h5">
<p>Caro André, obrigado pela função, funcionou perfeitamente.</p>
<p> </p>
<p>Explicar melhor a questão da média, espero que não fique cansativo.</p>
<p> </p>
<ol style="margin-top:0cm" type="1" start="1">
<li>Mudei o nome coleta para talhão, que é o real, assim adaptei o CRM para talhão.</li>
<li>Adicionei no final os comandos que fazem o que quero, mas como verá, seria muito trabalhoso.</li>
<li>Se não houver jeito, aceito outras sugestões.</li>
</ol>
<p> </p>
<p>Objetivo:</p>
<ol style="margin-top:0cm" type="1" start="1">
<li>Estudo de previsão de safra em café, contam-se frutos, altura e diâmetro do pé.</li>
<li>Estas medidas são inseridas em quatro funções resultando nos modelos 1 e 2, só postei duas para simplificar.</li>
<li>Serão avaliados em dois níveis, ponto (ponto composto de 5 plantas) e talhão.</li>
<li>Foram amostrados muitos pontos, (200 pontos)</li>
<li>Havendo boa correlação no nível ponto terei a questão: quantos pontos são necessários para que a previsão do talhão seja acurada?</li>
<li>Pesquisas citam dois pontos por hectare</li>
<li>Como super amostrei minha idéia e ir diminuindo o numero de pontos até que a correlação comece a ficar ruim, daí a idéia do sorteio dos pontos.</li>
</ol>
<p> </p>
<p>Se houver ideia estatística melhor, aceito sugestões...</p>
<p> </p>
<p># 1. Dados</p>
<p>talhao=c(rep(1,25),rep(2,25)) ;talhao</p>
<p> </p>
<p>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,</p>
<p>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,</p>
<p>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,</p>
<p>4.54,4.72,7.69,3.45)</p>
<p> </p>
<p>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,</p>
<p>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,</p>
<p>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,</p>
<p>4.05,5.02,5.65,3.04,3.61,9.48,4.13)</p>
<p> </p>
<p>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,</p>
<p>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,</p>
<p>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,</p>
<p>6.36,6.84,6.48,7.58,7.35,4.99,5.3,11.99,6.21)</p>
<p> </p>
<p>dados=cbind(talhao,colheita,modelo1,modelo2);dados</p>
<p>x=cbind(colheita,modelo1,modelo2)</p>
<p>length(colheita)</p>
<p>cor(x)</p>
<p> </p>
<p> </p>
<p>#2. Separando dados da talhao 1 da 2</p>
<p>tal1=subset(dados,talhao==1) #talhão 1</p>
<p>tal2=subset(dados,talhao==2) #talhão 2</p>
<p> </p>
<p> </p>
<p>#3.1 Sorteando a talhao 1</p>
<p>L1</p>
<p> </p>
</div></div><p><span lang="EN-US">for (i in 1:100)</span></p>
<p><span lang="EN-US">{</span></p>
<p><span lang="EN-US"> stal1 </span></p>
<p><span lang="EN-US"> </span>L1[i]</p><div class="im">
<p>}</p>
<p> </p>
<p>#verificando sorteio do talhão 1</p>
<p>L1</p>
<p>#L1[[1]][,1] #dados do talhão 1, sorteio1, colheita</p>
<p>#L1[[1]][,2] #dados do talhão 1, sorteio1, modelo 1</p>
<p>#L1[[1]][,3] #dados do talhão 1, sorteio1, modelo 2</p>
<p> </p>
<p># 3.1.1 Correlação dos dados do talhão 1, sorteio 1</p>
<p>t1_s1=cor(cbind( L1[[1]][,1],L1[[1]][,2],L1[[1]][,3]))</p>
<p> </p>
<p>#verificando a correlação</p>
<p>t1_s1</p>
<p> </p>
<p># 3.1.2 Correlação dos dados do talhão 1, centésimo sorteio</p>
<p>t1_s100=cor(cbind( (L1[[100]][,1]),(L1[[100]][,2]),(L1[[100]][,3])) )</p>
<p>#verificando correlação do centésimo</p>
<p>t1_s100</p>
<p> </p>
<p># 3.1.3 média da correlação do modelo1, só fiz a média do primeiro sorteio e do centésimo</p>
<p>media_t1_s1_m1=(t1_s1[[2]]+t1_s100[[2]])/2 ;media_t1_s1_m1</p>
<p> </p>
<p># 3.1.4 média da correlação do modelo 2, primeiro e centésimo sorteio</p>
<p>media_t1_s1_m2=(t1_s1[[3]]+t1_s100[[3]])/2 ;media_t1_s1_m2</p>
<p> </p>
<p> </p>
<p># 3.2 Sorteio da talhao 2</p>
<p><span lang="EN-US">L2 </span></p>
<p><span lang="EN-US"> </span></p>
</div><p><span lang="EN-US">for (i in 1:100)</span></p>
<p><span lang="EN-US">{</span></p>
<p><span lang="EN-US"> stal2 </span></p>
<p><span lang="EN-US"> </span>L2[i]</p><div class="im">
<p>}</p>
<p> </p>
<p>L2</p>
<p> </p>
<p># 3.2.1 Correlação dos dados do talhão 2, primeiro sorteio</p>
<p>t2_s1=cor(cbind( L2[[1]][,1],L2[[1]][,2],L2[[1]][,3]))</p>
<p>#verificando a correlação</p>
<p>t2_s1</p>
<p> </p>
<p># 3.2.2 Correlação dos dados do talhão 1, centésimo sorteio</p>
<p>t2_s100=cor(cbind( (L2[[100]][,1]),(L2[[100]][,2]),(L2[[100]][,3])) )</p>
<p>#verificando correlação do centésimo</p>
<p>t2_s100</p>
<p> </p>
<p># 3.2.3 média da correlação do modelo1, só fiz a média do primeiro sorteio e do centésimo</p>
<p>media_t2_s1_m1=(t2_s1[[2]]+t2_s100[[2]])/2 ;media_t2_s1_m1</p>
<p># 3.2.4 média da correlação do modelo 2, primeiro e centésimo sorteio</p>
<p>media_t2_s1_m2=(t2_s1[[3]]+t2_s100[[3]])/2 ;media_t2_s1_m2</p>
<p> </p>
<p>###fim</p>
</div></div>
<div class="gmail_extra"><br><br>
<div class="gmail_quote">Em 18 de junho de 2013 21:59, andrebvs [via R-br] <span dir="ltr"><<a href="http://../../../undefined/compose?to=ml-node" rel="nofollow" link="external" target="_blank">ml-node+s2285057n4659691h89@...</a>></span> escreveu:<br>
<blockquote style="border-left:2px solid #cccccc;padding:0 1em" class="gmail_quote">
<div><div class="im">Olá Hélio, o sorteio das 100 vezes com 20 amostras:<br><br>L <br></div><div class="im">for (i in 1:100)<br>{<br></div> r L[i] }</div><div><div class="h5">
<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 < <a href="http://user/SendEmail.jtp?type=node&node=4659691&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a> ></strong> escreveu:</div>
<div>
<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>
</div>
</div>
_______________________________________________ <br>R-br mailing list <br><a href="http://user/SendEmail.jtp?type=node&node=4659691&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a> <br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível. <br> <br><hr noshade size="1">
<div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
<a href="http://r-br.2285057.n4.nabble.com/R-br-sortear-100x-e-correlacionar-tp4659687p4659691.html" rel="nofollow" link="external" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-sortear-100x-e-correlacionar-tp4659687p4659691.html</a></div>
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">To unsubscribe from R-br, <a rel="nofollow" link="external">click here</a>.<br> <a style="font:9px serif" href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" link="external" target="_blank">NAML</a></div>
</div></div></blockquote>
</div><div><div class="h5">
<br><br>
<div> </div>
-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho</div></div></div><div><div class="h5">
<div> </div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4659719&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr noshade size="1" color="#cccccc">
</div></div><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div><div class="h5">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div></div><a href="http://r-br.2285057.n4.nabble.com/R-br-sortear-100x-e-correlacionar-tp4659687p4659719.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-sortear-100x-e-correlacionar-tp4659687p4659719.html</a>
</div><div class="HOEnZb"><div class="h5">
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>