开发者社区

文档建议反馈控制台

最新优惠活动

文章/答案/技术大牛

发布

在混合效应模型中，如何将时间=0时的结果基线值作为固定效应添加？

在混合效应模型中，将时间=0时的结果基线值作为固定效应添加的方法是通过引入一个虚拟变量来表示时间。虚拟变量通常是一个二进制变量，用来表示时间是否等于0。当时间等于0时，虚拟变量取值为1，否则取值为0。

在混合效应模型中，固定效应可以通过以下步骤添加：

创建一个虚拟变量：根据时间是否等于0，创建一个二进制虚拟变量。例如，可以使用dummy变量表示时间是否等于0，dummy=1表示时间等于0，dummy=0表示时间不等于0。
添加固定效应：将虚拟变量作为自变量添加到混合效应模型中。固定效应表示时间=0时的结果基线值。通过将虚拟变量与其他自变量一起作为固定效应添加到模型中，可以控制时间=0时的结果。
解释固定效应：根据模型的系数估计值，解释固定效应的影响。系数估计值表示时间=0时的结果相对于其他时间点的差异。

需要注意的是，以上方法是一种常见的处理时间=0时结果基线值的方法，具体的实施方式可能因具体的混合效应模型而有所不同。在实际应用中，可以根据具体情况进行调整和优化。

腾讯云相关产品和产品介绍链接地址：

腾讯云混合云解决方案：https://cloud.tencent.com/solution/hybrid-cloud
腾讯云云服务器（CVM）：https://cloud.tencent.com/product/cvm
腾讯云数据库（TencentDB）：https://cloud.tencent.com/product/cdb
腾讯云人工智能（AI）：https://cloud.tencent.com/product/ai
腾讯云物联网（IoT）：https://cloud.tencent.com/product/iot
腾讯云移动开发（移动推送、移动分析等）：https://cloud.tencent.com/product/mobile
腾讯云对象存储（COS）：https://cloud.tencent.com/product/cos
腾讯云区块链（BCS）：https://cloud.tencent.com/product/bcs
腾讯云元宇宙（Tencent XR）：https://cloud.tencent.com/product/xr

相关搜索:如何使用UseEffect()钩子将数据从后台加载到材料表？如何使用GmailR包发送多个附件和图片(邮件中)正文？在Vue组件中导入外部类时出错对Typescript中可为空的类型进行Nullcheck 如何根据特定的索引将多个字符数组移到数组的后面？"PassthroughSubject“似乎是线程不安全的，是错误还是限制？将字典列表中的字典值替换为另一个字典列表中的值 Dockerized Python Project -管理依赖项的最佳方式: Dockerfile或requirements.txt 如何根据当前文件夹名称对旧目录进行排序和删除即使所有字段都有值，django is_valid也会为模型表单返回false

相关搜索:

页面内容是否对你有帮助？

有帮助

没帮助

相关·内容

广义估计方程和混合线性模型在R和python中的实现

，可以得到回归系数及其方差的一致性估计混合线性模型(mixed linear model，MLM)：它是一类对误差进行精细分解成对固定效应和随机效应等误差的广义线性模型的方法，相比广义线性模型而言，它能处理纵向数据...P*P维作业相关矩阵（自变量X），用以表示因变量的各次重复测量值（自变量）之间的相关性大小求参数$\beta$的估计值及其协方差矩阵混合线性模型(mixed linear model，MLM)：构建包含固定因子和随机因子的线性混合模型...区分混合线性模型中的随机效应和固定效应是一个重要的概念。固定效应是具有特定水平的变量，而随机效应捕捉了由于分组或聚类引起的变异性。比如下方正在探究尿蛋白对来自不同患者的GFR的影响。...固定效应：具有特定的水平或值需要进行研究的主要变量，如尿蛋白等随机效应：患者分层结构：尿蛋白嵌套在患者内模型方程：GFR = 尿蛋白 + 患者 + 误差解释：解释固定效应，以了解尿蛋白的变化如何与GFR...- 实例操作及结果解读（R、Python、SPSS实现）混合线性模型介绍--Wiki广义估计方程中工作相关矩阵的选择及R语言代码在Rstudio 中使用pythonAn Introduction to

2480 0

R语言nlme、nlmer、lme4用（非）线性混合模型non-linear mixed model分析藻类数据实例|附代码数据

