পিসিএ বিপ্লট চারটি অক্ষ কি?
উত্তর:
আপনি কি বোঝাতে চান, উদাহরণস্বরূপ, নিম্নলিখিত কমান্ডটি যে প্লটটিতে ফিরে এসেছে?
biplot(prcomp(USArrests, scale = TRUE))
যদি হ্যাঁ, তবে উপরে এবং ডান অক্ষগুলি প্লটটিতে লাল তীরগুলি (ভেরিয়েবলগুলি চিত্রিত করার পয়েন্টগুলি) ব্যাখ্যার জন্য ব্যবহার করা উচিত।
যদি আপনি জানেন যে কীভাবে মূল উপাদান বিশ্লেষণ কাজ করে এবং আপনি আর কোডটি পড়তে পারেন তবে নীচের কোডটি আপনাকে দেখায় যে চূড়ান্ত চক্রান্ত করার আগে prcomp()প্রাথমিকভাবে ফলাফলগুলি কীভাবে আচরণ করা biplot.prcomp()হয় biplot.default()। আপনি যখন প্লট করবেন তখন এই দুটি ফাংশনটিকে পটভূমিতে ডাকা হয় biplot()এবং নিম্নলিখিত সংশোধিত কোড উদ্ধৃতিটি আসে biplot.prcomp()।
x<-prcomp(USArrests, scale=TRUE)
choices = 1L:2L
scale = 1
pc.biplot = FALSE
scores<-x$x
lam <- x$sdev[choices]
n <- NROW(scores)
lam <- lam * sqrt(n)
lam <- lam^scale
yy<-t(t(x$rotation[, choices]) * lam)
xx<-t(t(scores[, choices])/lam)
biplot(xx,yy)শীঘ্রই, উপরের উদাহরণে, ভেরিয়েবল লোডিংয়ের ম্যাট্রিক্স ( x$rotation) x$sdevপর্যবেক্ষণের সংখ্যার বর্গমূলের মূল উপাদানগুলির ( ) গুণমানের বিচ্যুতি দ্বারা মাপানো হয় । এটি প্লটটিতে যা দেখা যায় তার উপরে এবং ডান অক্ষগুলির জন্য স্কেল সেট করে।
পরিবর্তনশীল লোডিংগুলি স্কেল করার জন্য অন্যান্য পদ্ধতি রয়েছে are এগুলি দেওয়া হয় যেমন আর প্যাকেজ ভেগান দ্বারা।
আমার কাছে বিপ্লটের জন্য আরও ভাল ভিজ্যুয়ালাইজেশন রয়েছে। নিম্নলিখিত চিত্র পরীক্ষা করুন।
পরীক্ষায়, আমি 3 ডি পয়েন্টগুলিকে 2 ডি (সিমুলেটেড ডেটা সেট) ম্যাপ করার চেষ্টা করছি।
2 ডি-তে বিপ্লট বোঝার কৌশলটি 3 ডি তে একই জিনিসটি দেখার জন্য সঠিক কোণটি খুঁজেছে। সমস্ত ডেটা পয়েন্ট নম্বরযুক্ত, আপনি ম্যাপিং স্পষ্ট দেখতে পাবেন।
ফলাফলগুলি পুনরুত্পাদন করার কোড এখানে।
require(rgl)
set.seed(0)
feature1=round(rnorm(50)*10+20)
feature2=round(rnorm(50)*10+30)
feature3=round(runif(50)*feature1)
d=data.frame(feature1,feature2,feature3)
head(d)
plot(feature1,feature2)
plot(feature2,feature3)
plot(feature1,feature3)
plot3d(d$feature1, d$feature2, d$feature3, type = 'n')
points3d(d$feature1, d$feature2, d$feature3, color = 'red', size = 10)
shift <- matrix(c(-2, 2, 0), 12, 3, byrow = TRUE)
text3d(d+shift,texts=1:50)
grid3d(c("x", "y", "z"))
pr.out=prcomp(d,scale.=T)
biplot(pr.out)
grid()biplotআর- এর কমান্ড দ্বারা উত্পাদিত বাইপ্লট-এ এটি নয় এবং ডানদিকে আপনার চিত্রে পুনরুত্পাদন করা হয়েছে: সেখানে বিন্দুগুলির মেঘ মানক করা হয়েছে তবে তীরগুলি দৈর্ঘ্যের সাথে বৈচিত্রের সাথে সামঞ্জস্য করে।
biplotসাথে কী উত্পন্ন করে তার সাথে আরও মিলিয়েছে scale=0।