<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:12.727272033691406px">Caros listeiros e interessados em recomendação de adubação.</span><div style="font-family:arial,sans-serif;font-size:12.727272033691406px">
<br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px">Estou precisando de ajuda para realizar recomendação de adubação.</div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px">
Como exemplo utilizei a cultura do café com cálculos realizados de acordo com  a quinta aproximação para Minas Gerais,  que pode ser obtido no seguinte link, que no caso, é especifico para a cultura do café, e a tabela em questão se encontra na pag. 4.</div>
<div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><a href="http://www.dpv24.iciag.ufu.br/new/dpv24/Apostilas/Cafe%20-%205%20Aproximacao%2004.pdf" target="_blank">http://www.dpv24.iciag.ufu.br/new/dpv24/Apostilas/Cafe%20-%205%20Aproximacao%2004.pdf</a><br>
</div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px">Uma parte já está pronta, postada em msg anterior pelo mestre Elias.</div>
<div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px">
1. Classifica-se o teor de fósforo em cinco classes, de acordo com o teror (%) de argila, segundo este CRM:</div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px">
<div>## Dados</div><div>p <- c( 7, 8 , 7, 7, 10);p</div><div>arg <- c(10 , 20, 40, 70, 70); arg</div><div><br></div><div><br></div><div>## Classificação da argila</div><div>class_arg=c(0,15.01,35.01,60.01,100)</div>
<div>nivel_arg=c("Arenoso","Siltoso","Franco_Argiloso","Argiloso")</div><div>clas_arg=cut(arg,class_arg,nivel_arg)</div><div><br></div><div>## Visualizando a classificação das amostras de argila </div>
<div>data.frame(arg,clas_arg) </div><div><br></div><div># Classificação do Fósforo de acordo com nível de argila</div><div>escada <- function(x, y, <a href="http://x.br/" target="_blank">x.br</a>, <a href="http://y.br/" target="_blank">y.br</a>, ...) {</div>
<div>  ix <- findInterval(x, <a href="http://x.br/" target="_blank">x.br</a>, ...)+1</div><div>  lr <- lapply(1:nrow(<a href="http://y.br/" target="_blank">y.br</a>), function(j) {</div><div>    i <- which(ix==j)</div>
<div>    list(i=i, r=findInterval(y[i], <a href="http://y.br/" target="_blank">y.br</a>[j,], ...)+1)</div><div>  })</div><div>  i <- sort(unlist(lapply(lr, function(x) x$i)),</div><div>            index.return=TRUE, method='quick')$ix</div>
<div>  unlist(lapply(lr, function(x) x$r))[i]</div><div>}</div><div><br></div><div>## teores de fosforo (mg/dm3) para a cultura do café que serão classificados </div><div>## de acordo com quatro classes de  solo de acordo com teror (%) de argila</div>
<div>cafefosf <- rbind(c(7.5,15.1,22.6,33.81),c(5.01,9.1,15.1,22.51), c(3.1,6.1,9.1,13.51 ), c(2,4.1,6.1,9.01))</div><div><br></div><div><br></div><div>## Determinando a classe de ferfilidade(cf.cafe) do solo para a cultura do café</div>
<div>cf.cafe=escada(arg, p, c(15.001, 35.001, 60.001),cafefosf)</div><div><br></div><div>## verificando a classe de fertilidade de acordo com quantidade de fósforo</div><div>##(mg/dm3) e teor de argila no solo, onde arg = %arg, p = fósforo, cf.cafe = classe </div>
<div>## é um exemplo com cinco classes<br></div><div>cbind(arg,p,cf.cafe)</div></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px">
<br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px">2. Para fazer a recomendação de adubação é realizada de acordo com a produção esperada(sacas café benificiado ha-1), segundo <b>produ.cafe</b></div>
<div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="font-family:arial,sans-serif;font-size:12.727272033691406px"><div>produ.cafe <- c(0,20,30,40,50,60,100)          </div><div>
<br></div><div>3. O resultado é de acordo com <b>aduba.cafe</b> ( em kg/ha/ano)</div><div><br></div><div><br></div><div>aduba.cafe <- rbind(c(80,60,40,20,0),c(70,55,35,18,0),c(60,50,30,15,0),c(50,40,25,0,0),c(30,20,10,0,0));aduba.cafe</div>
<div><br></div><div><br></div><div>Para verificar se deu certo, a primeira linha do resultado de cbind(cf.cafe) que é 1, com uma produção estimada em 25 sacas daria como resultado a recomendação de 30 kg/ha/ano de fósforo.</div>
<div><br></div><div><br></div><div>Agradeço desde já.</div><div class=""><div id=":17l" class="" tabindex="0"><img class="" src="https://ssl.gstatic.com/ui/v1/icons/mail/images/cleardot.gif"></div></div></div><div><br></div>
-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>