有機(jī)酸摩爾質(zhì)量的測(cè)定

有機(jī)酸摩爾質(zhì)量的測(cè)定有機(jī)酸摩爾質(zhì)量測(cè)定實(shí)驗(yàn)的數(shù)據(jù)處理，分為兩部分：1/標(biāo)準(zhǔn)溶液的標(biāo)定2/有機(jī)酸摩爾質(zhì)量的測(cè)定計(jì)算公式：1/

為鄰苯二甲酸氫鉀的質(zhì)量，單位：2/

為氫氧化鈉標(biāo)準(zhǔn)溶液的體積，單位：3/

為空白實(shí)驗(yàn)氫氧化鈉溶液的用量，單位：4/

為每

的質(zhì)量數(shù)據(jù)記錄表格如上所示，因此，需要完成數(shù)據(jù)的快速計(jì)算。設(shè)計(jì)腳本：m

<-

c()

v

<-

c()

c1

<-

m

/

v

/

0.2042

c1

c1

%>%

mean

c1

%>%

rd

c1

%>%

rd_mean

c1

%>%

sd

c1

%>%

rsd

t.test(c1,

conf.level

=

0.9)為各次稱量的基準(zhǔn)物質(zhì)的質(zhì)量為

溶液所用的體積，暫且不考慮

的影響對(duì)應(yīng)濃度的計(jì)算公式，由此得出濃度值用于計(jì)算平均值，這里使用管道操作符，來(lái)更加簡(jiǎn)明地表示為相對(duì)偏差的計(jì)算為相對(duì)平均偏差的計(jì)算計(jì)算標(biāo)準(zhǔn)偏差計(jì)算相對(duì)標(biāo)準(zhǔn)偏差完成t檢驗(yàn)，輸出

置信水平的置信區(qū)間說(shuō)明：1/需要自行編寫2/

需要自行編寫3/

也是

語(yǔ)言默認(rèn)不存在的函數(shù)示例：>

m

<-

c(0.4803,0.4823,0.4813)

>

v

<-

c(

21.20,21.32,21.24)

>

c1

<-

m

/

v

/

0.2042

>

c1

[1]

0.11095

0.11078

0.11097

>

c1

%>%

mean

[1]

0.1109

>

c1

%>%

rd

[1]

0.0431%

-0.1057%

0.0626%

>

c1

%>%

rd_mean

[1]

0.6356%

>

c1

%>%

sd

[1]

0.00010212

>

c1

%>%

rsd

[1]

0.09209%

>

t.test(c1,

conf.level

=

0.9)

One

Sample

t-test

data:

c1

t

=

1881,

df

=

2,

p-value

=

2.8e-07

alternative

hypothesis:

true

mean

is

not

equal

to

0

90

percent

confidence

interval:

0.11073

0.11107

sample

estimates:

mean

of

x

0.1109數(shù)據(jù)不具有真實(shí)性，僅用于驗(yàn)證腳本可行性。編寫函數(shù)實(shí)際效果：>

x

<-

c(0.1104,0.1098,0.1114)

>

x

[1]

0.110

0.110

0.111

>

rd(x)

[1]

-0.1206%

-0.6634%

0.7841%

>

rd(x,3)

[1]

-0.121%

-0.663%

0.784%

>

rd(x,2)

[1]

-0.12%

-0.66%

0.78%需要兩個(gè)參數(shù)：待計(jì)算的數(shù)據(jù)，以及控制輸出的位數(shù)源代碼：rd

<-

function(x,digits

=

4)

{

x1

<-

mean(x)

t

=

(x

-

x1)/

x1

paste0('%0.0',digits,'f')

%>%

sprintf(t*100)

%>%

paste0('%')

%>%

print(quote

=

F)

}位數(shù)的輸出控制，由sprintf()完成，同時(shí)傳遞給paste0()函數(shù)，對(duì)數(shù)據(jù)添加

%,實(shí)際輸出的內(nèi)容僅僅是字符串！而quote=F，則是去除引號(hào)，使得最終輸出的內(nèi)容更像是

百分?jǐn)?shù)表示的數(shù)字rd_mean():>

rd_mean(x)

[1]

4.729%

>

rd_mean(x,3)

[1]

4.73%

>

rd_mean(x,2)

[1]

4.7%則類似于源代碼：rd_mean

<-

function(x,

digits

=

4){

x1

<-

mean(x)

t

=

(x

-

x1)/

x1

abs(t)

%>%

mean

->

d

d

/

x1

*

100

->

f

format(f,

digits

=

digits)

%>%

paste0(

'%')

%>%

print(

quote

=

F)

}相對(duì)應(yīng)的計(jì)算公式：

相對(duì)而言，format()的格式控制，比sprintf()更加簡(jiǎn)單，但對(duì)于更復(fù)雜的格式輸出，仍建議使用

sprintf()rsd:>

rsd(x)

[1]

0.7313%

>

rsd(x,3)

[1]

0.731%

>

rsd(x,2)

[1]

0.73%源代碼：rsd

<-

function(x,digits

=

4){

sd(x)

/

mean(x)

*

100

->

a

format(a,

digits

=

digits)

%>%

paste0(

'%')

%>%

print(

quote

=

F)

}計(jì)算公式：t分布>

t.test(c1,

conf.level

=

0.9)

One

Sample

t-test

data:

c1

t

=

1881,

df

=

2,

p-value

=

2.8e-07

alternative

hypothesis:

true

mean

is

not

equal

to

0

90

percent

confidence

interval:

0.11073

0.11107

sample

estimates:

mean

of

x

0.1109利用t.test()返回

置信水平的置信區(qū)間反應(yīng)式：計(jì)算方法：c2

<-

mean(c1)

m2

<-

c()

v2

<-

25

n

<-

v3

<-

c()

M

<-

m2

/

(c2

*

v3

*

0.001

/

n

*10)

M

%>%

mean

M

%>%

rd

M

%>%

rd_mean為測(cè)定的

溶液的濃度的平均值為有機(jī)酸的質(zhì)量為移取的有機(jī)酸溶液的體積：為多元酸對(duì)應(yīng)的

值為

溶液的體積，單位：為計(jì)算所得的有機(jī)酸的摩爾質(zhì)量示例：>

c2

<-

0.1002

>

m2

<-

c(0.4905,0.4906,0.4910)

#

有機(jī)酸質(zhì)量

>

v2

<-

25

#

有機(jī)酸體積

>

n

<-

2

#

多元酸

n

>

v3

<-

c(18.00,18.02,18.06)

#

NaOH

體積

>

M

<-

m2

/

(c2

*

v3

*

0.001

/

n

*

10)

>

M

[1]

54.391

54.342

54.266

>

M

%>%

mean

[1]

54.333

>

M

%>%

rd

[1]

0.1072%

0.0165%

-0.1237%

>

M

%>%

rd_mean

[1]

0.001518%數(shù)據(jù)僅用于驗(yàn)證腳本的可行性，不具有實(shí)際用途！Q檢驗(yàn)一組數(shù)據(jù)，按從小到大排列：

設(shè)

為可疑值，根據(jù)

確定可疑值的取舍：假設(shè)

為可疑值，

假設(shè)

為可疑值，則

通過(guò)與

表

進(jìn)行比對(duì)，若

值大于表中的值，則舍去該可疑值。計(jì)算方法：Q.test

<-

function(x)

{

x

<-

sort(x)

n

<-

length(x)

x1

<-

diff(x)

x1

/

(x[n]

-

x[1])

->

f

(x[2]

-

x[1]

)

/

(x[n]

-

x[1])

->

g

appen