请原谅我的问题很奇怪。这是最合理的框架。真的,我不在乎花朵。
我最终想做的是找到在花瓣宽度/花瓣长度的每个“临界值”以上的鸢尾花的萼片长度的总和。
考虑以下代码
library(tidyverse)
data("iris")
iris <- iris %>%
mutate(prop_width_length = Petal.Width/Petal.Length)
prop_width_length <- as.data.frame(iris$prop_width_length)
portion = as.data.frame(seq(0,1,0.001))
cumsum = NULL
for (i in 1:1001) {
cumsum[i] = sum(prop_width_length >= portion[i,1])
}
sigportion <- cbind(portion, cumsum)
这使我累积了多少个鸢尾花的宽度/长度比例大于或等于每个“临界值”的总和。然后最后将其放入数据帧中,这样我就可以制作一个漂亮的ggplot。基本上,它计算每个“临界值”上有多少朵花。
除了上面的代码外,我还想为每个虹膜的花瓣长度/长度比大于或等于存储在我的部分变量中的每个“临界值”求和所有的萼片长度。
类似
sum all the sepal lengths of iris flowers which have petal width/length >= critvalue
答案 0 :(得分:1)
您可以使用lapply
进行操作。一定有一种更整洁的方法,但这是我的版本:
mySeq<-seq(0,1,0.001)
df<-data.frame(mySeq=mySeq,cumsum=as.numeric(lapply(mySeq,function(X)sum(ifelse(iris$Petal.Width/iris$Petal.Length>=X,1,0)))))
在功能中,我正在检查条件并总结满足的情况
下面的代码应为您提供隔片长度的总和
dfSepalLengths<-data.frame(mySeq=mySeq,cumsum=as.numeric(lapply(mySeq,function(X)sum(ifelse(iris$Petal.Width/iris$Petal.Length>=X,iris$Sepal.Length,0)))))
您还可以如下修改代码以得到相同的结果:
library(tidyverse)
data("iris")
iris <- iris %>%
mutate(prop_width_length = Petal.Width/Petal.Length)
portion = as.data.frame(seq(0,1,0.001))
cumsum = NULL
for (i in 1:1001) {
cumsum[i] = sum(ifelse(iris$prop_width_length >= portion[i,1],iris$Sepal.Length,0))
}
sigportion <- cbind(portion, cumsum)
答案 1 :(得分:1)
使用data.table
library(data.table)
iris<-as.data.table(iris)
iris[,prop_width_length := Petal.Width/Petal.Length]
portion<-as.data.table(seq(from = 0,to = 1,by = 0.001))
cumsum<-vector()
for(i in 1:nrow(portion)){
cumsum[i]<-iris[prop_width_length >= portion[[1]][i],sum(Sepal.Length)]
}
sigportion<-cbind(portion,cumsum)
希望有帮助!