我想拍摄嘈杂的2d数据并使用LOESS平滑来找到更平滑的2D表面。数据是r行和c列的数据帧或矩阵的形式。 ' asbio'包将为n * n数据做表面图,我可以弄清楚如何用1D拟合'黄土'从统计学基础包,但我不知道如何做2D图。
library(rgl)
age <- 1:100
year <- 1959:2012
z<- data.frame(matrix(array(0, length(age)*length(year)),nrow=length(age), ncol=length(year)))
rownames(z) <- age
colnames(z) <- year
distYear <- dnorm(1:length(year), 26, 10)
distAge <- dnorm(age, 50, 15)
for(theYear in year){
#for(theAge in age){
z[,theYear-year[1]+1] <- as.numeric(distYear[theYear-year[1]+1]*distAge*runif(length(age))*50000)
#}
}
z<- data.matrix(z)
z <- data.frame(z)
colnames(z) <- year
# First define the colours for the Z values.
zlim <- range(z)
zlen <- zlim[2] - zlim[1]
scale <- 20/zlen
colorlut <- terrain.colors(20,alpha=1.0) # height color lookup table
col <- colorlut[ (t(as.matrix(z))-zlim[1])*scale+1 ] # assign colors to heights for each point.
surface3d(age,year,as.matrix(z), col=col)
loess(z~age+year,data=data.frame(z))
我收到一个错误,拒绝将data.frame作为列表,但我怀疑它还有更多。搜索信息时,我只能找到一维线性描述。
有人可以帮忙吗?
答案 0 :(得分:3)
解决了....
需要将数据帧或矩阵转换为三列数据帧。 创建上面的数据框z后,以下内容将平滑并显示它(使用rgl)。
open3d(useFreeType=par3d("useFreeType"))
surface3d(age,year,as.matrix(z), col=col)
zz <- z
zz$Ages <- rownames(z)
zzx <- melt(zz, id.vars="Ages", value.name="Vals")
theResult <- loess(Vals~as.numeric(Ages)+as.numeric(variable),data=data.frame(zzz))
resultTable <- matrix(theResult$fitted, nrow=length(age), length(year))
open3d(useFreeType=par3d("useFreeType"))
surface3d(age,year,as.matrix(resultTable), col=col)