63得票1回答
在R公式中,为什么我必须在幂项上使用I()函数,比如y ~ I(x^3)?

我试图理解波浪线运算符及其相关函数的使用。我的第一个问题是为什么需要使用I()来指定算术运算符?例如,这两个绘图生成不同的结果(前者生成直线,后者生成期望的曲线)。x <- c(1:100) y <- seq(0.1,10,0.1) plot(y~x^3) plot(y~I(x^...

43得票3回答
R线性回归公式中的大写字母"I"代表什么意思?

我一直无法找到这个问题的答案,主要是因为使用单独字母(例如“I”)进行任何谷歌搜索都会导致问题。 在这样的模型中,“I”有什么作用?data(rock) lm(area~I(peri - mean(peri)), data = rock) 鉴于以下内容无法正常工作:lm(area ~ (pe...

22得票5回答
无法理解sklearn中的PolynomialFeatures。

需要在sklearn的Polynomial Features中寻求帮助。它在一个特征方面运作得相当不错,但是每当我添加多个特征时,它还会输出一些除了值的幂次方之外的值。 例如:对于这个数组,X=np.array([[230.1,37.8,69.2]]) 当我尝试时X_poly=poly.fit...

13得票2回答
整数系数多项式的快速分解

我想在整数环上快速分解多项式(原多项式具有整数系数和所有因子具有整数系数)。 例如,我想将4*x^6 + 20*x^5 + 29*x^4 - 14*x^3 - 71*x^2 - 48*x分解为(2*x^4 + 7*x^3 + 4*x^2 - 13*x - 16)*(2*x + 3)*x。 应该...

13得票2回答
多项式在质数模意义下的根

我正在寻找一种快速算法,可以在有限素域内查找一元多项式的根。 也就是说,如果 f = a0 + a1x + a2x2 + ... + anxn (n > 0),那么对于给定的质数p,需要找到所有满足 f(r) = 0 mod p 的 r < p。 我发现了Chien搜索算法http...

12得票1回答
SymPy有限域多项式

import sympy as S F = S.FiniteField(101) 当我调用f = S.poly(y ** 2 - x ** 3 - x - 1,F)时,我会得到以下错误: 'FiniteField'对象没有属性'is_commutative' 但是有限域根据定义是...

11得票2回答
使用np.polyfit在三维空间中拟合多项式

我有一个数据数组,维度为(N,3),其中 N 是某个整数,它指定了在三维空间中粒子的轨迹,即每行条目是粒子的(x,y,z)坐标。该轨迹是平滑和简单的,我希望能够将多项式拟合到这些数据上。 我可以使用np.polyfit来处理只有(x,y)坐标的数据:import numpy as np #...

11得票2回答
SymPy:多项式中高阶项的删除

使用Sympy,假设我们有一个表达式f,它是符号"x"(和其他可能的符号)的多项式。 我想知道是否有一种有效的方法可以删除f中高于某个整数n阶的所有项。 作为一个特例,我有一个非常复杂的函数,但我只想保留x的二次项。有什么更有效的方法吗? 显然,不太有效的方法是对于每个小于n的m,取m阶...

11得票1回答
如何使用sympy将一个表达式重写为另一个表达式的形式

编辑:我不是在问如何解出给定变量的方程(如此所谓的重复问题),而是如何表示一个表达式以另一个表达式为基础,就像问题中指定的那样。 我认为这个“重复”的问题有一个误导性的标题。 我对SymPy非常陌生。我有一个表达式,一旦用另一个表达式表示出来,就会变得非常好看。 问题是我不知道如何“强制”使...

10得票2回答
多项式比值逼近

我正在尝试对我的数据集进行多项式拟合,看起来像这样(完整的数据集在帖子末尾): 理论预测曲线的公式为: 当x在0到1之间时,它看起来像这样: 但是当我尝试用R做一个线性模型时:mod <- lm(y ~ poly(x, 2, raw=TRUE)/poly(x, 2)) 我...