统计学中的纵向数据分析:M. Davidian的方法

需积分: 9 3 下载量 187 浏览量 更新于2024-07-22 收藏 5.97MB PDF 举报
"Applied Longitudinal Data Analysis - M. Davidian" 本书《Applied Longitudinal Data Analysis》由M. Davidian撰写,旨在概述在分析纵向数据(即,对同一实验单元随时间进行重复测量的数据)时所使用的统计模型和方法。目标读者群体包括已经学习过线性回归和方差分析应用方法的非统计学研究生,以及寻求应用入门课程的硕士水平的统计学学生。 在书中,作者首先介绍了课程的目的和动机,强调了纵向数据分析的重要性。通过实例,展示了这种类型数据在不同领域的应用,如医学研究、社会科学和工程学等。接着,书中详细讲解了适用于纵向数据的统计模型,包括混合效应模型、随机效应模型和动态模型等,这些都是处理随时间变化的变量关系的关键工具。 在回顾矩阵代数这一部分,书中涵盖了矩阵记法、基本操作,为后续的统计建模打下基础。然后,作者讨论了随机向量和多变量正态分布,这是理解随机过程和多元统计分析的基础。多变量正态分布在分析多个相关变量的联合分布时至关重要,而线性回归则用于研究一个或多个自变量与因变量之间的关系。 在介绍如何建模纵向数据时,书中提出基本的统计模型,并分析了纵向数据中的变异来源。这部分内容包括均值结构和协方差结构的探索,以及各种常用的协方差结构模型,如独立性模型、第一阶自相关模型和 Toeplitz 模型等。为了评估模型的适用性,作者还探讨了在假设平稳性条件下的诊断计算方法,并给出了使用 SAS 软件实现这些模型的指导。 此外,书中详细阐述了单变量重复测量的方差分析,这是一种经典的统计方法,用于检验不同处理组之间是否存在显著差异。这部分内容包括基本的情景设定和统计模型,以及如何处理缺失数据和多重比较问题。 总体来说,这本书为读者提供了深入理解和应用纵向数据分析的全面指南,是从事相关研究和数据分析工作的宝贵参考资料。它不仅涵盖了理论基础,还提供了实用的统计软件操作指导,有助于读者将理论知识转化为实际分析能力。

PCX1 = 1.5482 $Shape factor Cfx for longitudinal force PDX1 = 1.1632 $Longitudinal friction Mux at Fznom PDX2 = -0.11154 $Variation of friction Mux with load PDX3 = 0.94173 $Variation of friction Mux with camber squared PEX1 = 0.27 $Longitudinal curvature Efx at Fznom PEX2 = 0.011693 $Variation of curvature Efx with load PEX3 = 0.053303 $Variation of curvature Efx with load squared PEX4 = 0.59223 $Factor in curvature Efx while driving PKX1 = 32.9102 $Longitudinal slip stiffness Kfx/Fz at Fznom PKX2 = 12.7911 $Variation of slip stiffness Kfx/Fz with load PKX3 = -0.11254 $Exponent in slip stiffness Kfx/Fz with load PHX1 = -0.0017527 $Horizontal shift Shx at Fznom PHX2 = 0.00068824 $Variation of shift Shx with load PVX1 = 0.068079 $Vertical shift Svx/Fz at Fznom PVX2 = 0.0023429 $Variation of shift Svx/Fz with load PPX1 = -0.8733 $linear influence of inflation pressure on longitudinal slip stiffness PPX2 = 0.7035 $quadratic influence of inflation pressure on longitudinal slip stiffness PPX3 = -0.0080216 $linear influence of inflation pressure on peak longitudinal friction PPX4 = -0.47776 $quadratic influence of inflation pressure on peak longitudinal friction RBX1 = 18.3369 $Slope factor for combined slip Fx reduction RBX2 = 18.2559 $Variation of slope Fx reduction with kappa RBX3 = 607.8133 $Influence of camber on stiffness for Fx combined RCX1 = 0.96372 $Shape factor for combined slip Fx reduction REX1 = -1.2699 $Curvature factor of combined Fx REX2 = -0.14323 $Curvature factor of combined Fx with load RHX1 = 0.0037359 $Shift factor for combined slip Fx reduction帮我做成表格

2023-06-12 上传

它的具体实现是这样的,再详细解释一下 bool Spline2dConstraint::Add2dBoundary( const std::vector<double>& t_coord, const std::vector<double>& angle, const std::vector<Vec2d>& ref_point, const std::vector<double>& longitudinal_bound, const std::vector<double>& lateral_bound) { if (t_coord.size() != angle.size() || angle.size() != ref_point.size() || ref_point.size() != lateral_bound.size() || lateral_bound.size() != longitudinal_bound.size()) { return false; } Eigen::MatrixXd affine_inequality = Eigen::MatrixXd::Zero(4 * t_coord.size(), total_param_); Eigen::MatrixXd affine_boundary = Eigen::MatrixXd::Zero(4 * t_coord.size(), 1); for (uint32_t i = 0; i < t_coord.size(); ++i) { const double d_lateral = SignDistance(ref_point[i], angle[i]); const double d_longitudinal = SignDistance(ref_point[i], angle[i] - M_PI / 2.0); const uint32_t index = FindIndex(t_coord[i]); const double rel_t = t_coord[i] - t_knots_[index]; const uint32_t index_offset = 2 * index * (spline_order_ + 1); std::vector<double> longi_coef = AffineCoef(angle[i], rel_t); std::vector<double> longitudinal_coef = AffineCoef(angle[i] - M_PI / 2, rel_t); for (uint32_t j = 0; j < 2 * (spline_order_ + 1); ++j) { // upper longi affine_inequality(4 * i, index_offset + j) = longi_coef[j]; // lower longi affine_inequality(4 * i + 1, index_offset + j) = -longi_coef[j]; // upper longitudinal affine_inequality(4 * i + 2, index_offset + j) = longitudinal_coef[j]; // lower longitudinal affine_inequality(4 * i + 3, index_offset + j) = -longitudinal_coef[j]; } affine_boundary(4 * i, 0) = d_lateral - lateral_bound[i]; affine_boundary(4 * i + 1, 0) = -d_lateral - lateral_bound[i]; affine_boundary(4 * i + 2, 0) = d_longitudinal - longitudinal_bound[i]; affine_boundary(4 * i + 3, 0) = -d_longitudinal - longitudinal_bound[i]; } // std::cout << affine_inequality << std::endl; return AddInequalityConstraint(affine_inequality, affine_boundary); }

2023-07-11 上传
2023-07-12 上传