measurment関数の技術的なエラーR

Question

私はまだ関数を書くスキルを習得していないので、助けを求めています！

私は、測定誤差（TEM）と2つの列にわたる相対TEMを計算する関数を作成したいと考えています。私はこのことについて、以前の記事を見てきました

は、しかし、彼らは私が彼らが何をする必要があるかやっていないようです...

最初の1：Technical error of measurement (TEM) for 3 or more participantsは3人のオブザーバーのためです。違う結果が出たら、2人のオブザーバーと一緒に試してみてください。

第2のもの：Technical error of measurement in between two columnsは私が探しているものですが、私はそれを関数でありマニュアルコードではないと思います。

誰かが私のデータのための機能にこれを書くのを手伝ってくださいすることができ：相対TEM用

n <- 23 
D <- sum(MHDo1Ana$MHD.o1m1 - MHDo1Ana$MHD.o1m2) 
D2 <- D^2 
TEM <- sqrt((D2)/2*n) 
TEM <- sqrt(((sum(MHDo1Ana$MHD.o1m1 - MHDo1Ana$MHD.o1m2)^2))/2 * n) 

などの機能を？

私のデータは次のとおりです。

MHD.o1m1 MHD.o1m2 
1  46.58 46.77 
2  50.66 50.50 
3  51.15 50.98 
4  46.54 46.16 
5  45.78 45.82 
6  47.83 47.98 
7  43.59 43.66 
8  40.27 40.59 
9  50.08 50.11 
10 47.66 47.69 
11 46.11 46.11 
12 42.90 43.34 
13 38.93 39.14 
14 43.89 44.05 
15 48.14 48.39 
16 46.20 46.11 
17 51.34 51.12 
18 43.48 43.41 
19 46.29 46.14 
20 42.26 42.28 
21 45.33 45.28 
22 47.43 47.23 
23 37.37 37.32 

はありがとうございました！

Answer 1

これは、ジョブ実行する必要があります。

TEM <- function(df, col1, col2){ 
    sqrt(((sum(df[, col1] - df[, col2])^2))/2 * nrow(df)) 
} 

まず、次に、あなたのコンソールでそれを貼り付けて、それを宣言し、それを使用する（データセットに適用）：下記

TEM(MHDo1Ana, "MHD.o1m1", "MHD.o1m2") 

Reproducible example：

MHDo1Ana <- read.table(text=" MHD.o1m1 MHD.o1m2 
1  46.58 46.77 
      2  50.66 50.50 
      3  51.15 50.98 
      4  46.54 46.16 
      5  45.78 45.82 
      6  47.83 47.98 
      7  43.59 43.66 
      8  40.27 40.59 
      9  50.08 50.11 
      10 47.66 47.69 
      11 46.11 46.11 
      12 42.90 43.34 
      13 38.93 39.14 
      14 43.89 44.05 
      15 48.14 48.39 
      16 46.20 46.11 
      17 51.34 51.12 
      18 43.48 43.41 
      19 46.29 46.14 
      20 42.26 42.28 
      21 45.33 45.28 
      22 47.43 47.23 
      23 37.37 37.32") 

TEM(MHDo1Ana, "MHD.o1m1", "MHD.o1m2") 
[1] 1.254731 

これで問題は解決しますか？