Yn(1)=[]; %从第二个数开始,即x(2),x(3)... Yn=Yn'; %Yn=A(:,2:m)'; Z=Z1(:,2:m)'; B=[-Z,X]; C=(inv(B'*B))*B'*Yn;%由公式建立GM(1,n)模型 a=C(1); C1=C'; b=C1(:,2:n+1); F=[]; F(1)=A(1); u=zeros(1,m); for i=1:m for j=1:n u(i)...
首先我们令测试的向量命为test_vector 如上图所示,我们将十个数值型的数据输入至向量test_vector中,分别为1-10共10个数,其中符号<-为数据赋值符号,这一点区别于其他编程语言的=号,如下所示,我们可以发现10个数值型数据已经完成了输入并正确显示: 接下来我们完成字符型数据的向量输入,步骤与语句和数值型数据输入类...
倾向性得分匹配是以倾向性得分为作为匹配条件,将倾向性得分相同或者相近的研究对象按照1:1或者1:N的比例匹配在一起,从而实现在原始数据库的基础上获得一个倾向性得分基本相近的子数据库,也就是基线基本可比的数据库。 匹配后的数据库在一定程度上控制了混杂,能获得更加真实可靠的结论。本讲我们通过两种软件SPSS和R...
首先为了建立数据集,引入生成因子水平的函数g1(),其调用格式为: gl(n, k, length=n*k,labels=1:n,ordered=FALSE) n是因子的水平个数;k表示每一水平上的重复次数;length=n*k表示总观测数;可通过参数labels对因子的不同水平添加标签;ordered为逻辑值,指示是否排序。 > x=c(20,12,20,10,14,22,10,20,...
1. 跨期相关系数计算原始代码 cor.time=function(x,y,n){ ###x is matrix and y is vector,and the lag term is y### ###n is total periods of lag, while first half is the LAG### ###row of x must be variates, column of x must be time### ...
sum(x)给出了x中所有元素的总和; prod(x)给出x中所有元素的乘积; mean(x)和var(x),分别计算样本均值和样本方差,这两个函数分别相当于sum(x)/length(x),sum((x-mean(x)) \^2)/(length(x) -1)。如果var()的参数是一个n*p的矩阵,那么函数的值是一个p*p的样本协方差矩阵,认为每行是一个p变量的...
R语言中，如何用循环函数生成一个五阶的Hilbert矩阵？输入如下命令即可得到五阶的Hilbert矩阵，以此类堆，若为k阶则将5换为相应数字即可。> n<-5;x<-array(0,dim=c(n,n))> for(i in 1:n){ + for(j in 1:n){ + x[i,j]<-1/(i+j-1)+ } + } > x 如图所示 ...
(1) 函数cat()输出到屏幕 x = "I love R!" y = "R is very powerful!" cat(x, y) cat(x, "\n", y) (2) 函数cat()输出到文件 cat("abc\n", file="./u.txt") cat("de\n", file="./u.txt", append=TRUE) 第一次调用函数cat()创建了文件u,包含一行内容abc;第二次调用追加了第...
cells <-(1,2,6,8) : 向量中有4个元素,分别为1,2,6,8 dim(x):获取向量的维度 dim(x)[1] :获取向量的行数 dim(x)[2]: 获取向量的列数 nrow (x):获取向量的行数 ncol(x): 获取向量的列数 03 04 6.转置矩阵: R语言矩阵提供了 t()函数,可以实现矩阵的行列转换。例如:可以把m行n列转换成...
mean(x)和var(x),分别计算样本均值和样本方差,这两个函数分别相当于sum(x)/length(x),sum((x-mean(x)) \^2)/(length(x) -1)。如果var()的参数是一个n*p的矩阵,那么函数的值是一个p*p的样本协方差矩阵,认为每行是一个p变量的样本向量。