[R-br] Dúvida lógica do código R.

Diogo Jerônimo diogojose21 em yahoo.com.br
Seg Fev 14 13:49:57 -03 2022


Boa tarde a todos!!! Dúvida boba, mas tenho de explicar antes: estou fazendo uma análise de cluster, para poder identificar resultados iguais de laboratórios analíticos. O resultado foi um banco parecido com esse reproduzível:
library(dplyr)
x01<-c(8.49,8.62,8.63,8.49,8.44,2)
x02<-c(9.18,9.21,9.15,8.86,9.36,3)x03<-c(8.49,8.62,8.63,8.49,8.44,2)x04<-c(9.74,9.69,9.81,9.25,9.98,3)x05<-c(9.74,9.69,9.81,9.25,9.98,3)x06<-c(8.15,9.52,9.89,9.95,8.33,2)x07<-c(8.32,8.73,9.60,9.04,9.97,3)x08<-c(7.42,7.01,7.31,7.20,7.61,1)x09<-c(9.79,8.76,9.68,8.31,9.37,3)x10<-c(7.42,7.01,7.31,7.20,7.61,1)
dados<-data.frame(rbind(x01,x02,x03,x04,x05,x06,x07,x08,x09,x10))
names(dados)<-c("v1","v2","v3","v4","v5","Cluster")
Para facilitar essa identificação, ordenei os resultados por cluster, pelos dados das variáveis, e então calculei a média em cada linha:
dados<-dados %>% arrange(Cluster, v1,v2,v3,v4,v5) 
for (i in (1:nrow(dados))){dados$Média[i]<-mean(dados$v1[i],dados$v2[i],dados$v3[i],dados$v4[i],dados$v5[i],na.rm = TRUE)}
MEU PROBLEMA: para finalmente identificar os resultados iguais, eu apliquei esse código abaixo, onde em resumo, se o cluster for igual e a média do resultado do laboratório da linha i for igual a média do resultado do laboratório da linha i+1, então ambos os laboratórios responderam resultados iguais. Apliquei o código abaixo, e ele está resultando em erro: 
i=1for (i in (1:nrow(dados))){if ((dados$Cluster[i]==dados$Cluster[i+1]) & (dados$Média[i]==dados$Média[i+1])){dados$sit[i] == "Resultado Igual"dados$sit[i+1] == "Resultado Igual"}else {}i=i+1}
Alguém poderia dizer qual o meu erro? Ou se faltou algo no código?
Muito agradecido pela atenção, e desculpa a extensão!!!
Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQIhttp://lattes.cnpq.br/8996149312896520


-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20220214/1593455d/attachment.htm>


Mais detalhes sobre a lista de discussão R-br