-
0 votesanswersviews
Matlab多指数曲线拟合到同一组系数
我'm trying to fit multiple exponential curves to pairs of data using the same set of coefficients for ' exp1'拟合方法,会产生这样的东西: 示例数据如下: timeToFit = [24 48;48 72]; dataToFit = [150000 230000;105000 430000]... -
0 votesanswersviews
LMFIT - Extrackt变量形成Fitfunction
我写了一个简短的代码,以适应光学发射规格的光谱 . 因此,我将一个VoigtModel安装到Peak中,将LinearModel安装到背景中 . 喜欢这样: mod=VoigtModel() pars = mod.guess(y, x=x) out = mod.fit(y, pars, x=x)` 和 mod=LinearModel() pars = mod.guess(y, x=x) fun... -
1 votesanswersviews
如何为nlm提供渐变(或粗麻布)?
在非线性最小化函数 stats::nlm 的帮助文件中,它声明其中一个参数 f 是(强调我的): 要最小化的函数,返回单个数值 . 这应该是一个函数,第一个参数是长度为p的向量,后跟由...参数指定的任何其他参数 . 如果函数值具有称为渐变的属性或梯度和粗麻布属性,则这些属性将用于计算更新的参数值 . 这是否意味着渐变和hessians由以下命令提供: attr(f, 'gradient') ... -
5 votesanswersviews
在python中为一对发行版生成MLE
好的,所以我当前的曲线拟合代码有一个步骤,使用scipy.stats根据数据确定正确的分布, distributions = [st.laplace, st.norm, st.expon, st.dweibull, st.invweibull, st.lognorm, st.uniform] mles = [] for distribution in distributions: par... -
12 votesanswersviews
矢量自回归模型拟合与scikit学习
我正在尝试使用scikit-learn中包含的广义线性模型拟合方法拟合向量自回归(VAR)模型 . 线性模型的形式为y = X w,但系统矩阵 X 具有非常奇特的结构:它是块对角线,并且所有块都是相同的 . 为了优化性能和内存消耗,模型可以表示为 Y = BW ,其中 B 是来自 X 的块, Y 和 W 现在是矩阵而不是向量 . LinearRegression,Ridge,RidgeC... -
2 votesanswersviews
如何在zelig中获得多次推算数据的模型拟合度量(AIC,F-统计量)?
继an earlier post之后,我有兴趣学习如何使用乘法估算数据(使用Amelia创建)在zelig中获得统计模型相对质量的常用度量 . require(Zelig) require(Amelia) data(freetrade) #Imputation of missing data a.out <- amelia(freetrade, m=5, ts="year&quo... -
0 votesanswersviews
R:差异变异函数拟合,克里金结果差
我试着在雅加达湾做一个克里金 . 我有一套测量点,有适当的坐标和属性(pH,盐度,......) 为了做克里金法,我首先需要为变差函数找到一个模型 . 当我使用"variogram"函数输出不完美但是应该没问题,但是当我尝试拟合变异函数时,我得到一个警告信息:在fit.variogram中(ph.vgm,model = vgm(0.12,"Sph",0.1)... -
0 votesanswersviews
mselect函数的drc包在函数内部不起作用
我想使用drc包中的mselec函数来选择剂量反应数据集的最佳模型 . 但是,如果在函数内部使用mselect,则mselect不起作用 . 以下代码有效: library(drc) ryegrass.m1 <- drm(rootl~conc, data = ryegrass, fct = LL.4()) mselect(ryegrass.m1,list(LL.5(), LN.4(), W...