##Read parameters set data 读取参数数据,参数数据是10000行*2列数据
ParamData <- read.csv("Coeffi.csv", header = TRUE)
##Read choice set data 读取样本参数,样本参数是283行*12列数据
Data<-read.csv("Sequence.csv",header=TRUE)
想要做的是参数数据每一行分别乘以样本数据某一列所有行的数据
因为是10000*1参数矩阵 乘以283*1样本数据矩阵,所以把样本数据矩阵转置了
##Transfer data for easily calculation 转置样本数据矩阵
TData <- t(Data)
这样就可以变成10000*1参数矩阵乘以1*283样本数据得到10000*283矩阵
##calculate utility of each alternative
d <- matrix(ParamData[,1]*TData[3,] + ParamData[,2]*TData[8,], nrow = 10000, ncol = 283)
e <- matrix(ParamData[,1]*TData[4,] + ParamData[,2]*TData[9,], nrow = 10000, ncol = 283)
f <- matrix(ParamData[,1]*TData[5,] + ParamData[,2]*TData[10,], nrow = 10000, ncol = 283)
g <- matrix(ParamData[,1]*TData[6,] + ParamData[,2]*TData[11,], nrow = 10000, ncol = 283)
h <- matrix(ParamData[,1]*TData[7,] + ParamData[,2]*TData[12,], nrow = 10000, ncol = 283)
print(e)
问题:发现e矩阵每一行都比前一行结果大,貌似结果逐行累加了,求破解!哪错了?
万分感谢~