3

ここに投稿するのは初めてです。

次の例のようなデータがあります。

C1 C2 C3 C4
1 1.33 1.86 1.91 2.51
2 1.54 1.74 1.98 2.25
3 1.75 1.97 2.20 2.59

行あたりの x 未満の値の平均数が 1 に最も近くなるようなパラメーター値 (x) を見つけたいと考えています。ただし、1,000 の行と 100 の列がある場合は、R を使用して検索したいと考えています。

x = 1.75
dataset.counts <- rowSums(dataset[,] < x)
dataset.counts
[1] 1 2 0
mean(dataset.counts)
[1] 1

optim() を使用して x を見つけることはできますか?

4

2 に答える 2

2

@Glen_b、ありがとう。uniroot を使用してうまく機能したコードのコンソールを投稿しています。ご助力ありがとうございます。

データセット
R1 R2 R3 R4
1 1.33 1.86 1.91 2.51
2 1.54 1.74 1.98 2.25
3 1.75 1.97 2.20 2.59

f <- function(x){
+ 1 - mean(rowSums(dataset[,] < x))
+ }

uniroot(f, c(0,6), tol= 1e-10)
$root
[1] 1.747159

$f.root
[1] 0

$iter
[1] 8

$estim.prec
[1] 0.008522727

于 2013-07-19T11:56:41.103 に答える