简单的说，混合模型中把研究者感兴趣的自变量对因变量的影响称为固定效应，把其他控制的情景变量称为随机效应。由于模型中包括固定和随机效应，故称为混合线性模型。...结果不理想 fixef(nlmerfit2) range(predict(nlmerfit2)) 我不能确定，在nlmer中是否有更简单的方法来做固定效果。..., run.opts) 结果不错固定效应模型现在尝试用固定效应分组，使用上面构建的虚拟变量（也可以使用if语句，或者用R[Group[i]]的for循环中的R值向量，或者（最佳选择..." 当然，在参数相同的情况下，也可以工作。随机效应现在添加随机效应。...我们可以尝试在xmid和scale参数中加入随机效应。在组间或作为X的函数的方差（无论是残差还是个体间的方差）中可能有额外的模式。

8240 0

R语言广义线性混合模型GLMMs在生态学中应用可视化2实例合集|附数据代码

展示了如何模拟固定效应和混合效应模型的数据，并比较了这两种模型的拟合结果。您的目标是展示固定效应模型和混合效应模型之间的区别，并通过绘图来可视化这些差异。...# 结论： # 在这个模拟例子中，固定效应和混合效应模型的拟合结果没有明显的视觉差异， # 这可能是因为随机效应的影响较小。...代码从mod3_lmer（只包含NAP作为固定效应的模型）中模拟新的观测值，然后拟合零模型和替代模型，并计算似然比检验统计量。最后，基于模拟的统计量计算p值，以评估固定效应Exposure是否显著。...拟合GLMM模型： glmer 函数用于拟合广义线性混合效应模型，这里以物种丰富度（Richness）作为响应变量，NAP和Exposure作为固定效应，同时考虑到Beach作为随机效应。...这两项分析从我们的原始模型 1 开始，该模型已有 10 年的研究时间。添加更多组我们可以像为_x _添加额外值一样为_g _添加额外级别。

4661 0

Nature neuroscience：结构束的改变预示着淀粉样蛋白阳性老年人的下游tau蛋白累积

采用基于最大似然估计的逐步假设驱动的线性混合效应(LME)建模方法进行纵向分析，其中包括固定效应的兴趣预测器、每个受试者的随机截距和随机斜率(基线和随访之间的年数)。...在第四步中，我们通过在模型中添加三向交互作用（基线时的HCB扩散率×淀粉样状态×时间），检查了由HCB扩散率预测的tau变化对于淀粉样蛋白阳性和淀粉样蛋白阴性个体是否有所不同。...然而，神经退化相关连通性缺失是否预示着HCB解剖端区域内的累积增加。为此，线性混合效应模型包括PCC中tau蛋白的年积累作为结果的衡量指标，白质扩散率作为基线的预测指标。...临界值为1.28，采用高斯混合建模方法确定(补充图6)。将这些模型与PCC tau结合作为连续测量，得到相似的结果。记忆表现评分是基于整个HABS的因素评分。...研究者以基线HCB扩散率为预测因子，以记忆性能的年变化为结果，进行线性混合效应模型。

7083 0

R语言线性混合效应模型（固定效应&随机效应）和交互可视化3案例|附代码数据

在本文中，我们将用R语言对数据进行线性混合效应模型的拟合，然后可视化你的结果线性混合效应模型是在有随机效应时使用的，随机效应发生在对随机抽样的单位进行多次测量时。...这两个来源中的哪一个对应于"（截距）"，哪一个对应于 "残差"？同时检查固定效应结果的输出。模型公式中唯一的固定效应是所有长度测量的平均值。它被称为"（截距）"，但不要与随机效应的截距相混淆。...从保存的lmer对象中提取参数估计值。检查固定效应的结果。给出的系数与使用lm分析的分类变量的解释相同。检查随机效应的输出。我们的混合效应模型中再次出现了两个随机误差的来源。它们是什么？...在每个实验和持续时间水平的组合中没有很多数据点，所以按组画条形图可能比按组画箱形图更好。添加线段来连接成对的点。拟合一个线性混合效应模型使用的是什么类型的实验设计？...*这将决定对数据的线性混合模型的拟合。在没有实验和持续时间之间的交互作用的情况下，对数据进行线性混合模型拟合。使用酚类物质的对数作为因变量，因为对数转换改善了数据与线性模型假设的拟合。

