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.
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
|
|
|
Fiz o seguinte:
solo=c(1:50) # resultado da
análise do solo
argila=c(60,35,15,0) # teor de argila
p1=c(0,2.7,5.4,8,12) # fósforo com + de 60%
de argila
p2=c(0,4,8,12,18) # fósforo com 35 a
60% de argila
p3=c(0,6.6,12,20,30) # fósforo com 15 a 35%
de argila
p4=c(0,10,20,30,45) # fósforo com < 15%
de argila
res=c("Muito.Baixo","Baixo", "Medio", "Alto",
"Muito.alto")
Assim se o resultado de P é 3.5 :
60% de arg, seria classificado como Baixo
35% de arg. seria classificado como muito baixo
andei dando uma olhada na solução do post
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
R-br@listas.c3sl.ufpr.br
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.