KS 代码翻译(Kernel Density Estimation and KernelDiscriminant Analysis for Multivariate Data inR)

library("ks")
samp <- 200
mus <- rbind(c(-2,2), c(0,0), c(2,-2))
Sigmas <- rbind(diag(2), matrix(c(0.8, -0.72, -0.72, 0.8), nrow = 2),diag(2))
cwt <- 3/11
props <-c((1-cwt)/2, cwt, (1-cwt)/2)
x<-rmvnorm.mixt(n = samp, mu = mus, Sigma = Sigmas, props = props)
Hpi1 <- Hpi(x = x, pilot = "amse", pre = "scale")
Hpi2 <- Hpi(x = x, pilot = "samse", pre = "scale")
Hpi3 <- Hpi(x = x, pilot = "amse", pre = "sphere")
Hpi4 <- Hpi(x = x, pilot = "samse", pre = "sphere")
Hpi5 <- Hpi.diag(x = x, pilot = "amse", pre = "scale")
Hpi6 <- Hpi.diag(x = x, pilot = "samse", pre = "scale")

Hpi1_1<-kde(x = x, H = Hpi1)
Hpi2_1<-kde(x = x, H = Hpi2)
Hpi3_1<-kde(x = x, H = Hpi3)
Hpi4_1<-kde(x = x, H = Hpi4)
Hpi5_1<-kde(x = x, H = Hpi5)
Hpi6_1<-kde(x = x, H = Hpi6)

> # plot(Hpi1_1,xlab="x",ylab="y",main="amse scale")

> # plot(Hpi2_1,xlab="x",ylab="y",main="samse scale")

plot(Hpi3_1,xlab="x",ylab="y",main="amse sphere")

plot(Hpi4_1,xlab="x",ylab="y",main="samse sphere")

plot(Hpi5_1,xlab="x",ylab="y",main="Hpi.diag amse scale")

plot(Hpi6_1,xlab="x",ylab="y",main="Hpi.diag samse scale")

Hlscv1 <- Hlscv(x = x)
Hlscv2 <- Hlscv.diag(x = x)
Hbcv1 <- Hbcv(x = x, whichbcv = 1)
Hbcv2 <- Hbcv.diag(x = x, whichbcv = 1)
Hscv1 <- Hscv(x = x, pre = "scale")
Hscv2 <- Hscv.diag(x = x, pre = "scale")

kde_Hlscv1<-kde(x = x, H = Hlscv1)
plot(kde_Hlscv1)

kde_Hlscv2<-kde(x = x, H = Hlscv2)
plot(kde_Hlscv2)

> # kde_Hbcv1<-kde(x = x, H = Hbcv1)
> # plot(kde_Hbcv1)

> # kde_Hbcv2<-kde(x = x, H = Hbcv2)
> # plot(kde_Hbcv2)

> # kde_Hscv1<-kde(x = x, H = Hscv1)
> # plot(kde_Hscv1)

> # kde_Hscv2<-kde(x = x, H = Hscv2)
> # plot(kde_Hscv2)

x <- rmvnorm.mixt(10000, mus, Sigmas, props)
Hpi(x = x, binned = TRUE, pilot = "samse")
Hdiag <- Hscv.diag(x = x, binned = TRUE)
kde1<-kde(x = x, H = Hdiag, binned = TRUE)
plot(kde1)

library("MASS")
library("rgl")
library("misc3d")
data("iris")
ir<-iris[,1:3]
ir.group <- iris[,5]
Hpi1 <- Hkda(x = ir, x.group = ir.group, bw = "plugin",pilot = "samse", pre = "sphere")
Hpi2 <- Hkda.diag(x = ir, x.group = ir.group, bw = "plugin",pilot = "samse", pre = "scale")
Hscv1 <- Hkda(x = ir, x.group = ir.group, bw = "scv", pre = "sphere")
Hscv2 <- Hkda.diag(x = ir, x.group = ir.group, bw = "scv", pre = "scale")
kda<-kda(x = ir, x.group = ir.group, Hs = Hpi1)#HS可以换成Hpi1、Hpi2、Hscv1等
plot(kda)

mycolors <- c(‘royalblue1’, ‘darkcyan’, ‘oldlace’)
irKaTeX parse error: Expected 'EOF', got '#' at position 43: …ic(ir.group) ] #̲绘制三维散点图 plot3d(…Sepal.Length, y=ir$‘Sepal.Width‘,z=ir$Petal.Length,
col = ir$color,
type = ‘s’,
radius = .1,
xlab=“Sepal Length”, ylab=“Sepal Width”, zlab=“Petal Length”)

compare.kda.cv(x = ir, x.group = ir.group, bw = “plugin”,pilot = “samse”, pre = “sphere”)
compare.kda.diag.cv(x = ir, x.group = ir.group, bw = “plugin”,pilot = “samse”, pre = “scale”)
compare.kda.cv(x = ir, x.group = ir.group, bw = “scv”, pre = “sphere”)
compare.kda.diag.cv(x = ir, x.group = ir.group, bw = “scv”,pre = “scale”)