我有一个dataframe
,我想申请一个函数,三列的值,并计算三个值之间的最小差值。
#dataset
df <- data.frame(a= sample(1:100, 10),b = sample(1:100, 10),c= sample(1:100, 10))
#function
minimum_distance <- function(a,b,c)
{
dist1 <- abs(a-b)
dist2 <- abs(a-c)
dist3 <- abs(b-c)
return(min(dist1,dist2,dist3))
}
我在寻找的东西,如:
df$distance <- sapply(df, function(x) minimum_distance(x$a,x$b,x$c) )
## errormessage
Error in x$a : $ operator is invalid for atomic vectors
虽然我可以用ddply:
df2 <- ddply(df,.(a),function(r) {data.frame(min_distance=minimum_distance(r$a,r$b, r$c))}, .drop=FALSE)
这不把所有的列。 有什么建议?
编辑:我最终使用:
df$distance <- mapply(minimum_distance, df$a, df$b, df$c)
尝试mapply():
qq <- mapply(minimum_distance, df$a, df$b, df$c)
试试这个:
do.call("mapply", c(list(minimum_distance), df))
但你可以写矢量版本:
pminimum_distance <- function(a,b,c)
{
dist1 <- abs(a-b)
dist2 <- abs(a-c)
dist3 <- abs(b-c)
return(pmin(dist1,dist2,dist3))
}
pminimum_distance(df$a, df$b, df$c)
# or
do.call("pminimum_distance", df)
我知道这已经回答了,但其实我采取了不同的方法,采用任意数量的列,并使用外部方法更普及:
vdiff <- function(x){
y <- outer(x, x, "-")
min(abs(y[lower.tri(y)]))
}
apply(df, 1, vdiff)
我认为这是一个有点清洁和灵活。
编辑:每扎克的意见,我提出这一点,通过删除它们,只对数字列作用于与非数字列数据帧作品,以及更为正式的功能。
cdif <- function(dataframe){
df <- dataframe[, sapply(dataframe, is.numeric)]
vdiff <- function(x){
y <- outer(x, x, "-")
min(abs(y[lower.tri(y)]))
}
return(apply(df, 1, vdiff))
}
#TEST it out
set.seed(10)
(df <- data.frame(a = sample(1:100, 10), b = sample(1:100, 10),
c = sample(1:100, 10), d = LETTERS[1:10]))
cdif(df)
它能够更好地编写一个函数,然后在矢量mapply使用:
f1 <- function(a,b,c){
d =abs(a-b)
e =abs(b-c)
f= abs(c-a)
return(pmin(d,e,f))
}
qq <- mapply(f1, df$a, df$b, df$c)