请给我Benjamini–Hochberg procedure方法校正p值的matlab代码

以下是使用Matlab实现Benjamini-Hochberg过程来校正p值的示例代码：

function [adjusted_p] = benjamini_hochberg(p_values, alpha)
    % 对p值进行Benjamini-Hochberg校正
    
    % 获取p值的个数
    m = length(p_values);
    
    % 根据p值排序，并获取排序后的索引
    [~, sorted_idx] = sort(p_values);
    
    % 计算校正因子
    q_values = (1:m) * alpha / m;
    
    % 初始化调整后的p值向量
    adjusted_p = NaN(1, m);
    
    % 逐个检查每个p值
    for i = 1:m
        % 计算校正后的p值
        adjusted_p(sorted_idx(i)) = min(1, p_values(sorted_idx(i)) * m / i);
        
        % 如果当前校正后的p值大于之前的校正后的p值，则使用之前的校正后的p值
        if i > 1
            adjusted_p(sorted_idx(i)) = min(adjusted_p(sorted_idx(i)), adjusted_p(sorted_idx(i-1)));
        end
    end
    
end

使用示例：

% 示例数据
p_values = [0.01, 0.02, 0.03, 0.04, 0.05];
alpha = 0.05;

% 进行校正
adjusted_p = benjamini_hochberg(p_values, alpha);

% 输出校正后的p值
disp(adjusted_p);

请注意，这只是一个简单的示例代码，你可以根据自己的需求进行修改和优化。