{"provider_name":"Hatena Blog","author_url":"https://blog.hatena.ne.jp/ryamada22/","type":"rich","categories":["R","\u6700\u5c0f\u4e8c\u4e57\u6cd5"],"provider_url":"https://hatena.blog","published":"2018-05-15 15:41:23","image_url":null,"description":"\u306e\u6700\u5c0f\u5316 \u306e\u6700\u5c0f\u5316 M\u306e\u6700\u5c0f\u56fa\u6709\u5024\u306b\u5bfe\u5fdc\u3059\u308b\u56fa\u6709\u30d9\u30af\u30c8\u30eb\u304c\u4e00\u756a\u62e1\u5927\u7387\u304c\u5c0f\u3055\u3044\u304b\u3089\u3001\u305d\u306e\u56fa\u6709\u30d9\u30af\u30c8\u30eb\u65b9\u5411\u306e\u304c\u6c42\u3081\u308b\u89e3 n <- 100 X <- matrix(rnorm(n*2),ncol=2) Y <- X %*% c(1,2) + rnorm(n,0,0.01) library(rgl) plot3d(X[,1],X[,2],Y) lm(Y ~ X-1) XY <- cbind(X,Y) M <- t(XY) %*% XY eigenout <- eigen(M) eigenout[[2]] eigenout[[2]][,3]/(-eigenout[[2]][3,3])","blog_title":"ryamada\u306e\u907a\u4f1d\u5b66\u30fb\u907a\u4f1d\u7d71\u8a08\u5b66\u30e1\u30e2","height":"190","author_name":"ryamada22","blog_url":"https://ryamada22.hatenablog.jp/","url":"https://ryamada22.hatenablog.jp/entry/20180515/1526366483","width":"100%","version":"1.0","title":"\u56fa\u6709\u5024\u5206\u89e3\u3059\u308b","html":"<iframe src=\"https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada22.hatenablog.jp%2Fentry%2F20180515%2F1526366483\" title=\"\u56fa\u6709\u5024\u5206\u89e3\u3059\u308b - ryamada\u306e\u907a\u4f1d\u5b66\u30fb\u907a\u4f1d\u7d71\u8a08\u5b66\u30e1\u30e2\" class=\"embed-card embed-blogcard\" scrolling=\"no\" frameborder=\"0\" style=\"display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;\"></iframe>"}