请教R编程问题，谢谢各位大侠！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

有一个6300*240的矩阵
要求算它的6300行的spatial depth
大概运算过程是这样的
1.
U[ij]=(Xj-Xi)/||Xj-Xi||, i,jЄ{1,...,n},i不等于j
2.
Vi=(1/(n-1))∑U[ij] for j from 1 to n, j不等于i
3.
The Sample depth at Xi is obtained as
D(Xi, Fn)=1-||Vi||
下面是我写的循环
大概运行到U[j,]=(T2-T1)/(sum((T2-T1)*(T2-T1))^0.5) 这步就显示阵里的下标数目不对
本人R初学者，希望各位大侠不吝赐教
若嫌问题过于简单，还望海涵！
library(Matrix)
DataSD=function(Data){
TimeSeries=Data
U=matrix(nrow=6300,ncol=240)
V=matrix(nrow=6300,ncol=240)
D=matrix(nrow=6300,ncol=1)
for(i in 1:nrow(TimeSeries)){
for(j in 1:nrow(TimeSeries)){
if(i==j){
U[j,]=rep(0,length=240)}
else{
T1=TimeSeries[i,2:ncol(TimeSeries)]
T2=TimeSeries[j,2:ncol(TimeSeries)]
U[j,]=(T2-T1)/(sum((T2-T1)*(T2-T1))^0.5)
}
}
V[i,]=(1/(n-1))*colSums(U$Vector[i,],na.rm=TRUE)
D[i,]=1-(sum(V[i,]*V[i,]))^0.5
}
return(list(U=U,V=V,D=D))
}

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：各位大侠 R编程 timeseries obtained function 请教 大侠

你的TimeSeries是240列
T1和T2都是239列，然后你又赋值给240列的U，所以就出错了。

我的数据Timeseires第一列是location所以是241行的 2# qoiqpwqr

library(foreach)
A <- 63
B <- 24
TimeSeries <- abs(matrix(rnorm(A*B),A,B)*10)+10
UU <- foreach(T2=iapply(TimeSeries,1),.combine="rbind") %:% foreach(T1=iapply(TimeSeries,1),.combine="rbind") %do% {
  T2-T1
}
U <- foreach(a=iapply(UU,1),.combine="rbind") %do% {
  a/crossprod(a)^0.5
}
K <- seq(1,A * A,A)
V <- foreach(a=K,.combine="rbind") %do% {
   colSums(U[a:(a+A-1),],na.rm=TRUE)
}
D <- 1-apply(V,1,function(x) crossprod(x)^0.5)