用分类变量拟合nls模型
分类变量拟合nls模型是一种统计建模方法，用于拟合非线性最小二乘模型。在这种模型中，自变量是分类变量，即具有有限个离散取值的变量。nls模型是一种广义线性模型，可以用于描述自变量和因变量之间的非线性关系。
分类变量拟合nls模型的优势在于可以捕捉到分类变量对因变量的影响，并且可以通过模型参数的估计来解释分类变量的效应。这种方法适用于具有多个分类变量的数据集，可以帮助我们理解分类变量在非线性模型中的作用。
分类变量拟合nls模型的应用场景包括但不限于以下几个方面：
- 医学研究：例如，研究某种药物对不同年龄组患者的治疗效果，将年龄作为分类变量拟合nls模型可以帮助确定不同年龄组的药物反应差异。
- 市场营销：例如，研究某种产品在不同地区的销售情况，将地区作为分类变量拟合nls模型可以帮助确定不同地区的市场需求差异。
- 社会科学：例如，研究某种政策对不同教育程度人群的影响，将教育程度作为分类变量拟合nls模型可以帮助确定不同教育程度人群的政策效果差异。
腾讯云提供了一系列与云计算相关的产品，其中与分类变量拟合nls模型相关的产品是腾讯云的机器学习平台（https://cloud.tencent.com/product/tiia）和数据分析平台（https://cloud.tencent.com/product/dla）。这些产品提供了强大的机器学习和数据分析功能，可以帮助用户进行分类变量拟合nls模型的建模和分析工作。
相关·内容
我有一些数据，我想用nls对数据的每个子集拟合一个非线性模型，然后用ggplot2将拟合的模型叠加到数据点的图上。具体而言，模型的形式是
y~V*x/(K+x)
你可能会认出是Michaelis-Menten。其中一种方法是使用geom_smooth，但是如果我使用geom_smooth，我就没有任何方法来检索模型拟合的系数。或者，我可以用nls对数据进行拟合，然后用geom_smooth绘制曲线，但如何知道geom_smooth绘制的曲线与我的nls拟合曲线相同？我不能将适合我的nls的系数传递给geom_smooth并告诉它使用它们，除非我可以让geom_smooth只使用数据的一个子集，
浏览 6提问于2017-10-17得票数 2
回答已采纳
我的问题与估计中的人口增长率有关。作为一个玩具示例，考虑一个玩具数据集df
structure(list(x= c(0L, 24L, 48L, 72L, 96L, 120L, 144L, 168L
), y = c(10000, 18744.0760659189, 35134.0387564953, 65855.509495469,
123440.067934292, 231377.002294256, 433694.813090781, 812920.856596808
)), .Names = c("x", "y"), row.names = c(NA,
浏览 7提问于2016-06-23得票数 3
回答已采纳
1回答
我对R相当陌生，我试图用nls函数拟合曲线。我首先用dgamma函数生成曲线y，然后我想用nls拟合它。这是我的玩具例子。
´´´
x <- 1:250
y <- dgamma(x,2,0.02)
df <- data.frame(x=x,y=y)
nls(y~ dgamma(x,a,b),data=df,start = list(a =2,b =0.4))
´´´
我所犯的错误是
numericDeriv中的错误(form[3L]，name(Ind)，env)：在计算模型时产生的缺失值或无穷大:警告消息:在dgamma中(x，a，b)：NaNs生成
我在这里做错什么了？
浏览 1提问于2020-12-13得票数 0
回答已采纳
1回答
我想用nls拟合一个非线性模型，但是我得到了一个non-numeric argument to binary operator错误：
df <- data.frame(task=rep(c("x", "y", "z"), each=6), time=rep(c(1:18)), y=c(1, 5, 9, 14, 17, 29, 2, 3, 1, 5, 2, 3, 29, 21, 18, 16, 14, 5))
OLS <- lm(y~ A + time*task, data=df)
我试图用nls函数得到类似的结果：
NLS <
浏览 10提问于2022-09-08得票数 1
回答已采纳
我想拟合一个渐近模型，如下所示：
y(x) = Asym + (R0 - Asym) * exp(-exp(lrc) * x)
可以使用nls()函数将模型拟合到R中：
y <- c(22.180,21.941,21.786,21.724,20.964,19.128,18.371,17.508,16.931)
x <- c(1,3,5,9,18,36,72,144,288)
nls(y ~ SSasymp(x, Asym, R0, lrc))
但是，如果我知道R0的值(例如R0 = 25)，那么我如何用自启动函数SSasymp来拟合这个模型？
浏览 0提问于2018-09-26得票数 0
我觉得我即将找到我的问题的答案，但不知何故，我就是做不到。我使用nls函数拟合了3个参数，使用了一个相当复杂的函数来描述精子浓度范围(x轴)内卵子受精成功率(y轴)(Styan的模型，)。拟合参数很好，但我无法使用predict函数绘制平滑的外推曲线(请参阅本文末尾的部分)。我猜这是因为我使用了一个不适合x轴的值。我的问题是，如何使用x轴上未拟合的参数，基于nls函数拟合的模型绘制平滑和外推的曲线？
下面是一个示例：
library(ggplot2)
data.nls <- structure(list(S0 = c(0.23298, 2.32984, 23.2984, 232.983
浏览 1提问于2013-04-06得票数 1
回答已采纳
我试图用R中的负指数模型来拟合一些数据，但它似乎不起作用，我也找不到错误或绕过它的方法。
t<-c(1,4,16)
y<-c(0.8,0.45,0.04)
tabla<-data.frame(t,y)
reg<-nls(y~exp(-b*t),tabla,start = list(t=0,b=0))
运行代码后，我会得到以下错误
Error in qr(.swts * attr(rhs, "gradient")) :
dims [product 2] do not match the length of object
浏览 0提问于2019-06-03得票数 0
回答已采纳
我有一个nls fitting任务，我想用R。我第一次尝试做这个，就像@Roland指出的
“关键是复杂的模型很难拟合。越是如此，数据就越不支持模型，直到它变得不可能。如果你有非常好的起始值，你也许能够适应这个模型。”
我同意@Roland的观点，但如果excel R 能做到这一点，为什么R不能这么做呢？
基本上，这种拟合可以用Excel的GRG非线性求解器进行，但过程非常耗时，有时拟合不好。(因为现实中有很多数据)。
这是我的样本data.frame。我想让每个set组与下面提供的模型相匹配，
set.seed(12345)
set =rep(rep(c("1",&#
浏览 2提问于2016-05-10得票数 1
回答已采纳
假设我有x值、y值和期望的y值f (来自一些非线性的最佳拟合曲线)。
如何在R中计算R^2？请注意，此函数不是线性模型，而是非线性最小二乘(nls)拟合，因此不是lm拟合。
浏览 1提问于2013-01-26得票数 17
回答已采纳
我使用R在对数曲线上用方程拟合数据：
Y=a* log(b * x)
我的数据如下所示：
#Creating example data
pre <- c(946116, 1243227, 1259646, 1434124, 1575268, 2192526, 3252832, 6076519)
post <- c(907355, 1553586, 1684253, 2592938, 1919173, 1702644,3173743, 3654198)
data <- data.frame(pre,post)
#Plotting data
ggplot(data,
浏览 0提问于2017-09-01得票数 1
我有一个每日降雨量(x)和观测值(y)的散点图，它看起来像x^-2图的右/正x值的一半，或者是以1/2为底的对数图。基本上，当x值非常低时，y值真的很高。X值越大，y就越小。但是y值下降的速度变慢了，并且y的值永远不会是负的。
下面是一个有代表性的示例：
rain <- c(1, 1.2, 1.3, 2.5, 3.2, 4.2, 5, 7, 7.5, 10.3, 11.7, 12.9, 14.1, 15, 15.5, 17.5, 18.3, 20, 20.2, 20.3, 25, 28, 30, 34, 40)
obs <- c(42, 44, 43.9, 43.5, 35,
浏览 18提问于2021-02-16得票数 0
回答已采纳
如何指定比简单的ID列更复杂的数据结构？
如果我有一个glmertree模型，我如何指定(例如)集群协方差检验中的交叉分类模型？
tree_1 <-
glmertree(
data = sim_dat,
formula =
performance ~ 1 + predictors |
(1 | student_id) + (1 | question_number) |
partitioning_variables,
family = 'binomial',
cluster = ???
浏览 6提问于2021-09-29得票数 1
我正在使用R中的nlstools包。我将一个模型拟合到我的数据中，并将其绘制得很好。我不确定我做错了什么，但我似乎不能从nlstools调用方法，因为我得到了一个错误，“对象'd‘找不到”。以下是代码，并注释了错误行：
Y=y_coll
X=x_ntips
d=data.frame(X,Y)
thisfit=nls(Y ~ a*X^b,data=d, start = list(a = .1, b = .1))
### Error in eval(expr, envir, enclos) : object 'd' not found
#a=nlsBoot (this
浏览 0提问于2012-12-12得票数 0
1回答
R中多条乙状结肠曲线的拟合与比较
、、、、
我想在一次拟合多个曲线，并比较他们的3个估计参数-渐近线，斜率和x0。下面是要建模的数据的理想图像：
大多数可搜索的页面都会使用方法来拟合一条曲线，如下所示：和here
在我的例子中，我想测试(统计)改变变量水平对乙状结肠的三个参数的影响。也就是说，当我拟合这个模型时：
model <- nls(y ~ asym / (1 + exp( -slope * (x – x0) ) ), start = c(…), data = my_data)
I想添加两个因素(例如，“factorA”和“factorB”)与lm()、lm()和x0每个术语的交互作用(lm()或glm() )，如下
浏览 3提问于2017-02-03得票数 6
2回答
考虑R中的非线性最小二乘模型，例如以下形式)：
y ~ theta / ( 1 + exp( -( alpha + beta * x) ) )
(我的实际问题有几个变量，外部函数不是逻辑的，而是涉及的更多；这个比较简单，但我认为如果我能做到，我的情况应该马上跟进)。
我想用(比如说)自然三次样条来代替术语"alpha + beta * x“。
下面是一些代码，用于在逻辑系统中创建一些具有非线性函数的示例数据：
set.seed(438572L)
x <- seq(1,10,by=.25)
y <- 8.6/(1+exp( -(-3+x/4.4+sqrt(x*1.1)*(1
浏览 4提问于2012-02-04得票数 6
回答已采纳
我注意到我可以使用nls拟合参数向量，如下所示。这让我来决定我想要的参数的数目。如下面的例子所示，我正在拟合y = k + a_1 x^2 + a_2 x^3 + a_3 x^3。我可以简单地改变初始值的数目，这改变了系数的数目来估计。
但是，这种方法不适用于nls2。它只是对待适合y = k + a_1 * x，三次！
我的问题是如何让nls2根据初始值(或类似的值)来确定参数的数量，就像nls的情况一样。
我对nls或类似的软件包没有太多的经验。所以，我正试着弥补这个问题。我猜nls2比nls有更多的能力。
数据与模型拟合
x <- c(32,64,96,118,126,144,15
浏览 2提问于2017-10-28得票数 0
回答已采纳
1回答
我是新的R，并试图拟合一个非线性最小二乘模型。功能是
R=a*SSB/(1+b*SSB)
变量SSB和R是数据集中的列。
初始参数为a = 0.005，b= 3.0e-10
浏览 4提问于2022-07-12得票数 -1
4回答
在R中，我使用nls进行非线性最小二乘拟合。那么，如何使用拟合提供的系数的值来绘制模型函数呢？
(是的，这是一个来自R相对新手的非常天真的问题。)
浏览 5提问于2012-03-29得票数 4
回答已采纳
我有一个年龄和长度的数据集，加上一些分类变量，包括性别和地点(2级因子)。我用nls()拟合了Gompertz模型
gompertz <- nls(Length~a*exp(-b*exp(-c*age)),
data=df,
start=list(a=155,b=0.4, c=0.1))
但也一直在努力加入一个分类变量--也就是说，我想比较两个地点之间的增长率。在此提出了一项解决办法：
这个模型工作-我得到了独立的a，b和c参数估计，对我的两个水平的位置因子。然而，我不知道如何绘制这个-我想绘制我
浏览 14提问于2022-07-06得票数 0
回答已采纳
1回答
我正在尝试使用Python进行分类。我有一些输入列(让k个变量)和一个输出列。
Let inputfeatures
= array([[ 0, 0, 0, ..., 0, 0, 0],
[ 1, 0, 0, ..., 0, 0, 0],
[ 1, 0, 0, ..., 0, 0, 0],
...,
[ 0, 0, 0, ..., 0, 0, 0],
[ 1, 0, 0, ..., 0, 0, 0],
[ 0, 0,
浏览 5提问于2021-03-06得票数 2
回答已采纳
我在R中使用lda函数来拟合模型。在对模型进行拟合后，我想使用拟合，然后将其输入到计算机程序中，根据输入进行分类。我只看到线性判别式和分组均值的系数。我想我需要截距和系数来得到每组的实际分数。
浏览 3提问于2021-07-09得票数 0
3回答
我有两个变量的散点图，例如：
x<-c(0.108,0.111,0.113,0.116,0.118,0.121,0.123,0.126,0.128,0.131,0.133,0.136)
y<-c(-6.908,-6.620,-5.681,-5.165,-4.690,-4.646,-3.979,-3.755,-3.564,-3.558,-3.272,-3.073)
我想找出更适合这两个变量之间关系的函数。
准确地说，我想比较三种模型的拟合情况：linear、exponential和logarithmic。
我在考虑将每个函数与我的值进行拟合，计算每种情况下的可能性，并比较AIC值。
浏览 1提问于2013-02-23得票数 4
回答已采纳
2回答
滑雪场随机森林
、、
我试着用sklearn的随机森林分类器包来拟合随机森林模型。但是，我的数据集由具有字符串值('country')的列组成。这里的随机林分类器不接受字符串值。它需要所有特征的数值。我想用一些虚拟变量来代替这样的列。但是，我搞不懂现在的特色重要情节会是什么样子。会有诸如country_India、country_usa等变量。如果我用R进行分析，如何才能得到国家变量的综合重要性。
浏览 0提问于2016-04-03得票数 0
回答已采纳
我是一个R新手，试图将植物光合作用的光响应曲线(饱和、曲线)拟合到专家接受的特定模型中。目标是获得Am、Rd和LCP的估计系数值。下面是我不断得到的错误：
Error in numericDeriv(form[[3L]], names(ind), env) :
Missing value or an infinity produced when evaluating the model
我已经多次改变了起始值，但仍然没有成功。帮助?提前谢谢你。下面是示例数据集。
photolrc= c(3.089753, 6.336478, 7.737142, 8.004812, 8.031599)
P
浏览 1提问于2015-10-22得票数 11
回答已采纳
1回答
在Windows10上使用RStudio v1.3.1093，R v.4.0.2。
我正在处理一个作业问题，并给出了一个数据集：
x <- cbind(c(8, 8, 14, 14, 14, 16, 16, 16, 18, 18, 20, 20, 20, 22, 22, 22, 24, 24, 24, 26, 26, 26, 28, 28, 30, 30, 30, 32, 32, 34, 36, 36, 42))
y <- cbind(c(0.49, 0.49, 0.45, 0.43, 0.43, 0.44, 0.43, 0.43, 0.46, 0.45, 0.42, 0.42,
浏览 4提问于2020-12-05得票数 0
回答已采纳
1回答
以下是我的数据集的摘录:相对频率和大小。从开环中可以看到，这是一个高斯分布。我使用R中的nls软件包来拟合一条非线性曲线。我的方程式是
或以纯文本形式：c * e^(-(x - z)^2/l)
我是怎么到这里来的
fit <- as.formula(y~c*(exp((-(x-z)^2/l))))
preview(fit_partial, data=mydata, start=list(c=0.005, x=mydata$x_values, z=130,l=2000))
起跑线似乎是合理的。所以我试着做一个非线性拟合
nls_fit <- nls(fit, data
浏览 2提问于2019-04-25得票数 1
回答已采纳
1回答
基于nls的指数回归
、、、
我很难将指数曲线与我的数据进行拟合。
这是我的密码：
x<-c(0.134,0.215,0.345,0.482,0.538,0.555)
y<-c(0,0,0.004,0.291,1.135,1.684)
plot(x,y)
estimates<- list(b1 = 0.1, b2 = 5e-7)
nlfit <- nls(y ~ b1 * (exp(x/b2)-1), start=estimates)
lines(x, predict(nlfit), col = 2)
但我得到了以下错误：
numericDeriv中的错误(form[3L]，name(Ind)，e
浏览 3提问于2021-04-18得票数 2
回答已采纳
1回答
我试图用ggplot2拟合一种蛋白质和两种浓度的小分子的缔合-解离自旋共振动力学数据。数据为。time变量表示时间(以秒为单位)，sample变量表示两种浓度(32 8nM和8nM)，values变量表示读出值。
我导入了data并运行了下面的代码来绘制：
# LINE PLOT
ggplot(data) +
geom_point(aes(x = time, y = values), size = 1, color = "black") +
geom_smooth(aes(x = time, y = values, color = sample), method = &
浏览 14提问于2022-10-10得票数 1
回答已采纳
我正在尝试使用mgcv包来拟合一个模型，并且我需要包含分类术语之间的交互。我知道如何实现连续变量和分类变量之间的交互，如下所示 s(x1, by=c1) 但我不知道如何将其更改为两个分类术语，并且我找不到任何东西。
浏览 72提问于2020-07-22得票数 0
1回答
我想用一个两项指数模型来拟合我的数据(x，y)，即
f(x) = a * exp(b * x) + c * exp(d * x)
本质上，我需要在R中复制Matlab的exp2模型类型，计算如下
f = fit(x, y, 'expo')
浏览 0提问于2016-03-08得票数 1
回答已采纳
我想用回归模型来拟合概率。我知道线性回归经常用于这个目的，但我有几个概率在0.0和1.0左右，并且想要拟合一个回归模型，其中输出被限制在0.0到1.0之间。我希望能够为模型指定正则化规范和强度，最好在python中这样做(但R实现也会有帮助)。我发现的所有逻辑回归包似乎只适合于分类，而这是一个回归问题(尽管我想使用logit链接函数)。我使用scikits为我的分类和回归需求学习，所以如果这个回归模型可以在scikits中实现，那就太棒了(在我看来，这是不可能的)，但是我很乐意在python和/或R中找到任何解决方案。
浏览 4提问于2015-11-21得票数 3
我使用nlsLM函数来拟合非线性回归。如何从nlsLM模型对象中提取HAT值和库克距离？
对于使用nls或nlreg函数创建的对象，我知道如何提取观察值的HAT值和库克距离，但我不知道如何使用nslLM获取它们。
有人能帮我解决这个问题吗？谢谢!
浏览 1提问于2016-02-28得票数 1
我有一个相当简单的问题。我一直在将一些统计分析代码从R转换到Python。到目前为止，我一直做得很好，但我被这条线卡住了：
nlsfit <- nls(N~pnorm(m, mean=mean, sd=sd),data=data4fit,start=list(mean=mu, sd=sig), control=list(maxiter=100,warnOnly = TRUE))
本质上，程序计算的是一组数据的非线性最小二乘拟合，即"nls“命令。在原文中，“代字号”看起来像一个"enye"，我不确定这是否有意义。
据我所知，Python的等价物是来自scipy.s
浏览 0提问于2019-03-10得票数 0
2回答
我试图通过使用ggplot2和包plotly来拟合指数模型，并进一步显示回归曲线，并得到一个R^2来检验模型假设。
这是我的资料
SR.irrig<-c(67.39368816,28.7369497,60.18499455,49.32404863,166.393182,222.2902192 ,271.8357323,241.7224707,368.4630364,220.2701789,169.9234274,56.49579274,38.183813,49.337,130.9175233,161.6353594,294.1473982,363.910286,358.32
浏览 0提问于2020-08-12得票数 0
回答已采纳
我正在尝试拟合模型y=c0*log(c1*x+1)，其中y和x是我的数据，c0和c1是模型的参数。我知道Matlab中的"fit“命令可以做到这一点，但我正在尝试在R中做到这一点。我想知道是否有人可以提供关于可以实现这一点的函数的建议。
model <- y ~ c0*log(c1*x + 1)
我期待一条对数曲线。
浏览 24提问于2019-11-02得票数 1
我正在处理文档分类问题，并尝试应用基本的RNN/LSTM模型。
我制作了多层LSTM模型，并用15个时期进行了训练，但我获得了大约53%的训练准确率。41%的测试访问率。当我用更大的时代训练时，训练acc。增加和测试访问权限。下降到33%左右。
我知道如果测试访问。低于训练ac.，它可能是过拟合的。而过拟合的原因是，模型甚至训练了训练数据集的噪声。但我认为，如果发生了过度拟合，训练是可行的。必须是相当高的(80-90%)和测试访问。这比这要低得多，因为过拟合意味着模型对于训练数据集来说是“过拟合”的。
现在我想知道为什么当我用15个时期训练时，我的准确率降低了10%，而训练准确率只有53%。它
浏览 0提问于2020-04-24得票数 0
1回答
因为我想单独使用nls模型，所以我对geom_smooth函数中的数据进行了拟合，并在ggplot之外进行了拟合：
library(ggplot2)
set.seed(1)
data <- data.frame(x=rnorm(100))
a <- 4
b <- -2
data$y <- with(data, exp(a + b * x) + rnorm(100) + 100)
mod <- nls(formula = y ~ (exp(a + b * x)), data = data, start = list(a = a, b = b))
data$fit &
浏览 2提问于2017-03-08得票数 3
回答已采纳
2回答
我想用指数衰减(或渐近曲线)的增长形式来拟合一个函数，这样：
Richness = C*(1-exp(k*Abundance)) # k < 0
我在上读过关于expn()函数的文章，但是根本找不到它(或者nls包)。我发现的只是一个nlstools包，但它没有expn()。我试过用通常的nls和exp函数，但我只得到了增加的指数.
我想像下面这样拟合这个图(用油漆绘制)，我不知道曲线应该稳定在哪里(丰富性= C)。提前谢谢。
浏览 3提问于2014-03-06得票数 1
回答已采纳
对于时间序列chuoie，我采用拟合模型AR(2)
e <- c(101,82,66,35,31,7,20,92,154,125,85,68,38,23,10,24,83,132,131,118,
90,67,60,47,41,21,16,6,4,7,14,34,45,43,48,42,28,10,8,2,0,1,5,12,14)
chuoie <- ts(e)
fit <- ar.yw(chuoie, order=2)
我想要模型拟合AR(2)的计算值：
fitted.values(fit)
然后我得到：
NULL
我不明白这个问题。
浏览 1提问于2016-02-20得票数 1
回答已采纳
2回答
我正在学习如何在R中使用nls函数，并且遇到了一些问题。目前，我只是试图重现一篇研究论文中的曲线。该模型拟合了1987年股市崩盘前的股市走势曲线。
我定义了一个函数func，如下所示：
func <- function(a,b,tc,t){
a+b*log(tc-t)
}
我这样调用nls：
nls1 <- nls(Y ~ func(a,b,tc,t), data2, start=list(a=0, b=1, tc=1466, t=1))
data2是由两列组成的数据框，一列是日期，另一列是值。共有1466行。
head(data2)
Date Y
1 1/4/82
浏览 0提问于2013-12-30得票数 2
1回答
对参数子集有二分影响的非线性模型
、、、、
我一直在研究R中的非线性模型，比如我们：
Y = Alpha1*time + Alpha2*sin(2*pi*time/Alpha3) + Alpha4*(-1)^time
我会列出伯努利变量是否会影响Alpha1*time和截取。这样的bernoulli变量可以是：
varia<-rep("BEFORE","AFTER"),each=30)
拥有
Y<- -2.5+rnorm(60)+2*sin(2*pi*time/8)+2.5*(-1)^time
time<-seq(1,60)
我发现nls可以拟合这个模型，但没有这个分类变量的影响，我还发
浏览 3提问于2014-04-28得票数 0
1回答
在滑雪板中是否支持直接使用熊猫的分类数据类型来拟合模型？据我所见，sklearn不支持这种数据类型，这是不幸的，因为分类数据类型既编码分类数据，又包含数据的映射方案。此外，分类编码纯粹是一个数据处理/处理问题，因此由Pandas来处理似乎更自然。
Note
我意识到在Pandas和sklearn中有几种编码分类变量的方法--这不是我要问的。
浏览 3提问于2015-06-15得票数 16
回答已采纳
1回答
我有一个具有数值和分类特性(变量)的数据集，我将所有的分类变量转换为虚拟数据，然后将训练和测试数据分割。
现在我的步骤，我想标准化的特点，然后拟合模型。
我应该把标准化应用于所有的特征，还是只适用于数字特征？
在这种情况下，在0到1的范围内使用MinMax标量更好吗？
浏览 0提问于2018-04-16得票数 1
1回答
我正在拟合一个具有6个预测值的多元线性回归模型(3个连续预测值和3个分类预测值)。残差vs.拟合图表明存在异方差，bptest()也证实了这一点。
我还计算了我的训练数据和测试数据的sqrt，如下所示：
sqrt(mean(sales_train_lm_pred-sales_train$SALES)^2)
3533.665
sqrt(mean(sales_test_lm_pred-sales_test$SALES)^2)
3556.036
我尝试拟合glm()模型，但仍然没有纠正异方差。
glm.test3<-glm(SALES~.,weights=1/sales_fitte
浏览 41提问于2019-11-11得票数 0
回答已采纳
3回答
我想对R中下列数据拟合等温线模型，最简单的等温线模型是这里给出的Langmuir模型，。下面给出了我的MWE，它会抛出错误。我想知道是否有关于等温线模型的R包。
X <- c(10, 30, 50, 70, 100, 125)
Y <- c(155, 250, 270, 330, 320, 323)
Data <- data.frame(X, Y)
LangIMfm2 <- nls(formula = Y ~ Q*b*X/(1+b*X), data = Data, start = list(Q = 1, b = 0.5), algorith = "port&#
浏览 4提问于2015-03-06得票数 3
回答已采纳
很简单的问题：
我想要拟合回归模型。我已经让熊猫知道这些变量是绝对的。样本代码：
bikes_reg['month'] = bikes_reg['month'].astype('category')
我是否需要创建虚拟变量，或者模型是否会在合适的情况下自动创建它们？
一般的经验法则是什么？每次使用分类变量时，是否都需要创建虚拟变量？
浏览 1提问于2017-08-25得票数 0
回答已采纳
1回答
我尝试将下面的函数y(x)=a*( 1+ (x^2)/(b^2) )^t拟合到一组特定的数据，其中a，b和t是希望通过拟合确定的常量。
例如，我尝试执行以下操作
len <- 24
x = runif(len)
y = x^3 + runif(len, min = -0.1, max = 0.1)
plot(x, y)
s <- seq(from = 0, to = 1, length = 50)
lines(s, s^3, lty = 2)
df <- data.frame(x, y)
m <- nls(y~a*( 1 + (x^2)/(b^2) )^t, data
浏览 1提问于2015-12-10得票数 1
1回答
我已经建立了一个机器学习分类器，使用Sklearn和熊猫作为我的主要工具。现在，模型的一个输入特性是country (给国家代码(如US)写信)。我用pd.get_dummies函数拟合了一个模型。
现在我想对数据进行推断，但是在过去的几个月里，我的数据集中还没有出现一些国家，所以pd.get_dummies函数与拟合的模型不一致。
这怎么处理呢？
浏览 0提问于2017-12-18得票数 4
回答已采纳
2回答
我正在使用scikit-学习做回归，我的问题是以下。我需要对几个参数(向量)进行回归。这在一些回归方法(如ensemble.ExtraTreesRegressor和ensemble.RandomForestRegressor )中很好。实际上，我们可以给出一个向量作为目标来拟合上述两种回归方法的模型(fit(X,y)方法)。
然而，当我尝试使用ensemble.GradientBoostingRegressor、ensemble.AdaBoostRegressor和linear_model.SGDRegressor时，分类器无法拟合模型，因为它期望一维值作为目标( fit(X,y)方法的y参数
浏览 7提问于2014-02-04得票数 21
1回答
我试图在非线性回归之后拟合一条曲线。
这是我的数据集：
stem_diameter <- c(15, 15, 16, 17, 19, 23, 23, 24, 24, 25, 25, 26, 27, 28, 29, 30, 30, 32, 32, 33, 34, 34, 35, 36, 36, 37, 38, 40, 41, 41, 42, 42, 46, 48, 48, 49, 51, 54, 55, 60)
total_biomass <- c(0.25, 0.65, 0.40, 0.40, 0.65, 0.60, 0.20, 0.60, 0.50, 0.35, 0.70, 0
浏览 0提问于2018-08-07得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
