分层回归模型MATLAB代码

clc;
clear all;

% 读取数据
data = readmatrix("副本赛题数据.xlsx", 'Sheet', 1);
x1 = data(:, 3);  % 学校类型
x2 = data(:, 9);  % 中考成绩
y = data(:, 15);  % 高考成绩

% 创建分层变量（假设学校类型的取值范围为0-5）
group = categorical(x1);

% 构建线性混合效应模型
formula = 'y ~ x1 + x2 + (1 | group)';  % 设置模型公式，考虑学校类型、中考成绩，并包含随机效应
model = fitlme(table(x1, x2, y, group), formula);  % 进行分层回归拟合

% 输出模型摘要
disp(model);

% 获取固定效应的估计值
fixedEffects = model.Coefficients;

% 输出固定效应估计值
disp('固定效应估计值:');
disp(fixedEffects);