文档建议反馈控制台
首页
学习
活动
专区
工具
TVP
最新优惠活动
文章/答案/技术大牛
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

如何将伽马分布拟合到R中的数据？

伽马分布是一种常见的概率分布，用于描述正偏斜的连续随机变量。在R中，可以使用fitdistr函数来拟合伽马分布到数据。

以下是将伽马分布拟合到R中数据的步骤：

  1. 导入所需的包：首先，确保安装了MASS包，该包提供了fitdistr函数。如果未安装，可以使用以下命令进行安装：
代码语言：txt
复制
install.packages("MASS")
  1. 准备数据：将要拟合的数据存储在一个向量或数据框中。
代码语言：txt
复制
data <- c(1.2, 2.3, 3.4, 4.5, 5.6, 6.7)  # 示例数据
  1. 拟合伽马分布：使用fitdistr函数来拟合伽马分布到数据。该函数需要指定要拟合的数据和要拟合的分布类型。
代码语言：txt
复制
library(MASS)
fit <- fitdistr(data, "gamma")
  1. 查看拟合结果：可以使用summary函数来查看拟合结果，包括估计的参数值和标准误差。
代码语言：txt
复制
summary(fit)
  1. 绘制拟合曲线：可以使用curve函数来绘制拟合的伽马分布曲线。
代码语言：txt
复制
curve(dgamma(x, shape = fit$estimate[1], rate = fit$estimate[2]), from = 0, to = 10, col = "red", lwd = 2)

在拟合伽马分布时，可以根据实际情况调整数据和分布类型。如果数据不符合伽马分布，可以尝试其他分布类型或进行数据转换。

腾讯云相关产品和产品介绍链接地址：

  • 腾讯云：https://cloud.tencent.com/
  • 云服务器（CVM）：https://cloud.tencent.com/product/cvm
  • 云数据库 MySQL 版：https://cloud.tencent.com/product/cdb_mysql
  • 人工智能平台（AI Lab）：https://cloud.tencent.com/product/ailab
  • 物联网开发平台（IoT Explorer）：https://cloud.tencent.com/product/explorer
  • 移动应用开发平台（MPS）：https://cloud.tencent.com/product/mps
  • 对象存储（COS）：https://cloud.tencent.com/product/cos
  • 区块链服务（TBC）：https://cloud.tencent.com/product/tbc
  • 腾讯元宇宙：https://cloud.tencent.com/solution/metaverse
相关搜索:将伽马分布和正态分布的混合拟合到R中的数据GLM伽马分布的SAS和R之间的标准误差差异将双模/双高斯分布拟合到R中的数据集将边缘分布拟合到R中的直方图的示例如何将beta分布拟合到具有给定均值的数据？如何将von mises分布拟合到我的数据中以生成随机样本如何将月份和年份组合到R中的单个列中如何将对数回归拟合到R中的“负指数”散点图从R中的遗传数据模拟正态分布数据如何使用R中的sd将曲线拟合到数据？如何将衰减的正弦波拟合到数据中？将R中的数据下采样到给定分布如何将验证数据整合到Keras的datagen.flow中？如何将数据点拟合到Python中的函数e**x？如何将部分余弦曲线拟合到Python中的数据？如何将Pandas DataFrame中的数据与多索引组合到列表中将多个样本的回归汇总输出组合到R中的单个数据帧中如何将t分布与R中的glm()一起使用？到达R中多个Postgres服务器的分片数据分布如何在R中绘制多个数据集之间的数据方差分布？
相关搜索:
页面内容是否对你有帮助？
有帮助
没帮助

相关·内容

Python用 PyMC3 贝叶斯推理案例研究：抛硬币和保险索赔发生结果可视化

应用贝叶定理从观察到的样本数据中推导出后验参数值。 重复步骤 1-4，以获取更多数据样本。 使用 PyMC3，我们现在可以简化和压缩这些步骤。 首先，我们设定先验信念和先验β-二项分布。...泊松分布由下式给出： 其中 lambda λ 是事件的“速率”，由事件总数 （k） 除以数据中的单位数 （n） 给出 （λ = k/n）。...数据形状的分布，但是伽马泊松最适合： 泊松可以取任何正数到无穷大（0，∞），而β或均匀是[0-100]。...伽马和泊松属于同一分布家族。 伽马的峰值接近于零。 伽马尾巴走向无穷大。 伽马泊松先验为： 其中 a 是伽马形状，b 是伽马速率参数。...伽马密度函数为： 其中 a>0 是形状参数，b>0 是速率参数，以及 和 注意在 scipy 中，伽马分布使用形状 a 和尺度参数化，其中速率 b 等于尺度的倒数（速率 = 1/尺度）。

