主题：求助：R中如何实现矩阵的开方

Python065

主题：求助：R中如何实现矩阵的开方,第1张

一种方法是先化到Jordan型，Jordan块 3 1 0 3 的一个平方根是 3^{1/2} 12^{-1/2} 0 3^{1/2} 另一种方法是假定B=uA+vI，B^2=A，然后利用Cayley－Hamilton定理A^2-6A+9解出待定系数u和v

在我们实际建模拟合时，若数据中某个变量过度离散，则需要做转化，常见的有做log和根号两种方式。 install.packages('gridExtra')  library(gridExtra)  #先画一个柱状图，观察friend_count这个变量的整体趋势，发现过度离散 p1<-ggplot(aes(x=friend_count),data=pf)+geom_histogram() #log转化 p2<-p1+scale_x_log10() #根号转化 p3<-p1+scale_x_sqrt() #将三张图放在一个面板里 grid.arrange(p1,p2,p3,ncol=1)最后我们发现，做了log转化的数据更接近正态分布。