Porque em vez de usar esse tipo escada de classificação você não usa um tipo suave? Para mim, este último faz mais sentido. Veja exemplo de como obter isso com a sua 'escada':
### Tabela da quinta aproximação
### Nível de Fósforo mg/dm3
###Argila% muito baixo Baixo Médio Bom Muito Bom
###60-100 <2.7(menor e igual) 2.8-5.4 5.5-8 8.1-12 >12
###35-60 <4 4.1-8 8.1-12 12.1-18 >18
###15-35 <6.6 6.7-12 12.1-20 20.1-30 >30
### 0-15 <10 10.1-20 20.1-30 30.1-45 >45
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]
}
### exemplo simulado
n <- 30000
x <- 90*runif(n)
y <- runif(n, 0, 100-x)*.7
cfosf <- rbind(c(10, 20, 30, 45), c(6.6, 12, 20, 30),
c(4, 8, 12, 18), c(2.7, 5.4, 8, 12))
cfosf
res <- escada(x, y, c(15, 35, 60), cfosf)
table(res)
plot(x, y, col=res, pch=19, las=1, asp=1, cex=0.3,
xlab="Argila", ylab="Fósforo")
abline(v=c(15, 35, 60), col=gray(.7))
for (j in 1:nrow(cfosf))
segments(c(0, 15, 35, 60)[j], cfosf[j,],
c(15, 35, 60, 100)[j], cfosf[j,], lty=3)
legend("topright", c("Muito baixo", "Baixo", "Médio", "Muito bom"),
pch=19, col=1:5, bty='n', title='Textura')
segments(0,100,100,0, lty=3)
### usando esses dados simulados
### (ideal e' usar dados reais)
### para obter uma funcao suave
require(MASS)
aj <- polr(resp ~ x*y, data.frame(resp=factor(res, ordered=TRUE)))
prd <- predict(aj)
plot(x, y, col=unclass(prd), pch=19, las=1, asp=1, cex=0.3,
xlab="Argila", ylab="Fósforo")
legend("topright", c("Muito baixo", "Baixo", "Médio", "Muito bom"),
pch=19, col=1:5, bty='n', title='Textura')
### classificando um novo solo (Argila=36, fosforo=19)
### usando a escada:
escada(36, 19, c(15, 35, 60), cfosf)
### suave:
predict(aj, newdata=data.frame(x=36, y=19))
Att.
Elias.
Em 18/05/2013 09:28, Hélio Gallo Rocha < [hidden email] > escreveu:Bom dia a todos,Para recomendar adubação fosfatada é necessário classificar a quantidade de fósforo da análise em conjunto a textura do solo.
Tabela da quinta aproximação
Nível de Fósforo mg/dm3
Argila%
muito baixo
Baixo
Médio
Bom
Muito Bom
60-100
<2.7 (menor e igual)
2.8-5.4
5.5-8
8.1-12
>12
35-60
<4
4.1-8
8.1-12
12.1-18
>18
15-35
<6.6
6.7-12
12.1-20
20.1-30
>30
0-15
<10
10.1-20
20.1-30
30.1-45
>45
Uso do ifelse
A saida desta função seria casada com uma recomendação,Grato a todos--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list [hidden email] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
_______________________________________________
R-br mailing list
[hidden email]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
http://r-br.2285057.n4.nabble.com/R-br-Classificar-fosforo-e-textura-do-solo-tp4659355p4659361.htmlIf you reply to this email, your message will be added to the discussion below: