如何在R中建立对数分配箱?

时间:2020-03-03 04:37:27

标签: r

我是R和数学的初学者。

我们正在运行一种设备,该设备可以按粒度测量细小的灰尘。仪器使用16个通道分析0.3-10 um。每个通道根据对数正态分布进行划分,如下表所示。

          Lower         Upper
1           0.3          0.374
2           0.374        0.465
3           0.465        0.579
4           0.579        0.721
5           0.721        0.897
6           0.897        1.117
7           1.117        1.391
8           1.391        1.732
9           1.732        2.156
10          2.156        2.685
11          2.685        3.343
12          3.343        4.162
13          4.162        5.182
14          5.182        6.451
15          6.451        8.032
16          8.032        10

我想再次将表格调整为0.3-2.5微米。如何均匀设置间隔?我尝试使用R编写代码,但这非常困难。

          Lower         Upper
1           0.3          
2                   
3                   
4                   
5                   
6                   
7                   
8                   
9                   
10                  
11                  
12                  
13                  
14                  
15                  
16                       2.5

假定图形形状为对数正态。

1 个答案:

答案 0 :(得分:1)

您可以使用seqlogexp创建以日志分隔的表

tt <- exp(seq(log(0.3), log(2.5), length.out=17))
data.frame(Lower=tt[-length(tt)], Upper=tt[-1])
#       Lower     Upper
#1  0.3000000 0.3425093
#2  0.3425093 0.3910422
#3  0.3910422 0.4464520
#4  0.4464520 0.5097133
#5  0.5097133 0.5819385
#6  0.5819385 0.6643980
#7  0.6643980 0.7585417
#8  0.7585417 0.8660254
#9  0.8660254 0.9887393
#10 0.9887393 1.1288415
#11 1.1288415 1.2887959
#12 1.2887959 1.4714155
#13 1.4714155 1.6799118
#14 1.6799118 1.9179517
#15 1.9179517 2.1897213
#16 2.1897213 2.5000000

或范围0.3到10:

tt <- exp(seq(log(0.3), log(10), length.out=17))
data.frame(Lower=tt[-length(tt)], Upper=tt[-1])
#       Lower      Upper
#1  0.3000000  0.3735091
#2  0.3735091  0.4650301
#3  0.4650301  0.5789766
#4  0.5789766  0.7208434
#5  0.7208434  0.8974719
#6  0.8974719  1.1173797
#7  1.1173797  1.3911716
#8  1.3911716  1.7320508
#9  1.7320508  2.1564557
#10 2.1564557  2.6848527
#11 2.6848527  3.3427230
#12 3.3427230  4.1617915
#13 4.1617915  5.1815565
#14 5.1815565  6.4511948
#15 6.4511948  8.0319330
#16 8.0319330 10.0000000