R语言一步到位求数据的均值,频数,标准差,标准误差

用R语言快速求数据各组别的均值,频数,标准差,标准误差。

使用R语言中的'plyr'包的ddply()函数,和常见的summarise()函数。

导入'plyr'包和你的数据。

library(plyr)
library(readxl)
NP <- read_excel("C:/Users/lenovo/Desktop/NP.xlsx")
View(NP)
head(NP)#查看数据前6行
str(NP)#查看数据形式

 数据为两组自变量: 添氮即SN(三个浓度水平:N1,N2,N3)、添磷SP(两个浓度水平:LP,HP)

因变量为:植物茎秆鲜重(sfw),植物根系鲜重(rfw)。

使用ddply()函数和summarise()函数求得sfw,rfw在不同处理水平下的的均值,频数,标准差,标准偏差。

zz<-ddply(NP,c("SN","SP"),summarise,psfw=mean(sfw,na.rm=TRUE),sd1=sd(sfw,na.rm = TRUE),n1=sum(!is.na(sfw)),se1=sd1/sqrt(n1),prfw=mean(rfw,na.rm=TRUE),sd2=sd(rfw,na.rm = TRUE),n2=sum(!is.na(rfw)),se2=sd2/sqrt(n2))
zz

这样就可以算出各组在不同水平下的sfw、rfw的均值(psfw,prfw)、标准差(sd1,sd2)、频数(n1,n2)、标准误差(se1,se2)。

 说明:na.rm=TURE为忽略缺失数据的意思,sum(!is.na())为计算非缺失数据的频数。

数据百度云链接:https://pan.baidu.com/s/1KpjncD26-MM_x3_E1HQ00A 
提取码:vvv2

  • 16
    点赞
  • 77
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
R语言,根据统计数据创建频数直方图(也称为频率分布直方图)是非常常见的数据可视化步骤。以下是如何使用`ggplot2`包来创建一个基本的频数直方图的步骤: 1. 首先,确保你已经安装了`ggplot2`和`tidyverse`(如果尚未安装,可以使用`install.packages("ggplot2")`和`install.packages("tidyverse")`),这两个包提供了强大的绘图功能。 2. 准备数据数据可以是一个向量,包含你要分析的数值。例如,假设我们有一个名为`data`的变量,它包含了样本数据。 ```R data <- c(2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30) # 示例数据 ``` 3. 使用`ggplot`创建基础图形: ```R library(ggplot2) library(tidyverse) # 创建数据框,方便ggplot处理 hist_data <- data.frame(value = data) # 创建直方图 hist_plot <- ggplot(hist_data, aes(x = value)) + geom_histogram(binwidth = 2, fill = "lightblue", color = "black") + labs(title = "Frequency Histogram", x = "Value", y = "Frequency") ``` 4. `binwidth`参数定义了每个柱子(或“箱”)的宽度,你可以根据数据的性质调整这个值。这里设置了`binwidth = 2`,表示每个箱间隔为2。 5. `fill`和`color`用于设置柱子的颜色。 6. `labs()`函数用于添加图表标题和轴标签。 现在你已经有了一个基本的频数直方图。如果你的数据更复杂,比如来自数据框或数据集,只需要将数据源换成相应的变量名即可。
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值