1.5K0 0

R语言线性混合效应模型（固定效应&随机效应）和交互可视化3案例|附代码数据

在本文中，我们将用R语言对数据进行线性混合效应模型的拟合，然后可视化你的结果线性混合效应模型是在有随机效应时使用的，随机效应发生在对随机抽样的单位进行多次测量时。...这两个来源中的哪一个对应于"（截距）"，哪一个对应于 "残差"？同时检查固定效应结果的输出。模型公式中唯一的固定效应是所有长度测量的平均值。它被称为"（截距）"，但不要与随机效应的截距相混淆。...从保存的lmer对象中提取参数估计值。检查固定效应的结果。给出的系数与使用lm分析的分类变量的解释相同。检查随机效应的输出。我们的混合效应模型中再次出现了两个随机误差的来源。它们是什么？...在每个实验和持续时间水平的组合中没有很多数据点，所以按组画条形图可能比按组画箱形图更好。添加线段来连接成对的点。拟合一个线性混合效应模型使用的是什么类型的实验设计？...*这将决定对数据的线性混合模型的拟合。在没有实验和持续时间之间的交互作用的情况下，对数据进行线性混合模型拟合。使用酚类物质的对数作为因变量，因为对数转换改善了数据与线性模型假设的拟合。

1.1K3 0

R语言线性混合效应模型（固定效应&随机效应）和交互可视化3案例

p=23050 在本文中，我们将用R语言对数据进行线性混合效应模型的拟合，然后可视化你的结果。线性混合效应模型是在有随机效应时使用的，随机效应发生在对随机抽样的单位进行多次测量时。...这两个来源中的哪一个对应于"（截距）"，哪一个对应于 "残差"？同时检查固定效应结果的输出。模型公式中唯一的固定效应是所有长度测量的平均值。它被称为"（截距）"，但不要与随机效应的截距相混淆。...从保存的lmer对象中提取参数估计值。检查固定效应的结果。给出的系数与使用lm分析的分类变量的解释相同。检查随机效应的输出。我们的混合效应模型中再次出现了两个随机误差的来源。它们是什么？...在每个实验和持续时间水平的组合中没有很多数据点，所以按组画条形图可能比按组画箱形图更好。添加线段来连接成对的点。拟合一个线性混合效应模型使用的是什么类型的实验设计？...*这将决定对数据的线性混合模型的拟合。在没有实验和持续时间之间的交互作用的情况下，对数据进行线性混合模型拟合。使用酚类物质的对数作为因变量，因为对数转换改善了数据与线性模型假设的拟合。

8.5K6 1

用SPSS估计HLM多层（层次）线性模型模型|附代码数据

接下来，单击Statistics以选择其他菜单以选择在输出中报告哪些结果。选择参数估计值报告固定效应的估计值。单击继续，然后单击确定。部分结果如下：这些结果对应于R＆B中的表4.2。...下一步是估计一种平均数- 结果模型。平均数之结果变项的回归模型在估计空模型之后，R＆B开发了一种“平均数结果变项的回归”模型，其中将学校级变量meanses添加到截距模型中。...在下一个菜单中，指定依赖变量和独立变量。因变量将是mathach，单个协变量将是均值。该meanses变量输入作为固定效应，所以点击固定按钮拉起固定效应菜单。...在“ 固定效应”菜单中，将grp_ses变量置于“ 模型”框中，并确保选中“ 包括截距”。单击继续，然后单击随机。在“ 随机效应”菜单中，将分组变量id放在“ 组合”框中。...可以放宽该假设，使得协方差是从数据估计的自由参数。为协方差类型指定Unstructured。单击继续。然后单击“ 统计”以指定输出中显示的内容。检查参数估计值以获得固定效应的结果。

2.2K1 0

Molecular Psychiatry：抑郁症患者经ECT治疗后齿状回增大的7T磁共振成像研究

