Question

我有一个数据框“蠕虫”，我对计算平均值有两个感兴趣的因素-长度相同的“ worm.density”和“ Vegetation”。植被有5个数据框，我平均需要2个-“ Meadow”和“ Grassland”，同时还要插入一个条件，该条件仅包含Area> 2.5和Soil.Ph> 3.5因子的结果。是否可以通过轻触来做到这一点，或者循环是实现此目的的最佳方法？我卡住了，只能用我的植被因子限制两个数据帧。

从www.bio.ic.ac.uk/research/crawley/statistics /

下载的

数据文件

worms<-structure(list(Field.Name = structure(c(8L, 17L, 10L, 16L, 7L, 
11L, 3L, 1L, 19L, 15L, 5L, 9L, 18L, 12L, 13L, 20L, 2L, 14L, 6L, 
4L), .Label = c("Ashurst", "Cheapside", "Church.Field", "Farm.Wood", 
"Garden.Wood", "Gravel.Pit", "Gunness.Thicket", "Nashs.Field", 
"North.Gravel", "Nursery.Field", "Oak.Mead", "Observatory.Ridge", 
"Pond.Field", "Pound.Hill", "Rookery.Slope", "Rush.Meadow", "Silwood.Bottom", 
"South.Gravel", "The.Orchard", "Water.Meadow"), class = "factor"), 
    Area = c(3.6, 5.1, 2.8, 2.4, 3.8, 3.1, 3.5, 2.1, 1.9, 1.5, 
    2.9, 3.3, 3.7, 1.8, 4.1, 3.9, 2.2, 4.4, 2.9, 0.8), Slope = c(11L, 
    2L, 3L, 5L, 0L, 2L, 3L, 0L, 0L, 4L, 10L, 1L, 2L, 6L, 0L, 
    0L, 8L, 2L, 1L, 10L), Vegetation = structure(c(2L, 1L, 2L, 
    3L, 5L, 2L, 2L, 1L, 4L, 2L, 5L, 2L, 2L, 2L, 3L, 3L, 5L, 1L, 
    2L, 5L), .Label = c("Arable", "Grassland", "Meadow", "Orchard", 
    "Scrub"), class = "factor"), Soil.pH = c(4.1, 5.2, 4.3, 4.9, 
    4.2, 3.9, 4.2, 4.8, 5.7, 5, 5.2, 4.1, 4, 3.8, 5, 4.9, 4.7, 
    4.5, 3.5, 5.1), Damp = c(FALSE, FALSE, FALSE, TRUE, FALSE, 
    FALSE, FALSE, FALSE, FALSE, TRUE, FALSE, FALSE, FALSE, FALSE, 
    TRUE, TRUE, TRUE, FALSE, FALSE, TRUE), Worm.density = c(4L, 
    7L, 2L, 5L, 6L, 2L, 3L, 4L, 9L, 7L, 8L, 1L, 2L, 0L, 6L, 8L, 
    4L, 5L, 1L, 3L)), class = "data.frame", row.names = c(NA, 
-20L))

> with(worms,tapply(Worm.density,list[Grassland,Meadow],mean))

sessioninfo（）R版本3.5.1（2018-07-02）平台：x86_64-apple-darwin15.6.0（64位）运行于：macOS Sierra 10.12.6

Matrix产品：默认 BLAS：/System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib LAPACK：/System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libLAPACK.dylib

语言环境： [1] C

附带的基本软件包： [1]统计图形grDevices utils数据集方法基础

其他附件包： [1] RevoUtils_11.0.1

通过名称空间（未附加）加载： [1]编译器_3.5.1工具_3.5.1

Answer 1

简单地subset tapply中的{em>蠕虫数据帧就会产生一个命名向量：

with(subset(worms, Area > 2.5 & Soil.pH > 3.5), 
            tapply(Worm.density, Vegetation, mean)
     )

#   Arable Grassland    Meadow   Orchard     Scrub 
# 6.000000  2.333333  7.000000        NA  7.000000

要仅返回草原和草地，请用[]

为项目建立索引

with(subset(worms, Area > 2.5 & Soil.pH > 3.5), 
            tapply(Worm.density, Vegetation, mean)
     )[c("Grassland", "Meadow")]

# Grassland    Meadow 
#  2.333333  7.000000

仅对于 Grassland 和 Meadow 类型均值，将其添加到子集中并取均值 Worm.density ：

sub_worms <- subset(worms, Area > 2.5 & Soil.pH > 3.5 &
                           Vegetation %in% c("Grassland", "Meadow"))
sub_worms
#       Field.Name Area Slope Vegetation Soil.pH  Damp Worm.density
# 1    Nashs.Field  3.6    11  Grassland     4.1 FALSE            4
# 3  Nursery.Field  2.8     3  Grassland     4.3 FALSE            2
# 6       Oak.Mead  3.1     2  Grassland     3.9 FALSE            2
# 7   Church.Field  3.5     3  Grassland     4.2 FALSE            3
# 12  North.Gravel  3.3     1  Grassland     4.1 FALSE            1
# 13  South.Gravel  3.7     2  Grassland     4.0 FALSE            2
# 15    Pond.Field  4.1     0     Meadow     5.0  TRUE            6
# 16  Water.Meadow  3.9     0     Meadow     4.9  TRUE            8

mean(sub_worms$Worm.density)
# [1] 3.5

Rextester demo

轻敲vs循环以获取因子子集和条件

1 个答案: