R语言
方月一 +1
方月皆半。
展开
-
矩阵归一化
矩阵归一化实际上就是向量的归一化,就是把向量的模变换为1,用公式表示即u∣∣u∣∣\frac{u}{||u||}∣∣u∣∣u。追究到基本的概念,实际上是个很简单的问题。...原创 2018-12-08 22:09:04 · 4889 阅读 · 0 评论 -
svd奇异值分解
svd奇异值分解R语言svdhilbert <- function(n) { i <- 1:n; 1 / outer(i - 1, i, "+") }X <- hilbert(9)[, 1:6](s <- svd(X))D <- diag(s$d)s$u %*% D %*% t(s$v) # X = U D V' 其中,%*% 表示矩阵相乘t(s...原创 2018-12-06 15:24:31 · 236 阅读 · 0 评论 -
齐一性检验
齐一性检验n1=2 #行n2=3 #列n=matrix(c(58,40,40,40,12,20),n1,n2) #按列写入 #63题nchi=matrix(0,n1,n2)for(i in 1:n1){ for(j in 1:n2){ chi[i,j]=(n[i,j]-sum(n[i,])*sum(n[,j])/sum(n))^2/(sum(n[i,])*sum...原创 2018-12-13 20:12:51 · 620 阅读 · 0 评论 -
R易错
诸如下面这种:dat = dat[-which(dat$ZDJ>1000),] 若which(dat$ZDJ>1000)的运行结果是integer(0)的话,程序会默认把所有行都删了。。。。。。。。原创 2018-12-13 20:29:10 · 118 阅读 · 0 评论 -
pca in R
用R写pca:x=data.frame(x1=c(1,0,6,2,7,8),x2=c(2,5,7,3,6,7),x3=c(3,6,7,2,9,0))x=scale(x,center = T,scale = T)xcov=cov(x)xspe = eigen(xcov) xval = xspe$valuesxvec = xspe$vecrval=xval[-3]rvec=xvec[...原创 2018-12-06 16:39:29 · 415 阅读 · 0 评论 -
R语言计算标准差
R语言计算标准差:sd(x)上面计算的是样本标准差s=1n−1∑i=1nxi−xˉs=\frac{1}{n-1}\sum_{i=1}^{n}x_i-\bar{x}s=n−11∑i=1nxi−xˉ.原创 2018-12-09 22:09:09 · 39055 阅读 · 3 评论