在为我的线性模型创建qqPlot时,我遇到了一些问题。我的其他“测试”都很好。
首先,下面是我使用的包:
library(astsa)
library(TSA)
library(RcmdrMisc)
library(forecast)
library(readxl)
这是我的数据:
> data
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
2018 42 31 34 44 34 58 29 57 33 60 43 24
2019 30 22 27 37 29 69 39 36 52
我有一个矩阵mEps,它的形状是(10, 1042),其中10是资产的数量,1042是数据点的数量。我想显示每个资产的Q-Q图,这样我就可以绘制:
for i in range(iN):
sm.qqplot((mEps[i,:]), fit = True, line='q')
然而,我得到了10张Q-Q情节的图片。我想把它们放在一个数字中,所以我有以下代码:
fig, axes = plt.subplots(nrows=4, ncols=3, figsize=(15,10))
ax= axes.flatten()
for i in range(iN):
sm.qq
我正在尝试创建一个Q-Q图,以测试我的数据是否可以使用下面的命令用Weibull发行版建模。
qqplot(x,'weibull')
中使用数据
x =c(3.367, 0.769,0.8,1,1.2)
我不断地被错误所困扰
"In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion"
却找不出原因。这是否意味着我不能将威布尔分布与我的数据相匹配?如果有人能帮我指出为什么这不起作用的话,我将非常感激。
我刚刚在运行Python3.6的IPython笔记本(Jupyter)中测试了statsmodels.api.qqplot函数,在一列中得到了两个相同的图(我只要求了一个)。怎么一回事?
import statsmodels.api as sm
test = np.random.normal(0,1, 1000)
sm.qqplot(test, line='45')
我一直在尝试用python做一个用于泊松分布的QQ图。这是我到目前为止所知道的:
import numpy as np
import statsmodels.api as sm
import scipy.stats as stats
pois = np.random.poisson(2.5, 100) #creates random Poisson distribution with mean = 2.5
fig =sm.qqplot(pois, stats.poisson, line = 's')
plt.show()
每当我这样做时,我都会得到"Attribute
这是我绘制并保存图形的代码的一个示例:
我正在使用Python3.7.4和matplotlib==3.0.3。
import matplotlib.pyplot as plt
import pandas as pd
from yahoo_fin import stock_info
import statsmodels.api as sm
brk_data = stock_info.get_data("BRK-A")
with plt.style.context('dark_background'):
fig, ax = plt.subplots(f
我正试着开始使用statsmodel包来生成qqplots。我是从源代码中安装的,使用的是使用python 3.6的主分支。对于我想要做的事情,我想做一个qq图,比较两个不同样本大小的数据分布。我只是试图运行文档中的示例代码,但是它会引发一个关于不同样本大小的错误。
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
import statsmodels.api as sm
from statsmodels.graphics.gofplots import qqplot
# example 6
x = n
我已经从一个包含两个组的数据集创建了一个qq图(带有beta分布的分位数)。为了形象化,哪个点属于哪个组,我想给它们着色。我尝试了以下几点:
res <- beta.mle(data$values) #estimate parameters of beta distribution
qqplot(qbeta(ppoints(500),res$param[1], res$param[2]),data$values,
col = data$group,
ylab = "Quantiles of data",
xlab = "
我正在尝试为dataset的各个列在for循环中创建子图。我正在使用的加州住房数据集从滑雪。因此,这里有4列,我想在一个子图中为每个列显示三个数字。我已经提供了我尝试过的代码。有人能帮我解决这个问题吗?我们能使它动态,以便如果我需要增加更多的数字,那么我们可以很容易地添加标题?
from sklearn.datasets import fetch_california_housing
import pandas as pd
import seaborn as sns
import matplotlib.pyplot as plt
california_housing = fetch_cali