protoclust{protoclust}
Verilerimi sınıflandırmak, sınıflara göre renklendirmek ve her sınıf için% 95 güven aralığı için elipsleri üst üste getirmek için kullanılan her değişken çifti için scater grafikleri oluşturarak (ile üretilen ) bir kümelenmenin sonuçlarını görselleştirmek istiyorum elipses sınıfları her değişken çifti altında örtüşür).
Elipslerin çizimini iki farklı şekilde uyguladım ve ortaya çıkan elipsler farklı! (ilk uygulama için daha büyük elipsler!) Eksenlerin merkezleri ve açıları her ikisinde de benzer gibi göründükleri için sadece boyut olarak farklılık gösterdikleri (bazı farklı ölçekleme?). Sanırım bunlardan birini kullanarak yanlış bir şey yapmalıyım (her ikisiyle de değil!) Ya da argümanlarla.
Biri bana neyi yanlış yaptığımı söyleyebilir mi?
Burada iki uygulamanın kodu; her ikisi de bir veri elipsinin ggplot2 dağılım grafiğine nasıl bindirilebileceği yönündeki cevaplara dayanır.
### 1st implementation
### using ellipse{ellipse}
library(ellipse)
library(ggplot2)
library(RColorBrewer)
colorpal <- brewer.pal(10, "Paired")
x <- data$x
y <- data$y
group <- data$group
df <- data.frame(x=x, y=y, group=factor(group))
df_ell <- data.frame()
for(g in levels(df$group)){df_ell <- rbind(df_ell, cbind(as.data.frame(with(df[df$group==g,], ellipse(cor(x, y),scale=c(sd(x),sd(y)),centre=c(mean(x),mean(y))))),group=g))}
p1 <- ggplot(data=df, aes(x=x, y=y,colour=group)) + geom_point() +
geom_path(data=df_ell, aes(x=x, y=y,colour=group))+scale_colour_manual(values=colorpal)
### 2nd implementation
###using function ellipse_stat()
###code by Josef Fruehwald available in: https://github.com/JoFrhwld/FAAV/blob/master/r/stat-ellipse.R
p2 <-qplot(data=df, x=x,y=y,colour=group)+stat_ellipse(level=0.95)+scale_colour_manual(values=colorpal)
İşte iki grafik birlikte (sol grafik p1
uygulama ( ellipse()
):
Verilere buradan ulaşabilirsiniz: https://www.dropbox.com/sh/xa8xrisa4sfxyj0/l5zaGQmXJt
Warning message: In cov.trob(cbind(data$x, data$y)) : Probable convergence failure
, bu da kodu çalıştırdığınızda mı oluyor?