我在Matlab中的曲线拟合工具箱中使用由2个求和高斯组成的拟合拟合了一些数据 . Matlab能够非常好地为我绘制95%的上下置信区间 .

现在,我有一些新数据最适合由3个高斯组成的曲线,但在这种情况下工具箱不会绘制95%CI . 它沿曲线的所有点返回NaN,而不是95%CI的数值 .

以下是一些代码:

X = [0 1 2 3 4 5 6 7 8 9 10]';

Y = [0.009 0.014 0.031 0.044 0.062 0.065 0.052 0.037 0.037 0.027 0.026]';

STDV = [0.0028 0.0013 0.0047 0.00356 0.0077 0.0141 0.0072 0.0055 0.0055 0.0040 0.0055]';

[a1,b1,c1] = deal(0.036, 4.34, 2.02);

[a2,b2,c2] = deal(0.034, 6.78, 6.02);

fo1 = fitoptions('method','NonlinearLeastSquares','Lower',[-Inf -Inf 0]);

ok1 = isfinite(X) & isfinite(Y);

if ~all( ok1 )

warning( 'GenerateMFile:IgnoringNansAndInfs',...

'Ignoring NaNs and Infs in data.' );

end

st1 = [a1 b1 c1 a2 b2 c2];

set(fo1,'Startpoint',st1);

ft1 = fittype('gauss2');

% Fit this model using new data

cf1 = fit(X(ok1),Y(ok1),ft1,fo1);

% Plot this fit and the calculated dose points

h1 = plot(cf1,'fit',0.95);

set(h1(1),'Color',[0 0 0],...

'LineStyle','-', 'LineWidth',1,...

'Marker','none', 'MarkerSize',6);

hold on;

errorbar(X,Y,STDV,'*k')

%Break dose range into many intervals

Dose_range = 0:0.01:10;

%Plot confidence interval fits

load('ratio_lower.mat');

load('ratio_upper.mat');

plot(Dose_range,ratio_lower,'k:');

plot(Dose_range,ratio_upper,'k:');

'ratio_lower'和'ratio_upper'是从cftool中的分析GUI生成的CI . 它适用于2个求和的高斯,但不适用于3个求和的高斯 .

区间求和Java,matlab中95%置信区间的求和高斯分布相关推荐

  1. Java区间拆分子集求和,对列表中的数字子集求和

    Is there a way in Kotlin for doing the sum() operation on a filtered list of numbers, without actual ...

  2. matlab 把区间等分分,matlab中怎么把坐标轴刻度等分

    1.matlab怎么设置横坐标刻度 选中坐标轴,单击鼠标右键,选择properties选项,在定义的坐标轴的数值,可以随便调整坐标轴范围. 坐标轴设置 范围设置: a. axis([xmin xmax ...

  3. 如何理解95%置信区间

    如何正确理解95%置信区间 科技论文里经常会出现[95%CI]的评价,这个评价到底有什么意义,他和68-95-99.7法则的关系是什么,可能很多人没有清楚的理解,包括之前写论文评价95%CI的自己. ...

  4. java8 bigdecimal_在java 8中并行使用BigDecimal进行分组和求和

    我有一个产品列表,其中包含其中一个属性.并且列表可以包含公共产品名称,其他属性不同.所以我想按产品和使用分组和求和在java 8中共享公用名的产品数量的总和对列表进行分组. Example: [ { ...

  5. matlab中调用java代码_Matlab中调用第三方Java代码

    在Java中采用Matlab JA Builder可以实现调用m文件,采用这样的方式,可在Matlab的M文件中,直接调用Java类.这种方式可以表示为Java--> Matlab( m, Ja ...

  6. R语言生成数据及其95%置信区间数据(或者其它区间)、使用geom_ribbon函数为ggplot2图像手动添加置信区间的阴影区域using geom_ribbon to create shaded

    R语言生成数据及其95%置信区间数据(或者其它区间).使用geom_ribbon函数为ggplot2图像手动添加置信区间的阴影区域(using geom_ribbon to create shaded ...

  7. matlab矩阵对某一列求和,将矩阵中的每一列与另一列中的对应行相乘,然后在Matlab中求和...

    可以说我有矩阵A = [1 2; 3 4], B = [4 3; 2 1]. 我想将矩阵A ([1; 3], [2; 4])中的每一列乘以矩阵B ([4 3], [2 1])中的相应行,并对得出的矩阵 ...

  8. Java实验9 T3.对二进制数据文件中的所有数据求和

    题目要求 (对二进制数据文件中的所有数据求和)假设已经使用DataOutputStream中的writeInt(int)创建一个名为Exercise 19_3.dat的二进制数据文件,文件包含数目不确 ...

  9. matlab怎么做到使用java,如何在Matlab中使用java.nio?

    我的目标是检查在网络上的文件夹中是否找到具有特定名称(名称的一部分)的文件,同时还要考虑其下的所有文件夹.为此,我需要一种有效获取给定文件夹内和下的所有文件和文件夹的列表的方法.我的递归函数在本地驱动 ...

最新文章

  1. wamp环境下安装imagick扩展
  2. 磁盘镜像工具Guymager
  3. 只不过是R.java文件的特性-----出错信息:R.java was modified manually! Reverting to generated version!...
  4. 滴滴笔试准备 项目分配利益最大化
  5. linux系统怎样写单片机程序,单片机知识是Linux驱动开发的基础之一以及如何学单片机...
  6. MFC开发-树形菜单
  7. 电商后台管理项目知识点总结
  8. 浅谈低成本视频动作捕捉技术
  9. 软件开发实训(720科技)――第五课:前端css规范
  10. 成功解决data too long for column 'name' at row 1
  11. DBCC CHECKDB用法 手工修复数据库
  12. 极客日报第 39 期:小米 11 不送充电器;阿里达摩院 2021 十大科技趋势发布!
  13. Bellman-Ford算法图解及手算过程详解 —— C++代码实现
  14. Android App架构设计
  15. csdn怎么了???
  16. 【CH9434】使用CH9434移植到32单片机扩展4串口
  17. Eclipse Other Projects
  18. 扫地机器人杂牌的怎么样_扫地机器人贵的和便宜的之间有何区别?
  19. docker常用命令-docker start
  20. http://blog.csdn.net/guolin_blog/article/details/17482095/

热门文章

  1. 2022年中式烹调师(中级)最新解析及中式烹调师(中级)找解析
  2. org.apache.zookeeper.KeeperException$ConnectionLossException: KeeperErrorCode = ConnectionLoss for /
  3. git push和pull时遇到的两个问题及解决办法
  4. 产品体验营销或成广告业主们的新选择
  5. JavaScript 个人笔记3(详细BOMDOM)
  6. 跳出阿里之外做阿里之一------发布产品之外做阿里
  7. python爬虫 requests+lxml爬取前程无忧网之模拟浏览器登录
  8. 我来学网络——IIS中出现无效的应用程序池名称
  9. 拼多多根据关键词获取商品列表接口调用展示
  10. 关于安卓保存网络图片并更新图库