统计分析对于每个亚区的时间和组别的交互作用分别通过R软件包的线性混合模型进行重复测量，其中时间、组别、年龄和性别作为固定因素，半球（左右半球建模为为个体被试的斜率）和个体被试作为随机因素（建模为截距）...将对时间和组别交互作用有显著影响的线性混合模型进一步分为患者模型和对照组模型，分别测试哪一组驱动效果。...如果患者和/或对照组在后一项分析中显示出明显的时间效应，则分别对左侧和右侧齿状回建立线性混合模型，以观察哪个亚区驱动了这种效应。...并且在线性回归中，将抑郁评分的降低作为因变量，将线性混合模型重要亚区基线时的体积作为预测变量，年龄和性别则作为协变量。并进行事后t检验。...虚线表示了B相关3D海马的位置。线性混合模型表明在齿状回出现明显的时间组别交互效应(t = −2.57, p = 0.0138)，其它亚区并未出现。

8132 0

R语言建立和可视化混合效应模型mixed effect model|附代码数据

本文的重点是如何建立和_可视化_ 混合效应模型的结果设置本文使用数据集，用于探索草食动物种群对珊瑚覆盖的影响。...urchinden，c.fishmass，c.maxD）和调查地点作为随机效应（地点）。。注意：由于食草动物种群的测量规模存在差异，因此我们使用标准化的值，否则模型将无法收敛。...title="草食动物对珊瑚覆盖的影响")模型结果表输出：创建模型摘要输出表。这将提供预测变量，包括其估计值，置信区间，估计值的p值以及随机效应信息。...注意：数据已标准化以便在模型中使用，因此我们绘制的是标准化数据值，而不是原始数据步骤1：将效应大小估算值保存到data.frame中# 使用函数。term=固定效应，mod=你的模型。...＃基本步骤：＃1创建空图 #2 从数据中添加geom_points（） #3 为模型估计添加geom_point。

1.5K0 0

R语言建立和可视化混合效应模型mixed effect model|附代码数据

p=20631 最近我们被客户要求撰写关于混合效应模型的研究报告，包括一些图形和统计输出。我们已经学习了如何处理混合效应模型。本文的重点是如何建立和_可视化_ 混合效应模型的结果。...<- read_csv("mixede.csv") 创建一个基本的混合效应模型：该模型以珊瑚覆盖层为因变量（elkhorn_LAI），草食动物种群和深度为固定效应（c。...urchinden，c.fishmass，c.maxD）和调查地点作为随机效应（地点）。。注意：由于食草动物种群的测量规模存在差异，因此我们使用标准化的值，否则模型将无法收敛。...title="草食动物对珊瑚覆盖的影响") 模型结果表输出：创建模型摘要输出表。这将提供预测变量，包括其估计值，置信区间，估计值的p值以及随机效应信息。...注意：数据已标准化以便在模型中使用，因此我们绘制的是标准化数据值，而不是原始数据步骤1：将效应大小估算值保存到data.frame中 # 使用函数。term=固定效应，mod=你的模型。

1.5K2 0

R语言建立和可视化混合效应模型mixed effect model

urchinden，c.fishmass，c.maxD）和调查地点作为随机效应（地点）。。注意：由于食草动物种群的测量规模存在差异，因此我们使用标准化的值，否则模型将无法收敛。...title="草食动物对珊瑚覆盖的影响") 模型结果表输出：创建模型摘要输出表。这将提供预测变量，包括其估计值，置信区间，估计值的p值以及随机效应信息。...注意：数据已标准化以便在模型中使用，因此我们绘制的是标准化数据值，而不是原始数据步骤1：将效应大小估算值保存到data.frame中 # 使用函数。term=固定效应，mod=你的模型。...＃基本步骤：＃1创建空图 #2 从数据中添加geom_points（） #3 为模型估计添加geom_point。...R语言估计多元标记的潜过程混合效应模型（lcmm）分析心理测试的认知过程 R语言因子实验设计nlme拟合非线性混合模型分析有机农业施氮水平 R语言非线性混合效应 NLME模型(固定效应&随机效应)

2392 0

SAS，Stata，HLM，R，SPSS和Mplus分层线性模型HLM分析学生受欢迎程度数据|附代码数据

