我正在使用GAM来模拟逻辑回归中的时间趋势.然而,我想从中提取拟合样条曲线,将其添加到另一个模型中,不能用于GAM或GAMM.
因此我有两个问题:
随着时间的推移,我怎样才能更顺畅,以便在让模特找到其他结点的同时强迫一个结处于特定位置?
如何从拟合的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的大量文档,但我还不确定.任何建议都非常感谢.
如何将具有特定扩展名的所有文件.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)
我认为这是一项微不足道的任务,但我无法正确理解语法.
非常感谢,
安东尼奥佩德罗.
我有一个如下所示的数据集:
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)
谢谢,
安东尼奥佩德罗.
我想重新编码一个变量,以便,例如,我们可以通过以下方式获得变换这个向量:
> 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"表示转换后四年.有什么建议可以做到最好吗?谢谢!安东尼奥.