[R-br] Recomendação de adubação

[Hélio Gallo Rocha]

Caros listeiros e interessados em recomendação de adubação.

Estou precisando de ajuda para realizar recomendação de adubação.
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.

http://www.dpv24.iciag.ufu.br/new/dpv24/Apostilas/Cafe%20-%205%20Aproximacao%2004.pdf

Uma parte já está pronta, postada em msg anterior pelo mestre Elias.


1. Classifica-se o teor de fósforo em cinco classes, de acordo com o teror
(%) de argila, segundo este CRM:

## Dados
p <- c( 7, 8 , 7, 7, 10);p
arg <- c(10 , 20, 40, 70, 70); arg


## Classificação da argila
class_arg=c(0,15.01,35.01,60.01,100)
nivel_arg=c("Arenoso","Siltoso","Franco_Argiloso","Argiloso")
clas_arg=cut(arg,class_arg,nivel_arg)

## Visualizando a classificação das amostras de argila
data.frame(arg,clas_arg)

# Classificação do Fósforo de acordo com nível de argila
escada <- function(x, y, x.br, y.br, ...) {
  ix <- findInterval(x, x.br, ...)+1
  lr <- lapply(1:nrow(y.br), function(j) {
    i <- which(ix==j)
    list(i=i, r=findInterval(y[i], y.br[j,], ...)+1)
  })
  i <- sort(unlist(lapply(lr, function(x) x$i)),
            index.return=TRUE, method='quick')$ix
  unlist(lapply(lr, function(x) x$r))[i]
}

## teores de fosforo (mg/dm3) para a cultura do café que serão
classificados
## de acordo com quatro classes de  solo de acordo com teror (%) de argila
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))


## Determinando a classe de ferfilidade(cf.cafe) do solo para a cultura do
café
cf.cafe=escada(arg, p, c(15.001, 35.001, 60.001),cafefosf)

## verificando a classe de fertilidade de acordo com quantidade de fósforo
##(mg/dm3) e teor de argila no solo, onde arg = %arg, p = fósforo, cf.cafe
= classe
## é um exemplo com cinco classes
cbind(arg,p,cf.cafe)



2. Para fazer a recomendação de adubação é realizada de acordo com a
produção esperada(sacas café benificiado ha-1), segundo *produ.cafe*

produ.cafe <- c(0,20,30,40,50,60,100)

3. O resultado é de acordo com *aduba.cafe* ( em kg/ha/ano)


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


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.


Agradeço desde já.

-- 
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140710/4677a90a/attachment.html>