仅截距模型（无条件模型）无条件混合模型规范类似于单因素方差分析，其总体均值和类效应。但是，我们将其视为随机效应（均值为零的正态分布变量），而不是像方差分析中那样的固定因子效应。...混合模型看起来像是基于带有协变量Extrav的类的ANCOVA，但请记住，我们仍然认为这是随机效应，而不是固定效应。因此，估算值与ANCOVA程序所得出的估算值不同。 ...对Sex的固定估计意味着，在Extrav不变的情况下，女学生（Sex = 1）的普遍得分比男学生（基线组，Sex = 0）高1.244。...Mplus结果也显示出比以前的模型更大的差异。此模型的ICC为：同样，当我们在模型中添加另一个学生级别的效果（包括随机斜率）时，ICC略有增加。 ...这与混合模型中的交互项有关，即通过外向的教师体验和按性别的教师体验。 SAS结果在固定效果表中，有两个交互作用项，其中一个（）远不重要，p值> 0.5。

1.4K1 0

SAS，Stata，HLM，R，SPSS和Mplus分层线性模型HLM分析学生受欢迎程度数据

仅截距模型（无条件模型）无条件混合模型规范类似于单因素方差分析，其总体均值和类效应。但是，我们将其视为随机效应（均值为零的正态分布变量），而不是像方差分析中那样的固定因子效应。...混合模型看起来像是基于带有协变量Extrav的类的ANCOVA，但请记住，我们仍然认为这是随机效应，而不是固定效应。因此，估算值与ANCOVA程序所得出的估算值不同。...对Sex的固定估计意味着，在Extrav不变的情况下，女学生（Sex = 1）的普遍得分比男学生（基线组，Sex = 0）高1.244。...Mplus结果也显示出比以前的模型更大的差异。此模型的ICC为：同样，当我们在模型中添加另一个学生级别的效果（包括随机斜率）时，ICC略有增加。...这与混合模型中的交互项有关，即通过外向的教师体验和按性别的教师体验。 SAS结果在固定效果表中，有两个交互作用项，其中一个（）远不重要，p值> 0.5。

2.5K1 0

SAS，Stata，HLM，R，SPSS和Mplus分层线性模型HLM分析学生受欢迎程度数据|附代码数据

仅截距模型（无条件模型）无条件混合模型规范类似于单因素方差分析，其总体均值和类效应。但是，我们将其视为随机效应（均值为零的正态分布变量），而不是像方差分析中那样的固定因子效应。...混合模型看起来像是基于带有协变量Extrav的类的ANCOVA，但请记住，我们仍然认为这是随机效应，而不是固定效应。因此，估算值与ANCOVA程序所得出的估算值不同。 ...对Sex的固定估计意味着，在Extrav不变的情况下，女学生（Sex = 1）的普遍得分比男学生（基线组，Sex = 0）高1.244。...Mplus结果也显示出比以前的模型更大的差异。此模型的ICC为：同样，当我们在模型中添加另一个学生级别的效果（包括随机斜率）时，ICC略有增加。 ...这与混合模型中的交互项有关，即通过外向的教师体验和按性别的教师体验。 SAS结果在固定效果表中，有两个交互作用项，其中一个（）远不重要，p值> 0.5。

1.7K2 0

使用SAS，Stata，HLM，R，SPSS和Mplus的分层线性模型HLM

仅截距模型（无条件模型）无条件混合模型规范类似于单因素方差分析，其总体均值和类效应。但是，我们将其视为随机效应（均值为零的正态分布变量），而不是像方差分析中那样的固定因子效应。...混合模型看起来像是基于带有协变量Extrav的类的ANCOVA，但请记住，我们仍然认为这是随机效应，而不是固定效应。因此，估算值与ANCOVA程序所得出的估算值不同。 ...对Sex的固定估计意味着，在Extrav不变的情况下，女学生（Sex = 1）的普遍得分比男学生（基线组，Sex = 0）高1.244。...Mplus结果也显示出比以前的模型更大的差异。此模型的ICC为：同样，当我们在模型中添加另一个学生级别的效果（包括随机斜率）时，ICC略有增加。 ...这与混合模型中的交互项有关，即通过外向的教师体验和按性别的教师体验。 SAS结果在固定效果表中，有两个交互作用项，其中一个（）远不重要，p值> 0.5。

2.9K2 0

Nature子刊：生命的最初几个小时内可快速学习音位辨别—基于fNIRS研究

