ฉันเชื่อว่าสิ่งที่คุณได้รับในคำถามของคุณเกี่ยวข้องกับการตัดข้อมูลโดยใช้ส่วนประกอบหลัก (PC) จำนวนน้อยลง สำหรับการดำเนินการดังกล่าวฉันคิดว่าฟังก์ชั่นprcompเป็นตัวอย่างที่ชัดเจนว่ามันง่ายต่อการมองเห็นการคูณเมทริกซ์ที่ใช้ในการสร้างใหม่
ก่อนอื่นให้ชุดข้อมูลสังเคราะห์Xtคุณดำเนินการ PCA (โดยทั่วไปคุณจะอยู่กึ่งกลางตัวอย่างเพื่ออธิบายพีซีที่เกี่ยวข้องกับเมทริกซ์ความแปรปรวนร่วม:
#Generate data
m=50
n=100
frac.gaps <- 0.5 # the fraction of data with NaNs
N.S.ratio <- 0.25 # the Noise to Signal ratio for adding noise to data
x <- (seq(m)*2*pi)/m
t <- (seq(n)*2*pi)/n
#True field
Xt <-
outer(sin(x), sin(t)) +
outer(sin(2.1*x), sin(2.1*t)) +
outer(sin(3.1*x), sin(3.1*t)) +
outer(tanh(x), cos(t)) +
outer(tanh(2*x), cos(2.1*t)) +
outer(tanh(4*x), cos(0.1*t)) +
outer(tanh(2.4*x), cos(1.1*t)) +
tanh(outer(x, t, FUN="+")) +
tanh(outer(x, 2*t, FUN="+"))
Xt <- t(Xt)
#PCA
res <- prcomp(Xt, center = TRUE, scale = FALSE)
names(res)
ในผลลัพธ์หรือprcompคุณสามารถดูพีซี ( res$x), ค่าลักษณะเฉพาะ ( res$sdev) ที่ให้ข้อมูลเกี่ยวกับขนาดของพีซีแต่ละเครื่องและโหลด ( res$rotation)
res$sdev
length(res$sdev)
res$rotation
dim(res$rotation)
res$x
dim(res$x)
ด้วยการหาค่าลักษณะเฉพาะคุณจะได้รับความแปรปรวนที่อธิบายโดยพีซีแต่ละเครื่อง:
plot(cumsum(res$sdev^2/sum(res$sdev^2))) #cumulative explained variance
สุดท้ายคุณสามารถสร้างข้อมูลที่ถูกตัดทอนโดยใช้พีซีที่สำคัญ (สำคัญ) เท่านั้น:
pc.use <- 3 # explains 93% of variance
trunc <- res$x[,1:pc.use] %*% t(res$rotation[,1:pc.use])
#and add the center (and re-scale) back to data
if(res$scale != FALSE){
trunc <- scale(trunc, center = FALSE , scale=1/res$scale)
}
if(res$center != FALSE){
trunc <- scale(trunc, center = -1 * res$center, scale=FALSE)
}
dim(trunc); dim(Xt)
คุณจะเห็นว่าผลลัพธ์เป็นเมทริกซ์ข้อมูลที่ราบรื่นกว่าเล็กน้อยโดยมีคุณสมบัติขนาดเล็กที่กรองออกมา:
RAN <- range(cbind(Xt, trunc))
BREAKS <- seq(RAN[1], RAN[2],,100)
COLS <- rainbow(length(BREAKS)-1)
par(mfcol=c(1,2), mar=c(1,1,2,1))
image(Xt, main="Original matrix", xlab="", ylab="", xaxt="n", yaxt="n", breaks=BREAKS, col=COLS)
box()
image(trunc, main="Truncated matrix (3 PCs)", xlab="", ylab="", xaxt="n", yaxt="n", breaks=BREAKS, col=COLS)
box()
และนี่คือวิธีการพื้นฐานที่คุณสามารถทำได้นอกฟังก์ชั่น prcomp:
#alternate approach
Xt.cen <- scale(Xt, center=TRUE, scale=FALSE)
C <- cov(Xt.cen, use="pair")
E <- svd(C)
A <- Xt.cen %*% E$u
#To remove units from principal components (A)
#function for the exponent of a matrix
"%^%" <- function(S, power)
with(eigen(S), vectors %*% (values^power * t(vectors)))
Asc <- A %*% (diag(E$d) %^% -0.5) # scaled principal components
#Relationship between eigenvalues from both approaches
plot(res$sdev^2, E$d) #PCA via a covariance matrix - the eigenvalues now hold variance, not stdev
abline(0,1) # same results
ตอนนี้การตัดสินใจว่าจะเก็บพีซีเครื่องใดเป็นคำถามที่แยกต่างหาก - คำถามหนึ่งที่ฉันสนใจในขณะนี้ หวังว่าจะช่วย