Ipfp函数 - 如何指定target.list参数

时间:2018-03-14 20:47:29

标签: r

如果有人明白下面粗体句的含义,请帮助我理解这意味着什么?它来自Ipfp函数的R文档。 (https://www.rdocumentation.org/packages/mipfp/versions/3.1/topics/Ipfp

" target.list - target.data中边际目标约束的维度列表。 列表的每个组件都是一个数组,其单元格指示相应边距与哪个维度相关。

" target.data - 包含目标边际表数据的列表。列表的每个组件都是一个存储边距的数组。列表顺序必须遵循target.list中定义的顺序。请注意,数组的单元格必须是非负数。"

作为一个例子,让我们说我有这张表:

seed.ex <- array(1,dim=c(3,4))
seed.ex
     [,1] [,2] [,3] [,4]
[1,]    1    1    1    1
[2,]    1    1    1    1
[3,]    1    1    1    1

这些目标存储在&#39; target.data&#39;:

target.row <- c(50,50,100)
target.col <- c(50,50,50,50) 
tgt.data.ex <- list(target.row, target.col) 

我应该如何指定&#39; target.list&#39;?

tgt.list.ex <- list(?,?)

那将进入Ipfp功能......

res.ex <- Ipfp(seed.ex, tgt.list.ex, tgt.data.ex, print = TRUE, iter = 1000)

1 个答案:

答案 0 :(得分:0)

我非常确定答案是list(1,2)。我从来没有使用过这个软件包,但是使用二维种子的两个例子都有这个答案。这是我的考试:

tgt.list.ex <- list(1,2)
also installing the dependency ‘cmm’

trying URL 'https://cloud.r-project.org/bin/macosx/el-capitan/contrib/3.4/cmm_0.8.tgz'
Content type 'application/x-gzip' length 197741 bytes (193 KB)
==================================================
downloaded 193 KB

trying URL 'https://cloud.r-project.org/bin/macosx/el-capitan/contrib/3.4/mipfp_3.1.tgz'
Content type 'application/x-gzip' length 333155 bytes (325 KB)
==================================================
downloaded 325 KB


The downloaded binary packages are in
    /var/folders/yq/m3j1jqtj6hq6s5mq_v0jn3s80000gn/T//Rtmp0NHfpB/downloaded_packages
> library(mipfp)
Loading required package: cmm
Loading required package: Rsolnp
Loading required package: numDeriv
> res.ex <- Ipfp(seed.ex, tgt.list.ex, tgt.data.ex, print = TRUE, iter = 1000)
Margins consistency checked!
... ITER 1 
       stoping criterion: 24 
... ITER 2 
       stoping criterion: 0 
Convergence reached after 2 iterations!
> res.ex

Call:
Ipfp(seed = seed.ex, target.list = tgt.list.ex, target.data = tgt.data.ex, 
    print = TRUE, iter = 1000)

Method:  ipfp - convergence:  TRUE 

Estimates:
     [,1] [,2] [,3] [,4]
[1,] 12.5 12.5 12.5 12.5
[2,] 12.5 12.5 12.5 12.5
[3,] 25.0 25.0 25.0 25.0