Δ[HbO]的平均振幅和峰值潜伏期采用线性混合效应回归建模，通过R包lme4进行单次试验的振幅或潜伏期分析，建立为三个中心化和编码的固定效应（刺激类型2*被试组3*测试阶段3）及其相互作用的函数模型。...然后将r值转换为Fisher的z分数，以进行进一步的统计分析，然后拟合了用于振幅和延迟分析的线性混合效应回归的相同简化模型，即包括每个受试者和每一对通道的最大随机效应结构。...回归模型忽略了刺激类型及其相互作用的固定和随机效应，因为测试前阶段没有包括这两种类型的刺激，但作为结构考虑，在主动对照vs实验的对比中保留了刺激类型的交互作用。...结果氧合血红蛋白浓度振幅分析对刺激开始后6-16秒的每次试验的平均氧合血红蛋白浓度[HbO]振幅进行线性混合效应回归分析确定刺激类型，组间（主动控制和实验和第二阶段）和测量时间 (T1和T2）存在显著的三者交互作用...我们首先在每次实验中确定[HbO]max的潜伏期，然后将相同的线性混合效应回归模型拟合到我们之前拟合到平均振幅数据的潜伏期数据中，发现在刺激类型间，组间（主动控制vs实验）和测试时间对比（T0 vs平均值

4582 0

R语言广义线性混合模型（GLMM）bootstrap预测置信区间可视化

library(lme4) # 加载lme4包，用于线性混合效应模型的分析 # 第一个案例：简单的线性混合效应模型，从10个组中模拟100个数据点，具有一个连续的固定效应变量...x <- runif(100, 0, 10) # 生成100个在0到10之间的均匀分布随机数，作为固定效应变量x # 固定效应系数 fixed <- c(1, 0.5) # 设定固定效应系数...在newdat数据框中添加预测值、预测区间的下限和上限、置信区间的下限和上限 newdat <- data.frame( newdat, plo = newdat$y -...：使用模型矩阵和固定效应系数来计算新数据点的预测值。...此外，bootMer函数可能需要较长时间来执行，特别是当模型复杂或自助法抽样次数较多时。在上述代码中，模拟数据的生成和模型的拟合都是基于线性混合效应模型（LMM）的。

1931 0

R语言估计多元标记的潜过程混合效应模型（lcmm）分析心理测试的认知过程

使用线性混合模型根据时间对定义为潜过程的感兴趣量进行建模：其中： X(t) 和 Z(t) 是协变量的向量（Z(t) 包含在 X(t) 中； β是固定效应（即总体平均效应）； ui 是随机效应（即个体效应...这里第一个随机效应 ui的方差设置为 1，平均截距（在 β 中）设置为 0。...为了进一步研究性别的影响，包括对共同因素的平均效应和对每个标志的差异效应（对比）（在这个例子中不与时间相互作用）。...为此，可以从估计向量（此处为第 21 个参数）中识别参数的位置： best 并且可以根据这些估计值和新固定的参数重新拟合模型： # 样条曲线 mult(B=mp$best) 有了这个约束，模型就可以正确收敛...summary 从估计结果来看，基础认知随着时间的推移有一个二次方的轨迹，基线时年龄较大的受试者的认知水平系统地较低。根据性别没有差异。

1764 0

R语言混合效应逻辑回归（mixed effects logistic）模型分析肺癌数据|附代码数据

p=22302在混合效应逻辑回归用于建立二元结果变量的模型，其中，当数据被分组或同时存在固定和随机效应时，结果的对数几率被建模为预测变量的线性组合（点击文末“阅读原文”获取完整代码数据*******...#估计模型并将结果存储在m中# 输出结果，固定效果之间不相关print(m, corr = FALSE)第一部分告诉我们，估计值是基于自适应高斯-赫米特的似然性近似。...在我们的案例中，这包括病人（8,525）和医生（407）的总数。最后一节是固定效应估计值的表格。这些估计值代表回归系数。这些是未标准化的，而且是在对数尺度上。估计值后面是它们的标准误差（SE）。...tmp <- sample(hdp, "DID", reps = 100)接下来，我们在重新取样的数据上重新拟合模型。首先，我们存储原始模型的估计值，我们将用它作为自助模型的起始值。...在普通逻辑回归中，你可以保持所有预测因子不变，只改变你感兴趣的预测因子。然而，在混合效应逻辑模型中，随机效应也对结果产生影响。

7960 0

点击加载更多

扫码

添加站长进交流群

领取专属 10元无门槛券

手把手带您无忧上云

扫码加入开发者社群

相关资讯

热门标签

活动推荐

运营活动

活动名称

广告关闭