小编Tom*_*Tom的帖子

如何从GAM中提取拟合样条(`mgcv :: gam`)

我正在使用GAM来模拟逻辑回归中的时间趋势.然而,我想从中提取拟合样条曲线,将其添加到另一个模型中,不能用于GAM或GAMM.

因此我有两个问题:

  1. 随着时间的推移,我怎样才能更顺畅,以便在让模特找到其他结点的同时强迫一个结处于特定位置?

  2. 如何从拟合的GAM中提取矩阵,以便我可以将其用作不同模型的推算?

我正在运行的模型类型如下:

gam <- gam(mortality.under.2~ maternal_age_c+ I(maternal_age_c^2)+
           s(birth_year,by=wealth2) + wealth2 + sex +
           residence + maternal_educ + birth_order,
           data=colombia2, family="binomial")
Run Code Online (Sandbox Code Playgroud)

我已经阅读了GAM的大量文档,但我还不确定.任何建议都非常感谢.

r gam mgcv

13
推荐指数
1
解决办法
9229
查看次数

使用mac中的终端跨目录复制文件

如何将具有特定扩展名的所有文件.SAV从一个目录(包括该目录中的所有子目录)复制到另一个目录?我正在使用mac的终端.我试过了

    /Users/tournillon/my_directory/research_projects/dissertation_related/inequality
    Antonio-P-Tournillon-Ramoss-iMac:inequality tournillon$ ls
    Brazil              dhs_dados           imr201101app.pdf
    GHME_Education_final.pptx   dhs_sav_files
    Antonio-P-Tournillon-Ramoss-iMac:inequality tournillon$ dhs_dados /*.SAV/ dhs_sav_files/ 
Run Code Online (Sandbox Code Playgroud)

我认为这是一项微不足道的任务,但我无法正确理解语法.

非常感谢,

安东尼奥佩德罗.

macos shell

4
推荐指数
1
解决办法
9650
查看次数

如何在R中对变量名称进行子集化?

我有一个如下所示的数据集:

head(data)
                              country      fixef
1             as.factor(country.x)Albania -0.4446439
2             as.factor(country.x)Algeria -0.3400060
3             as.factor(country.x)Andorra -1.0455948
4              as.factor(country.x)Angola  0.7477114
5 as.factor(country.x)Antigua and Barbuda -0.1996655
6           as.factor(country.x)Argentina -0.3404206
Run Code Online (Sandbox Code Playgroud)

我如何删除(在R中)除了国家名称之外的所有其他内容,以便它看起来像这样:

head(data)
                     country      fixef
1                    Albania -0.4446439
2                    Algeria -0.3400060
3                    Andorra -1.0455948
4                    Angola   0.7477114
5        Antigua and Barbuda -0.1996655
6                  Argentina -0.3404206
Run Code Online (Sandbox Code Playgroud)

谢谢,

安东尼奥佩德罗.

r

1
推荐指数
1
解决办法
183
查看次数

重新编码变量

我想重新编码一个变量,以便,例如,我们可以通过以下方式获得变换这个向量:

> a <- c(0,0,0,0,0,1,1,1,1,1) # original 
> b <- c(-5,-4,-3,-2,-1,0,1,2,3,4) # transformed
> cbind(a,b)
  a  b
 [1,] 0 -5
 [2,] 0 -4
 [3,] 0 -3
 [4,] 0 -2
 [5,] 0 -1
 [6,] 1  0
 [7,] 1  1
 [8,] 1  2
 [9,] 1  3
[10,] 1  4
>
Run Code Online (Sandbox Code Playgroud)

这些变量遵循一个顺序,恰好是一个时间顺序.在原始数据集中,我有一个编码为"0"或"1"的变量,例如这里的例子中的"a".它是每年的分类指标.在某些时候,存在从"0"到"1"的转变,就像在这些示例中的行号6中一样.然后我想重新编码原始变量,创建一个新变量,它实际上告诉我从"0"变为"1"之前和之后多少年.因此," - 5"表示转换前五年,"0"表示转换年份,例如,"4"表示转换后四年.有什么建议可以做到最好吗?谢谢!安东尼奥.

r data-manipulation

0
推荐指数
1
解决办法
244
查看次数

标签 统计

r ×3

data-manipulation ×1

gam ×1

macos ×1

mgcv ×1

shell ×1