MATLAB的答案

计算数据集的置信区间。

조회수:3(최근30일)
표시이전댓글
Rajvi Amle
Rajvi Amle 2021년월23일
  • 0
댓글: dpb 2021년월23일
我有两个变量（SYSTEM.m_BW和SYSTEM.w_L）N = 10和计算出的随机值。2个常微分方程和提及一个循环来为每10次迭代运行（对于i = 1：10）和绘制为相同的曲线图。现在我想计算平均从模拟数据的曲线图的所有行的（X（：，2）+ X（：，1））在回路中所提到的，然后保存X＆T（时间跨度）是一个矩阵，计算和从代码中给出的模拟数据图的置信区间为每x = 1至10行（X（：，1）：+ X（2））。如何计算的平均在所有的情节线，如何计算基于数据点的平均值和标准偏差的置信区间？我希望图形的东西，与从与上边界和下边界示出平均值计算的95％置信区间在附图中给出。
有人能帮我一下吗，这将是非常有帮助和感激的。以下是给定的代码，我已经工作到目前为止:
函数toy_example
清除所有;
关闭所有;
clc;
%体重和器官重量百分比[Min - Max](10只小鼠随机数字)
闵％= A
% Max = b
% = + (b) * r (n, 1);
n = 10;
系统。m_BW = 24 +(28-24)。*兰德(10,1);％的身体重量[g]
系统。w_L = 0.0417 +(0.0681 - -0.0417) *兰德(10,1);％肝[克/克]
系统。w_S = 0.019 +(0.0051 - -0.0019) *兰德(10,1);%脾(g / g)
系统。w_K = 0.015 +(0.0184 - -0.015) *兰德(10,1);%肾脏(g / g)
系统。w_Lu = 0.0065 +(0.0081 - -0.0065) *兰德(10,1);%肺(g / g)
。SYSTEM.w_BR = 0.0139+（0.0191-0.0139）*兰特（10,1）;％脑[克/克]
。SYSTEM.w_Blood = 0.049+（0.049-0.049）*兰特（10,1）;%的血液(g / g)
系统。w_Plasma = 0.029 +(0.029 - -0.029) *兰德(10,1);％血浆[克/克]
EXP.t = [0 6 12 18 24];％H（时间血液动力学）
药物。w_Au_iv = [15E-03 183E-03 693E-03 1058E-03];%(毫克/克)
药物。m_Au_iv = DRUG.w_Au_iv。* SYSTEM.m_BW (1);%(毫克)
EXP.Blood_15 =(((16.32 3.8 - 2.87 2.33 - 2.06) / 100)。* * SYSTEM.w_Plasma SYSTEM.m_BW。* DRUG.w_Au_iv(1)。* SYSTEM.m_BW)。/ (SYSTEM.m_BW。* SYSTEM.w_Plasma);
EXP.Blood_183 =(((31.3 1.48 - 0.98 0.84 - 0.8) / 100)。* * SYSTEM.w_Plasma SYSTEM.m_BW。* DRUG.w_Au_iv(2)。* SYSTEM.m_BW)。/ (SYSTEM.m_BW。* SYSTEM.w_Plasma);
EXP.Blood_1058 =(((57.82 0.42 - 0.35 0.35 - 0.35) / 100)。* * SYSTEM.w_Plasma SYSTEM.m_BW。* DRUG.w_Au_iv(4)。* SYSTEM.m_BW)。/ (SYSTEM.m_BW。* SYSTEM.w_Plasma);
数字
持有在
为i = 1:10
tspan = [0:24];
垫= [SYSTEM.m_BW (i) SYSTEM.w_L(我)';
x0 = 0 (2, 1);
抽搐;
[t,x]=ode15s(@ode_toy, tspan, x0， []， MAT);
TOC;
情节（T，（X（：，2）+ X（：，1）），“-”。，'颜色'，“蓝”，“线宽”,2)
结束
结束
功能dxdt = ode_toy (t x系统)
分裂状态向量
m_BW =系统(1);
w_L =系统(2);
m_Au_A = x (1);
m_Au_V = x (2);
% % ODE体系
%剂量
dm_Au_A_dt = m_BW - w_L * m_Au_A;
dx (1) = dm_Au_A_dt;
dm_Au_V_dt= w_L - m_BW*m_Au_V;
DX（2）= dm_Au_V_dt;
dxdt = dx (:);
结束
댓글수:3
표시安康이전댓글수:2
dpb
dpb 2021년월23일
我时间很紧，所以这里输入的是未经测试的“air code”…希望这足以让你明白……
= t细胞(10);y = t;%将输出单元格数组预分配为容纳列向量的行
数字
持有在
tspan = [0:24];%不变，不需要在循环中
为i = 1:10
垫= [SYSTEM.m_BW (i) SYSTEM.w_L(我)';
x0 = 0 (2, 1);
[t{i}，x]=ode15s(@ode_toy, tspan, x0，[]， MAT);％分配输出到单元阵列
y {} = ((:, 2) + x (: 1);％计算，节省另外，Y
情节(y t{我},{},“-”。，'颜色'，“蓝”，“线宽”,2)
结束
敌我识别 返回的 t Vector对于每个迭代都是相同的，然后你可以将它转换成一个单独的向量，所有的都是多余的 y 转换为普通的2D数组。
然后均值sd是
mny =意味着(cell2mat (y), 2);将单元格数组转换为矩阵，平均按行
sdmny = mny /大小(y, 2);％平均值（Y）的标准偏差是平均值（Y）/ NOBS

댓글을달려면로그인하십시오。

이질문에답변하려면로그인하십시오。

답변(0개)

범주

태그

社区寻宝

在MATLAB中心找到宝藏，并发现社区如何可以帮助你!

开始狩猎！

翻译的