23820

Python用 PyMC3 贝叶斯推理案例研究：抛硬币和保险索赔发生结果可视化

应用贝叶定理从观察到的样本数据中推导出后验参数值。 重复步骤 1-4，以获取更多数据样本。 使用 PyMC3，我们现在可以简化和压缩这些步骤。 首先，我们设定先验信念和先验β-二项分布。...泊松分布由下式给出： 其中 lambda λ 是事件的“速率”，由事件总数 （k） 除以数据中的单位数 （n） 给出 （λ = k/n）。...数据形状的分布，但是伽马泊松最适合： 泊松可以取任何正数到无穷大（0，∞），而β或均匀是[0-100]。...伽马和泊松属于同一分布家族。 伽马的峰值接近于零。 伽马尾巴走向无穷大。 伽马泊松先验为： 其中 a 是伽马形状，b 是伽马速率参数。...伽马密度函数为： 其中 a>0 是形状参数，b>0 是速率参数，以及 和 注意在 scipy 中，伽马分布使用形状 a 和尺度参数化，其中速率 b 等于尺度的倒数（速率 = 1/尺度）。

18430

    • 图像的表示（2）：YCbCr 怎么来的？必看这篇颜色空间发展简史丨音视频基础

    用于 NTSC 制式的 YIQ 可以与 RGB 之间进行转换（这里的 R、G、B 是经过伽马校正的，表示为 R’、G’、B’）。...2.2 次幂，即在上述公式中 γ 为 2.2，这个称为显示伽马（display gamma）。...这个处理过程通常是在图像采集设备的电路中完成的，比如对于 TV 摄像机，需要将感知到的亮度 Y 通过以典型的数值 1/γ = 0.45 做逆伽马重新映射： image.png 这里伽马校正的非线性转换的过程除了解决显示伽马的问题外...因为我们的视觉系统对相对亮度差别是敏感的，如下图所示，经过伽马校正后的非线性梯度明显对人眼感知来说更均匀： 在彩色电视发明后，R、G、B 信号会被分别做完伽马校正后再合起来编码（这也就是我们前文中表示的...到如今，尽管在传输系统中我们不再有模拟噪声，但信号压缩时仍然需要量化，因此在传感数据上做伽马校正仍然是有用的。

    1.7K11

    R语言非线性回归和广义线性模型：泊松、伽马、逻辑回归、Beta回归分析机动车事故、小鼠感染、蛤蜊数据、补剂钠摄入数据|数据分享

    这使得GLM成为处理非正态数据和非线性关系的强大工具。 泊松回归和伽马回归 - 探索联系 如果我们查看火车与机动车碰撞数据（查看文末了解数据免费获取方式），我们会发现一个有趣的模式。...部分原因是这里的响应变量在残差中不是正态分布的，而是泊松分布，因为它是计数数据。 泊松回归 具有泊松误差的广义线性模型通常具有对数链接，尽管也可以具有恒等链接。...summary(clam_gamma) 我们可以重新参数化伽马分布，使得均值=形状/速率。在这种情况下，我们使用该均值和形状参数化伽马分布。离散参数是1/形状。...但是，为了更容易理解，伽马的方差随均值的平方成比例地扩展。离散参数越大，方差扩展得越快。 最后，我们可以使用纳吉尔克计的伪R2来计算R2。...# fit r2(clam_gamma) 这是正态的吗？ 你可能会问为什么这里使用伽马分布而不是正态分布？我们可以用正态误差和对数链接进行glm拟合。

    84920

    概率论基础 - 15 - 伽马分布

    本文记录伽马分布。 整数次数的伽马分布 若事件服从泊松分布，泊松分布参数为\lambda，则事件第i 次发生和第i+k 次发生的时间间隔t的分布为伽玛分布。...更一般的伽马分布 事实上，若随机变量 X 服从伽马分布，则其概率密度函数为： p(X ; \alpha, \beta)=\frac{\beta^{\alpha}}{\Gamma(\alpha)} X...整数次数伽马分布的理解 已知Gamma分布的密度函数为: f(x, \alpha, \lambda)=\frac{\lambda^{\alpha} x^{\alpha-1} e^{-\lambda x...也可以反过来说，伽马分布是n个独立的指数分布随机变量的和。...伽马函数 伽玛函数（Gamma Function）作为阶乘的延拓，是定义在复数范围内的亚纯函数，通常写成 \Gamma(x) 。 在x取值为正整数时与阶乘是统一的。

    7.2K20

    【收藏】万字解析Scipy的使用技巧！

    ，泊松分布，伽马分布 二项分布 泊松分布 伽马分布 学生分布（t-分布）和t检验 卡方分布和卡方检验 数值积分 球的体积 解常微分方程 ode类 常数和特殊函数 物理常量 from scipy import...例如伽马分布可用于描述等待K个独立随机事件发生所需要的时间，k就是伽马分布的形状参数 print(stats.gamma.stats(1)) print(stats.gamma.stats(2.0))...伽马分布的尺度参数 和随机事件发生的频率有关，由scale参数指定： stats.gamma.stats(2.0,scale=2.0) (array(4.), array(8.))...二项分布足够大时，将会无限接近泊松分布 伽马分布 观察相邻两个事件之间的时间间隔的分布情况，或者隔k个时间的时间间隔的分布情况，根据概率论，事件之间的间隔应该符合伽马分布，由于时间间隔可以是任意数值的，...因此伽马分布是连续分布。

    4.1K20

    卡方分布、方差分析

    大家好，又见面了，我是你们的朋友全栈君。 卡方分布： 首先我们先把现代数学中的数理统计中的卡方分布已经烂大街的定义先放下来，我先回到卡方检验的诞生的之地。...（说实话每当我去了解数学史时由衷的钦佩这些数学家）于是得出了大名鼎鼎的伽马函数。...，每个字段的“理论次数”（或期望次数）为： 我们之前在文章中是提出了一下两个公式的 所以（参考维基百科上如下得出了一个卡方的统计值） 自由度=(r-1)(c-1) 那我们有了卡方分布的概率密度曲线可以用来假设检验了...) plt.show() # def gammaIntegral(s): # ”’ # 计算伽马函数积分值 # :param s: 伽马含参量变量 # :param integralCeiling: 积分上限...(gama, 0, sy.inf))#返回积分值 返回误差 # # #第二种计算伽马函数的办法 # # print(mellin_transform(exp(-x), x, s)) # #第三种直接调用伽马函数库

    1.5K31

    Scipy使用简介

    ，泊松分布，伽马分布 二项分布 泊松分布 伽马分布 学生分布（t-分布）和t检验 卡方分布和卡方检验 数值积分 球的体积 解常微分方程 ode类 常数和特殊函数 物理常量 from scipy import...例如伽马分布可用于描述等待K个独立随机事件发生所需要的时间，k就是伽马分布的形状参数 print(stats.gamma.stats(1)) print(stats.gamma.stats(2.0))...伽马分布的尺度参数和随机事件发生的频率有关，由scale参数指定： stats.gamma.stats(2.0,scale=2.0) (array(4.), array(8.))...二项分布足够大时，将会无限接近泊松分布 伽马分布 观察相邻两个事件之间的时间间隔的分布情况，或者隔k个时间的时间间隔的分布情况，根据概率论，事件之间的间隔应该符合伽马分布，由于时间间隔可以是任意数值的，...因此伽马分布是连续分布。

    2.1K20

    gamma correction什么意思_伽马校正计算方法

    在标准显示器上面，如果没有伽马，数码相机拍摄到的阴影内容便会跟我们实际看到的有所差异。平时我们所说的伽马校正、伽马编码、伽马压缩，都是伽马曲线的各种应用场景，属于相似的概念。...从图中可以看到，使用伽马编码后，整个亮度范围内的灰阶基本呈现均匀分布，从而确保了后续的图像编辑，使得颜色和直方图都能够基于自然且感知均匀的色调。...图像文件的伽马 精准的图像伽马通常会采用配置文件的形式嵌入在图像文件中。...相关知识内容 动态范围：伽马编码，一方面可以确保图像数据的有效使用；另一方面，对于给定位深，它可以增加记录的动态范围。伽马还可以通过改变图像对比度，帮助显示设备或打印机管理其有限的动态范围。...伽马校正：将伽马应用于抵消某些前期操作的效果。 伽马压缩/扩展：分别指应用的伽马小于或者大于1.0的情况。因此，文件伽马可被视为伽马压缩，显示伽马可被视为伽马扩展。

    88510

    MIT：自动驾驶系统也能拨开迷雾，看清物体

    而基于可见光的系统在能见度偏低的驾驶条件中，会受到严重的局限—— 晴朗的天气里，光线从射出到返回的时间可以准确反映物体的距离；但在雾中，传感器收到的光线很可能是经水滴反射而来，不一定是从汽车需要避让的障碍物身上返回...△ 伽马分布，OT=optical thickness MIT团队利用统计学来解决这一问题。雾中水滴反射光线的形态取决于雾有多浓，平均来说，光线进入浓雾的程度要低于其进入薄雾的程度。...不过，研究人员证明了，不论是多么重的雾，反射光线到达传感器所需的时间都与伽马分布相贴合。 与钟形的高斯分布相比，伽马分布要复杂一些，呈现的形状更为多样，不对称的情况很多。...但与高斯分布相似的是，伽马分布同样可以用两个变量就完全表示出来。MIT团队通过估算这些变量得出分布，用以将被雾反射的光线过滤出来。这样一来，物体距离测定的准确度，便不会受到大雾天气的过度影响。 ?...△不同浓度雾中的人形成像 关键的一点是，MIT系统会对传感器的1,024枚像素做出1,024个伽马分布。系统可以在不同浓度的雾中稳定发挥的原因，便是每一枚像素看到的并不完全是同一片雾。 ?

    36240

    ​常用的连续概率分布汇总

    这是伽马分布的一个特殊情况。它是几何分布的连续模拟，它具有无记忆的关键性质。除了用于分析泊松过程外，还可以在其他各种环境中找到。...伽玛分布 伽玛分布（Gamma Distribution）是统计学的一种连续概率函数，是概率统计中一种非常重要的分布。“指数分布”和“χ2分布”都是伽马分布的特例。...泊松过程中连续出现之间的时间具有指数分布。 对时间序列进行建模预测接下来发生 n 个事件时就会出现伽马分布。...它在机器学习中被当作“共轭先验”使用 Gamma 函数 当形状参数α=1时，伽马分布就是参数为γ的指数分布，X~Exp（γ） 当α=n/2，β=1/2时，伽马分布就是自由度为n的卡方分布，X^...由于它可以利用概率值很容易地推断出它的分布参数，被广泛应用于各种寿命试验的数据处理。 作者：Junlin Liu

    1.7K30

    基于R语言的lmer混合线性回归模型

    如果您有一个变量将您的数据样本描述为您可能收集的数据的子集，则应该使用混合模型而不是简单的线性模型。 什么概率分布最适合数据？ 假设你已经决定要运行混合模型。...接下来你要做的是找到最适合你的数据的概率分布。 ?...＃lnorm表示对数正态 qqp （recog $ Aggression.t，“lnorm” ） ＃qqp要求估计负二项式，泊松＃和伽玛分布的参数。 可以使用fitdistr ＃函数生成估计值。 ?...如何将混合模型拟合到数据 数据是正态分布的 如果你的数据是正态分布的， 你可以使用线性混合模型（LMM）。您将需要加载lme4软件包并调用lmer函数。...如果你的数据不正态分布 用于估计模型中效应大小的REML和最大似然方法会对数据不适用正态性假设，因此您必须使用不同的方法进行参数估计。

    4.2K30

    深度学习必须掌握的 13 种概率分布

    在贝叶斯概率论中，如果后验分布 p（θx）与先验概率分布 p（θ）在同一概率分布族中，则先验和后验称为共轭分布，先验称为似然函数的共轭先验。...它的形式与伯努利分布的负对数相同。 3. 二项分布（离散） 参数为 n 和 p 的二项分布是一系列 n 个独立实验中成功次数的离散概率分布。...: r = min(r, n-r) numer = reduce(op.mul, range(n, n-r, -1), 1) denom = reduce(op.mul, range...多伯努利分布，分类分布（离散） 多伯努利称为分类分布。 交叉熵和采取负对数的多伯努利分布具有相同的形式。 5. 多项式分布（离散） 多项式分布与分类分布的关系与伯努尔分布与二项分布的关系相同。...伽马分布（连续） 如果 gamma(a,1)/gamma(a,1)+gamma(b,1)与beta(a,b)相同，则 gamma分布为β分布。 指数分布和卡方分布是伽马分布的特例。 9.

    24920

    使用Numpy和Opencv完成图像的基本数据分析（Part III）

    伽马校正| Gamma Correction 伽马校正，或通常简称为伽玛，是用于对视频或静止图像系统中的亮度或三刺激值进行编码和解码的非线性操作，伽玛校正也称为幂律变换。...然后，通过应用以下等式获得伽马校正后的输出图像： Vo = Vi ^(1 / G) 其中Vi是我们的输入图像，G是设置的伽玛值，然后将输出图像Vo缩放回0-255范围。...对于伽马值而言，G 1被称为解码伽马，并且膨胀幂律非线性的应用被称为伽马展开。Gamma值大于1将使图像显得更亮。...伽马校正的原因|Reason for Gamma Correction 我们应用伽马校正的原因是，由于我们的眼睛感知颜色和亮度这一过程与数码相机中的传感器的工作原理不同。

    77620

    数据挖掘学习小组之（概率分布）

    指数分布 在概率理论和统计学中，指数分布（也称为负指数分布）是描述泊松过程中的事件之间的时间的概率分布，即事件以恒定平均速率连续且独立地发生的过程。这是伽马分布的一个特殊情况。...伽玛分布 伽玛分布（Gamma Distribution）是统计学的一种连续概率函数，是概率统计中一种非常重要的分布。“指数分布”和“χ2分布”都是伽马分布的特例。...偏态分布 偏态分布是与“正态分布”相对，分布曲线左右不对称的数据次数分布，是连续随机变量概率分布的一种。可以通过峰度和偏度的计算，衡量偏态的程度。...威布尔分布在可靠性工程中被广泛应用，尤其适用于机电类产品的磨损累计失效的分布形式。由于它可以利用概率值很容易地推断出它的分布参数，被广泛应用于各种寿命试验的数据处理。...F分布 F分布是1924年英国统计学家R.A.Fisher提出，并以其姓氏的第一个字母命名的。它是一种非对称分布，有两个自由度，且位置不可互换。

    70610

    Matlab系列记录之图像处理（结束篇）

    ','bmp'); 结果 从结果中的1、2、3，可以看到图像按矩阵形式导入到工作区，然后显示再窗口中，并以新的格式写入当前的运行路径~ 图像运算 直方图 图像灰度分布情况的信息对图像处理来说，而直方图就是直观的一个方法...灰度非线性变换 灰度非线性变换则是使用非线性函数进行变换来实现增强图像俩高度的目的，如：对数变换和伽马变换。...： 其中，c是尺度比例常数，y是伽马常数。...；伽马变换，则是先取了原图中灰度值为51~204的数据，然后再做的变换，对比度也相对原图强了一些。...图像变换 图像变换有种把原图像的在对应维度，一个起始坐标上的图像数据往维度中的其他坐标位置重新放置，类似在运动的感觉，这之中将用到插值的方法，对变换后的整个空间坐标中的新的图像数据值进行估计，Matlab

    1.6K20

    【opencv实践】图像增强基本操作

    直方图均衡化 直方图均衡化是常见的一种图像增强技术，直方图均衡可以让像素值由狭小区域扩大到整个像素区域，如下图，横轴为像素灰度值，竖轴为该像素值在图片中的比例： 左图在直方图均衡化之后，像素值有明显的扩散分布...而本例中就是因为图像太黑，也就是像素值集中分布在黑色部分，均衡化可以提高对比度，我们可以看下均衡化的效果： //直方图均衡化equalizeHist(img_input, img_output); 可以看到效果并不是很好...对像素进行log变换和伽马变换 我们还可以尝试对图像进行变换，比如log函数变换和伽马变换。以伽马变换为例： 先看一下伽马变换的公式： 输入r就是图片某点像素值，输出s为变换后的像素值。...} normalize(img_Gamma, img_Gamma, 0, 255, CV_MINMAX); //将像素值映射到0-255范围 return img_Gamma; } //伽马变换...img_output = gammaTrans(img_input,0.4,3); 伽马变换只是方式之一，可以用来像素映射的函数多种多样，一般都视实际情况而定。

    66620

    【R语言进行数据挖掘】回归分析

    10 11 12 -0.66666667 0.44583333 0.37916667 0.41250000 -0.05416667 除了将数据代入建立的预测模型公式中...上图中红色的三角形就是预测值。 2、Logistic回归 Logistic回归是通过将数据拟合到一条线上并根据简历的曲线模型预测事件发生的概率。...model, GLM)是简单最小二乘回归（OLS)的扩展，响应变量（即模型的因变量）可以是正整数或分类数据，其分布为某指数分布族。...这个建立模型的分布参数包括binomaial（两项分布）、gaussian（正态分布）、gamma（伽马分布）、poisson(泊松分布)等。...4、非线性回归 如果说线性模型是拟合拟合一条最靠近数据点的直线，那么非线性模型就是通过数据拟合一条曲线。在R中可以使用函数nls()建立一个非线性回归模型，具体的使用方法可以通过输入'?

    1.1K30

    图像伽马校正_自动梯形校正

    3.Gamma校正 定义 伽马是显示器电光传递函数的一种，是指对线性三色值和非线性视频信号之间进行编码和解码的操作。...编码 + 显示器显示 = 结果 左图为存在硬盘中，将捕获到的物理数据做一次gamma值约为0.4的映射 中间为显示图像时，需要为每一个像素做一次gamma值约为2.2的校正，来使的最终结果为正确的物理数据...伽马2.2的存在和显示器本身没有任何因果关系，伽马2.2是为了让8位图 有限的存储空间可以表达适合人眼的色彩感知特性而不产生色彩断层（因为人眼对亮部和暗部的感知是非线性的，对暗部更敏感），伽马2.2的目的是将有限的色彩信息存储空间更多地留给暗部...均匀记录美术灰阶的样本分布 就可以解决这个问题 所以为了用有限的灰阶采集到人眼感兴趣的信息，我们需要使用伽马校正在美术灰阶和物理灰阶之间转换 三、应用 unity中的伽马设置 Unity...Unity Gamma校正 转为线性空间 – 知乎 【技术美术百人计划】图形 2.6 伽马校正_哔哩哔哩_bilibili Gamma校正与线性工作流入门讲解_哔哩哔哩_bilibili 色彩校正中的

    1.4K20

    音视频知识图谱 2022.11

    VBO 和 EBO 的作用是在 GPU 显存中开辟一块存储空间来缓存顶点数据或者图元索引数据，避免每次绘制时 CPU 内存到 GPU 显存的数据拷贝，从而提升渲染性能。...编码伽马（Encoding Gamma） 修正显示伽马过程增加的伽马则叫做『编码伽马』。 增加编码伽马通常是在图像采集设备的电路中完成的。...端到端伽马（End-to-End Gamma） 编码伽马和显示伽马的乘积就是整个图像系统的『端到端伽马』。 如果端到端伽马乘积为 1，那么显示出来的图像亮度与捕捉到的真实场景的亮度就是成线性比例的。...这样一来，我们在采集电路中采集到光信号向电信号转换时，通常会将其转换为非线性信号，以利于我们做编码，因此在传感数据上做伽马校正仍然是有用的。...HLG 曲线最早在 ARIB STD-B67 中进行了标准化，后面也进入了 ITU-R BT.2100。

    93930
    点击加载更多

    扫码

    添加站长 进交流群

    领取专属 10元无门槛券

    手把手带您无忧上云

    扫码加入开发者社群

    相关资讯

    热门标签

    更多标签

    活动推荐

      运营活动

      活动名称
      广告关闭

      腾讯云开发者

      扫码关注腾讯云开发者

      扫码关注腾讯云开发者

      领取腾讯云代金券

      热门产品

      热门推荐

      更多推荐

      Copyright © 2013 - 2024 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有

      深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 深公网安备号 44030502008569

      腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号 | 京公网安备号11010802020287

      领券