e****e 发帖数: 677 | 1 关于卡方拟合度检验 chi-square goodness of fit
我有几万个数据,要和另外一组数据进行卡方拟合度检验
请问怎样确定bin size啊?
且这些bin的边缘或者说中心如何确定阿?
谢谢 |
|
b*****e 发帖数: 2511 | 2 对一组数据进行了Lorentz fitting, 拟合的误差应该怎么算的? Lorentz
distribution里的Gamma应该只是和这个distribution有关和拟合 误差没有任何关系吧
? |
|
c******x 发帖数: 350 | 3 对的,y_i属于某个区间。比如说y_i是收入,x是学历,家庭背景等等。对应于各个学
历,收入应该是个区间,比如说100个高中学历,收入区间为20K到30K,100个大学学历
,收入区间为30K到40K,等等。如果只考虑学历的话,那就是一元函数拟合,没有问题
。但是对于多元函数,比如加上家庭背景,好像无法简单拟合。 |
|
e*********g 发帖数: 77 | 4 不了解PCA
想用一堆变量拟合一堆数据,结果没办法拟合。
于是考虑减少变量个数,
但一般求特征值只是在变量空间里找正交主轴,特征值是variance,并不是按照定义域
所引起的值域变化量排序。
所以问题是,
怎样才能剔除假的不作用的变量,而找到函数值对其依赖最大的那些变量?
谢谢。 |
|
t**m 发帖数: 2262 | 5 英媒披露中国与多国拟合建欧亚高铁网
http://www.chinareviewnews.com 2010-03-11 09:24:38
中评社北京3月11日电/据英国每日电讯报9日报道,中亚高铁最终可能会与欧洲铁
路网连接,届时旅客从伦敦王十字火车站只需要48小时就可以抵达北京。报道说,印度
和欧洲的几个国家已经与中国政府就上述项目展开谈判,预计未来10年内中国时速高达
350公里的高铁动车就能在欧亚大陆上疾驶。
多国铁路网将相连
长江商报报道,中国工程院院士、隧道及地下工程试验研究中心主任王梦恕日前接
受媒体采访时说,中国计划筹建一个连接中国与中亚国家的高速铁路网,并最终可能与
欧洲铁路网连接。他透露说,中国目前已经开始与相关国家展开谈判,这一高速铁路网
将把中国乌鲁木齐与中亚的哈萨克斯坦、乌兹别克斯坦和土库曼斯坦等国连接起来,预
计覆盖17个国家。王梦恕提到:“我们的目标是使列车运行速度达到每小时350公里。
”
目前,伊朗、巴基斯坦和印度正就铺设高铁线与中国谈判。中俄已就铺设横跨西伯
利亚的高铁路线达成一致。王梦恕说,“随着快捷的高速铁路运输的引进,人们将在西
部建立矿山 |
|
d******a 发帖数: 32122 | 6 重磅!首都医科大学与北京航空航天大学拟合并
2015-12-26 来源:健康界 A- A+
编者按:12月26日,一则“首都医科大学与北京航空航天大学合并”的消息刷爆医疗圈
。其导火索是一份网传公告,“公告”称首医与北航已于12月23日召开会议,并签署紧
密合作协议。但细观网传公告原文,疑点重重。其中最明显的一条是,“公告”中出现
的北航党委书记周天行已于1992年离世,现任北航党委书记为张军。健康界立即对此事
进行核实。
据健康界核实,首医与北航合并商谈接近尾声,与此同时,北航还在与北京其他一些市
属大学商谈合并事宜。北航此举旨在快速成长为综合性大学,从而尽快实现习主席提出
的构建世界一流大学的目标。但网传消息中,两校于12月23日召开会议并签署紧密合作
协议之事并不属实。
另据消息人士透露,此公告之所以于12月26日被热炒,是因为前一日,教育部领导与北
航领导再次商谈合校事宜。
两校合并将如何演变?首都医科大学命运如何?健康界将持续关注!
网传公告如下:
12月23日,北京航空航天大学首都医学部筹建工作会议在北京航空航天大学沙河校区新
主楼会议中心隆重举行。原全国人大常委会副委员长... 阅读全帖 |
|
s********y 发帖数: 473 | 7 晕。。。钱老做了那么久还不知道怎么拟合数据么。。。 |
|
发帖数: 1 | 8 高阶拟合都会出现震荡,你不会连这都不知道吧
跟python没关系 |
|
|
发帖数: 1 | 9 搞多项式拟合也好意思说是机器学习吗? 而且
overfitting了
:polyfit.py:14: RankWarning: Polyfit may be poorly conditioned
: p = np.polyfit(x,y,27) |
|
发帖数: 1 | 10 而且还可以把所有武汉病患的微信购买记录从12中旬-1月中旬也做成地理坐标点状图,
再和早期患病家庭地理位置点状图做一个PCA拟合,多么好的一篇paper啊 |
|
f*****y 发帖数: 1997 | 11 我本来猜想这个形式
y=a(1- exp(-b x)),但是好像EXCEL拟合不出来这个公式?
你觉得是什么形式?
除了EXCEL还有什么软件好? |
|
z*****4 发帖数: 45 | 12 signal不是正态分布。如果把signal vs t 画出来,是exponential decay。这个问题
其实是把这12个数据点,fit成bi-exponential curve。
来做一下simulation: 假设a=0.5, b=0.5, R1=0.002, R2=0.1,根据已知的12个t,可以
算出12个Signal.然后Signal再加上随机的noise,在matlab可以用randn实现。这样我
们有了12个t,12个Signal(加了noise的),来fit a, b, R1和R2。在matlab里,可用fit
这个function来做。但每运行一次fit,得到的a,b,R1和R2的值,都是不一样的,经常
会和原始值a=0.5, b=0.5, R1=0.002, R2=0.1相差很远。
一篇文章上写着(是一幅图):用Bayesian方法,算出来的R1和R2是一个分布,像是正态
分布,R1的值是横坐标,次数是纵坐标。然后取 出现次数最多的R1最为最后的拟合值
。可惜这篇文章根本不告诉你这个过程是怎么实现的,开头是Bayesian概率论的基本公
式,后面直接上结果。所... 阅读全帖 |
|
M********t 发帖数: 5032 | 13 引力波
带人类倾听星辰大海之声
新浪科技
ugmbbc
8小时9分钟前
在一波又一波传言后,终于!北京时间2016年2月11日23:40左右,激光干涉引力波天
文台(LIGO)负责人、加州理工学院教授David Reitze宣布,人类首次发现了引力波。
发布会上究竟揭开了哪些谜团?与之前的传言有何不同吗?中国科学院国家天文台黑洞
来客团队(团队主要成员苟利军 研究员为国家天文台恒星级黑洞研究创新小组负责人
,而这次被探测到的引力波正是来自双黑洞系统),全面解读这场激动人心的发布会。
正如发布会所言,在被预言将近百年、苦苦追寻几十年之后,首个位于地球之外13亿光
年的引力波源GW150914被人类直接探测到,这是一个值得纪念的伟大时刻,一个新时代
的序幕正在拉开——地球人,欢迎你来到引力波时代!
激动人心的发布会结束了,其重点内容可以被简单归纳为三点:
(1)引力波终于被探测到了。
(2)引力波产生于两个恒星量级黑洞的合并(merger)。
(3)引力波是美国的激光干涉引力波天文台(英文简称LIGO)发现的。
接下来,让我们逐个分析和解释一下以上三点,从而对这项具有划时代意义的科学发现
... 阅读全帖 |
|
v***a 发帖数: 903 | 14 拟合了做什么用?
你知道数据之间大概是啥形式的关系不?我记得matlab 那个 fit 的tool 可以自己加
函数的。
你用样条的话也就能得到光滑的曲线s吧,还能干嘛用? |
|
T*********s 发帖数: 20444 | 15 Sigmaplot里面有不少内置函数
但是需要你选择
我觉得不会有软件能自己选择拟合函数 |
|
a***e 发帖数: 27968 | 16 没模型拟合出氢光谱线的那种奇葩是很少见
: 函数背后都有模型,
: 你模型都没有fit出来也没用啊
: 还不如搞你之前那个time series
: ,(
|
|
B******u 发帖数: 23763 | 17 知道一组数值,需要寻找其拟合样条曲线函数。
通过软件最好拉,
比如matlab,或者其他啥。
包子谢哈。
btw,
千万别发我信箱,打不开的。
发我邮箱c*******[email protected]可以
或者版面指导哈。
多谢多谢。 |
|
l*****d 发帖数: 359 | 18 名字不要管他, 道理很简单 ,一次导数为0的点就是最小值。那个E其实就是观测值z
和你拟合的f(x,y)之间的的均方差函数,把所有的x[m],y[n],z[m,n]带入E,你就得到
一个以系数aij为变量的多元函数E(a11,a12,a13,a21,...)。n个系数的话,这个函数就
有n个变量(为了方便,把你的系数aij换成x1,x2...,xn)。 分别对这n个变量求导,你
就得到n个线性方程。这个线性方程组中xi的系数可以写成n阶矩阵。 求xi向量就等于
求解一个矩阵方程。自己编程求矩阵的每个元素吧, 并不难。写个matlab function A
= get_coef(x[m],y[n],z[m,n])。A是N阶方阵,N为你系数个数,在你的case中N=9 |
|
d*****l 发帖数: 8441 | 19 伪逆(pseudo-inverse),因为精确拟合不可能,只能有最小二乘意义上的精确度。
确实是X^(+)=(X'X)^(-1)X'。 |
|
C*********m 发帖数: 213 | 20 把peak intensity or peak volume 读出来,把浓度换算出来,自己用数学软件拟合就
是了。他们的公式不是写在script里面就是另外有manual,自己找一下。每个体系不一
样,他们的公式也未必适用你的 |
|
k***n 发帖数: 383 | 21 知道MODEL是SPHERICAL HARMONIC(球谐振子),
现在有数据了,
谁知道有没有算法来拟合?
就一个参数,就是AMPLITUDE.
因为L,M都已知了. |
|
g******i 发帖数: 32 | 22 隔壁有人提到用excel, 我试了一下,发现可以在绘图结果上单击点,然后按右健选'add
trendline',从linear, exp, log等选项中选一种,并选择显示公式选项就可以得出公式
。只是拟合的不是很精确。有人提到使用datafit, 不知道效果如何。 |
|
g******i 发帖数: 32 | 23 汇报一下,datafit 真的很好用,是个不错的拟合软件。:)
我试了一下,发现可以在绘图结果上单击点,然后按右健选'add
式
with |
|
Y*e 发帖数: 63 | 24 我有几组数据,要拟合成f(x)=a*exp(-x/b)+c形式的函数,如何能在matlab中简单实现
?我找了半天没找到合适的函数。哪位能指导一下,多谢! |
|
g*********e 发帖数: 42 | 25 他是指用什么方法来优化函数中的参数吧?
如果是的话,建议你查查optimization;
以前常用matlab,现在忘了,不过记得有
个函数可以实习这个.
发信人: rossby (五十岚已夜), 信区: Computation
标 题: Re: matlab参数拟合
发信站: BBS 未名空间站 (Thu Apr 26 17:56:02 2007), 转信
log (f(x) -c) = log a - x/b |
|
c*u 发帖数: 916 | 26 我有一组实验数据,
y = y_high, if x > thr
y = y_low, if x < thr
y_high, y_low, thr 是未知的常数,thr是所需要得到的东西.好像3D bruteforce
search 寻找最低点不是很有效率.
举个通俗的例子,假设x是温度,y是水的体积,那么Thr是熔点.但是实验数据中
y的测量有很大noise, 所以必须拟合.
谢谢, |
|
c****l 发帖数: 14 | 27 小弟初学,许多地方不懂,想实现附件的积分并进行参数d,n拟合,不知哪位大侠能给出个
程序,谢谢。
x的范围[x,x+80]
x=[0;12;25;63]
y是积分的结果对应的项 [239;265; 276; 294]
多谢各位达人赐教 |
|
l****r 发帖数: 130 | 28 我这里有一个函数
a*y^(1/2)*x^(-1/2)+y+b*y^(3/2)*x^(1/2)=1.
然后,我还有一组(x,y)的数据,这组数据画出来类似抛物线。怎么通过拟合来求出
未知系数a,b呢。请求各位大师显灵点拨。多谢了。 |
|
s***m 发帖数: 6197 | 29 我换个问法吧
什么软件能做这种几个高斯分布叠加在一起的数据拟合?
谢谢! |
|
g******s 发帖数: 733 | 30 x=[0.1 0.2 0.4 0.6 0.8 1 2 4 6 8 10 20 40 60 80];
y=[200 100 50 35 28 22 12 6.6 5 4 3.6 2.5 1.9 1.8 1.8];
loglog(x,y,'o-')
set(gca,'FontSize',14)
请问什么函数(不能分段)可以拟合附图中的曲线?
先谢了! |
|
r**e 发帖数: 339 | 31 我在做一个吸附的问题,里面有各种平衡条件(平衡系数),有的数据点是(pH,加入酸碱量)
F(Ph,加入酸碱量)这个隐函数是由一堆平衡条件和物质守恒得到的方程组决定的. 而这
个方程组是没有解析解的. 很头疼不知道怎么拟合相应的系数(及平衡常数).
如果只是一个隐函数,我知道怎么办.但是现在是一个方程组导致的隐函数,没有任何形
式. (方程组里有其他变量,如果能解出,是应该被消掉的)
不知道我应该用到什么知识 ,-- 理论的方法和数字的软件等
Thanks a lot. |
|
f*******l 发帖数: 1239 | 32 小弟初学者, 用lsqcurvefit 拟合参数。但是总是出现Too many input arguments.
各种检查,都没有问题。不知道各位能不能给些意见。谢谢 |
|
c******k 发帖数: 1140 | 33 画的这一条曲线类似于三段直线连接在一起,但是他们有不同的斜率,就是想知道每一
段的斜率,我用的是Origin 8, 在 Analysis-->Fitting-->Fit Linear, 本想分段拟
合直线,但是每次都是把三段直线当成整体,作的直线拟合,指点一下,多谢了!!! |
|
x****g 发帖数: 2000 | 34 放在三个不同的column里呗,各自拟合,
最后再把它们的marker改成一样的,这样看起来就对了 |
|
m*1 发帖数: 5 | 35 ===问题是===
先考虑二维情况,已知一系列点(x1,y1),(x2,y2),...(xn,yn),但是不知道他们背后的
函数是什么。如何在某种误差限度内,找到一种带宽比较小的函数较好地拟合这些点。
所谓带宽较小,举例就是,(1,4) (3,-2) (4,10) 可以找到穿过这三个点的
曲线很多,希望找到穿过y=0坐标轴次数较少的解答。
=========
比较急,在线等,感谢先 :) |
|
|
g**r 发帖数: 48 | 37 几组数据,很多变量,要拟合比较复杂的函数,请推荐哪个软件最好用?
谢谢 |
|
o******e 发帖数: 4 | 38 如函数 y(x) = p1 * ( 1 - erf( x/p2)^(1/2) ) )
其中p1,p2为待定参数,erf是error function,已知若干实验数据点(x,y),问题出现在p1
数值约为10^8(真实值可能会相差两个数量级左右),而P2为10^(-8)(真实值与此相比可
能会相差两个数量级左右),我用matlab的curve fitting 工具箱拟合时都说出现错误,
给的初始值有问题,怎么办,是否考虑变量替换,如p1'=10^8/P1,p2=10^(-8)/p2??谢谢` |
|
c*u 发帖数: 916 | 39 我有一组实验数据,理论上应当符合以下规律:
y = y_high, if x > thr
y = y_low, if x < thr
y_high, y_low, thr 是未知的常数,thr是所需要得到的东西.好像3D bruteforce
search 寻找最低点不是很有效率.
举个通俗的例子,假设x是温度,y是水的体积,那么Thr是熔点.但是实验数据中
y的测量有很大noise, 所以必须拟合.
谢谢, |
|
h***s 发帖数: 226 | 40 看来这里的“数学家”是不会帮你解决这么具体的问题的。
方法可能有几个
1.弄一个数值计算工具,比如matlab等,里面有回归的部分,选择需要拟合的函数,和
xy数值,得到ab。但是否能有你需要的函数种类,不得而知;
2.弄一个已经编写好的非线性回归程序控件,然后写你的输入输出函数调用即可:http://download.csdn.net/source/142460;
3.自己写程序,比如利用非线性最小二乘方法编写程序,不过这个程序还不是那么容易
编写的,如果谁白手起家不看任何资料就能写出来,是比较牛的。 |
|
g****t 发帖数: 31659 | 41 matlab有专门的非线性最小二乘法的函数.
查一下help lsqnonlin
或者用fminsearch
看来这里的“数学家”是不会帮你解决这么具体的问题的。
方法可能有几个
1.弄一个数值计算工具,比如matlab等,里面有回归的部分,选择需要拟合的函数,和
xy数值,得到ab。但是否能有你需要的函数种类,不得而知;
2.弄一个已经编写好的非线性回归程序控件,然后写你的输入输出函数调用即可:http://download.csdn.net/source/142460;
3.自己写程序,比如利用非线性最小二乘方法编写程序,不过这个程序还不是那么容易
编写的,如果谁白手起家不看任何资料就能写出来,是比较牛的。 |
|
a***g 发帖数: 2761 | 42 【 以下文字转载自 Computation 讨论区 】
发信人: akoug (xuantai), 信区: Computation
标 题: 请教一个曲线拟合的问题
发信站: BBS 未名空间站 (Fri Dec 4 12:27:58 2009, 美东)
这个问题的表述是我自己想的:
假设有n个数据:x1,.....,xn;
现在把1,.....,n视作自变量,把x1,.....,xn视作因变量
可以假设xi=f(i)
如果x1,....,xn真的服从这样的曲线f
那么有什么方法可以比较准确的估计f在{1,....,n}这n个点上的导数呢
谢谢 |
|
e****e 发帖数: 677 | 43 我有几万个数据,要和另外一组数据进行卡方拟合度检验
请问怎样确定bin size啊?
且这些bin的边缘或者说中心如何确定阿?
谢谢 |
|
r**e 发帖数: 339 | 44 我在做一个吸附的问题,里面有各种平衡条件(平衡系数),有的数据点是(pH,加入酸碱量)
F(Ph,加入酸碱量)这个隐函数是由一堆平衡条件和物质守恒得到的方程组决定的. 而这
个方程组是没有解析解的. 很头疼不知道怎么拟合相应的系数(及平衡常数).
如果只是一个隐函数,我知道怎么办.但是现在是一个方程组导致的隐函数,没有任何形
式. (方程组里有其他变量,如果能解出,是应该被消掉的)
不知道我应该用到什么知识 ,-- 理论的方法和数字的软件等
Thanks a lot. |
|
U********a 发帖数: 75 | 45 最近在学习有理拟合,看到了Stoer-Bulirsch算法和vector fitting算法,版上有没有
对这两种算法了解的,解释下两种方法的优劣,谢谢! |
|
f*******r 发帖数: 204 | 46 excel可以做拟合,中文叫趋势线,在图上可以显示equation。 |
|
h****l 发帖数: 7290 | 47 理论上,拟合从来都是模拟原数据趋势,永远不可能推出正确公式,呵呵 |
|
g******i 发帖数: 32 | 48 多谢,我已找到你说的该功能,虽然不是很精确。没有可以更精确拟合的了么?有人提到
datafit,不知道效果如何。 |
|
p*****u 发帖数: 202 | 49 请教各位大侠:
我有四组数据,四个方程,一一对应。
但四个方程要同时share一个参数
如何同时拟合这四组数据求出这个参数?
比较傻,不会,请赐教,谢谢:) |
|
p*****u 发帖数: 202 | 50 谢谢大侠们的建议,我的确不是在“玩”数值计算,只是需要分析我们组的实验数据,
最后的目标是做蒙特卡罗simulation。我们是实验的group,平时也没有做过这样的拟
合,从头学起恐怕有点累,所以想走个捷径——像我们这样的恨不得有人给个
application我们直接把数据输进去得结果呢。只是老板雇了俺,那不管是啥,俺是会
也得会不会也得弄会了。
我记得有一次老板让我写程序做常微分方程计算,问我以前做没做过,我说没有(平时
尽和各种仪器搏斗了没怎么写程序),老板瞪大了眼睛问我:你没用Matlab解过常微分
方程?我很想翻白眼的说,没解过那也不是很稀奇啊,还不会学嘛。
好了说远了,哪位大侠能再给些提示就拜谢了,origin的拟合我们用在原始数据上了,
但是没法进一步产生随机数(可能是我不会用)进行蒙特卡罗simulation。所以还需要
用matlab的。
好。 |
|