0

私はマトリックスを持っています

結果[i]

いくつかのデータ（分布パラメータ）を使用：

list(structure(c(0.005, 0.004),
.Names = c("mean", "sd")))

例えば。fdistを使用したいので、results[i]に含まれる値を使用して分布パラメーターをシードします。

params = fitdist( data, dist,method="mle",
                    start=list(mean =mapply("[", results[1], 1),
                               sd=mapply("[", results[1], 2)))

次のエラーが発生しました：

the function mle failed to estimate the parameters, 
            with the error code 100

スタートリストは次のとおりです。

structure(list(mean = structure(0.005, .Names = "mean"), 
sd = structure(0.004, .Names = "sd")), .Names = c("mean","sd"))

そしてそれは次のようになります：

structure(list(mean = 0.005, sd = 0.004), .Names = c("mean","sd"))

最後の出力は以下から取得されます。

params = fitdist( data, dist,method="mle",
                    start=list(mean=0.005,
                               sd=0.004))

何か案が？

ありがとう！

4

1 に答える 1

2

「[」の代わりに「[[」を使用してみてください。これは、「[[」がリストノードで値を取得するのに対し、「[」は値をリスト内に残すためです。

res =list(structure(c(0.005, 0.004),
          .Names = c("mean", "sd")))
list(mean =mapply("[[", res, 1),
                                sd=mapply("[[", res, 2))
$mean
[1] 0.005

$sd
[1] 0.004

（私はsapplyを使用したでしょうが。）

> list(mean =sapply( res,"[[", 1),
+                                sd=sapply(res,"[[", 2))
$mean
[1] 0.005

$sd
[1] 0.004

> dput( list(mean =sapply( res,"[[", 1),
+                                sd=sapply(res,"[[", 2)) )
structure(list(mean = 0.005, sd = 0.004), .Names = c("mean", 
"sd"))
于 2012-08-23T17:31:08.867 に答える