R:如何获得与密度图的最后一个峰值相关的值？
获得与密度图的最后一个峰值相关的值，可以通过以下步骤实现：
- 密度图是一种图形表示方式，用于展示数据分布的相对密度。首先，需要明确密度图的定义和特点。
- 密度图通常通过核密度估计（Kernel Density Estimation，KDE）方法生成。核密度估计是一种非参数统计方法，用于估计概率密度函数。它基于一组观测值，在数据点周围放置核函数并进行平滑处理，得到连续的概率密度估计。
- 密度图的峰值表示数据中密度最高的区域。通常情况下，峰值与概率密度函数的极值相关联。要获得与密度图的最后一个峰值相关的值，可以按照以下步骤操作：
- a. 根据核密度估计方法计算得到密度图。
- b. 通过寻找概率密度函数的极值点，找到峰值所在的位置。
- c. 确定最后一个峰值的位置，可以采用以下方法之一：
- 寻找概率密度函数的最大值，即最后一个峰值。
- 根据峰值的高度阈值，确定峰值的位置，例如可以选择高于平均峰值高度的某个百分比。
- d. 在获得最后一个峰值的位置后，可以通过读取该位置的横坐标或纵坐标的值，获得与该峰值相关的数值。
- 在云计算领域中，密度图广泛应用于数据分析、数据挖掘、机器学习等领域。通过分析密度图的峰值，可以了解数据的分布情况、发现异常点、进行数据聚类等。
- 对于腾讯云相关产品，如果需要在云环境中进行密度图分析，可以使用腾讯云提供的弹性计算服务、大数据分析平台、人工智能服务等。具体推荐的产品包括腾讯云弹性MapReduce、腾讯云人工智能开放平台等。您可以通过访问腾讯云官方网站，获取更多关于这些产品的详细信息。
参考链接：
- 腾讯云弹性MapReduce：https://cloud.tencent.com/product/emr
- 腾讯云人工智能开放平台：https://cloud.tencent.com/product/tia
请注意，以上答案仅供参考，具体实施时需根据实际情况和需求进行调整。
相关·内容
1回答
我做了一个密度图，像这个，我会把轴线弄到很少的点，我的问题是，我需要知道，最高峰值的确切x值是多少。
我可以用循环找到它：
density = gaussian_kde(data)
aa = 0
bb = 0
for i in range(max value of data):
if density(i)[0]>aa:
aa = density(i)[0]
bb = i
在此之后，bb具有峰值的x值，但是执行这个循环的时间太长了。它现在大约需要25秒，而且在将来数据的大小将更大。
我希望这不是重
浏览 1提问于2013-07-29得票数 1
回答已采纳
有人能解释一下为什么以下两幅图会产生不同的结果吗？
prop.table(table(S))，其中'S‘是随机variable...representing轮盘轮的结果--在这里是d范数( S，mean(S)，sd(S))范围内的值列表。
下面是我的代码片段：
随机变量频率图(S)
plot(prop.table(table(S)), xlab = "Net Profit", ylab = "Probability", type = "h")
base <- seq(min(S),max(S),length = B)
pdf =
浏览 0提问于2020-04-13得票数 0
回答已采纳
1回答
我已经使用rbeta函数创建了一系列数字。
set.seed(123)
n = 100000
p1.12.2 = rbeta(n, 0.3225928, 1.2903712)
p4.7.2 = (rbeta(n, 0.3488823, 3.1399407)^2)
E2 = p4.7.2*p1.12.2
这运行得很好，但我想要找到E2的模式，所以我通过获得密度图的峰值来实现这一点。
d = density(E2)
i = which.max(d$y)
M2 = d$x[i]
M2
我一直得到模式的负值。但是beta分布被限制在0-1范围内。你知道负值是从哪里来的吗，或者有没有其他方法来获得bi
浏览 1提问于2015-09-02得票数 1
1回答
我正在将我的数据框的一列绘制为历史图。当我绘制它时，我可以看到有3个峰值，所以它是三峰。我如何注释我的图表，用它的值标记每个峰值？ 我下面的代码并排绘制了两个直方图，一个是有异常值的数据帧，另一个是没有异常值的数据帧。 fig, ax = plt.subplots(ncols=2, figsize=(10,4))
df['price'].hist(bins=40,ax=ax[0])
df_nooutlier['price'].hist(bins=40,ax=ax[1])
df_nooutlier['price'].value_counts()
浏览 24提问于2019-04-01得票数 0
请你向我解释什么是周期图，它的物理意义，以一个非常基本和简单的方式？
我在看这份文件()中的数字。它与傅里叶变换的幅度图有什么不同？
周期图和谱图之间有什么关系吗？
我能用matlab的周期图函数给出一个非常简单的例子吗？
谢谢
浏览 0提问于2012-11-26得票数 0
4回答
似乎已经有很多“密度函数中的峰值”线程了，但是我没有看到有一个专门解决这一点的。如果我错过了，很抱歉重复一下。
我的问题是:给定一个1000值的向量(样本附件)，我想识别数据的直方图或密度函数中的峰值。从下面样本数据的图像中，我可以看到直方图中的峰值，分别为0、6200和8400。但是我需要得到这些峰的确切值，最好是在一个简单的过程中，因为我有几千个这样的向量需要处理。
我最初开始使用柱状图输出本身，但无法使任何找峰命令正常工作(比如，一点也不)。我甚至不确定它将如何从splus2R包中获取splus2R命令来处理直方图对象或密度对象。这仍然是我的首选，因为我想确定每个峰值的最大频率
浏览 2提问于2012-10-30得票数 9
回答已采纳
1回答
我正在使用R语言，我必须绘制两个样本的密度比较图：
samples1：
id1 2
id2 2
id3 2
id4 2
id5 2
id6 2
id7 3
id8 3
id9 3
id10 2
和samples2：
id1 1
id2 1
id3 1
id4 1
id5 1
id6 1
id7 1
id8 2
id9 2
id10 1
第一列是ID，第二列是该ID的值。
前：
Id1 -值
Id2 -值
..。
请注意，这两个样本具有相同的大小和相同的ID，但值不同。
我使用了以下代码来绘制此图像：
library(ggplot2)
libra
浏览 0提问于2017-05-29得票数 0
回答已采纳
1回答
我可以使用提供的代码创建两个密度图。我希望做的是创建一个把两者结合在一起的图表。本质上，我想要覆盖黑线在情节的顶部，由条件着色。
d1 <- ggplot(df, aes(x=a)) + geom_density()
d2 <- ggdensity(biolReps, x = "a", color = "b",fill = "b"
浏览 0提问于2021-11-04得票数 0
回答已采纳
1回答
我想要查询一个区域中的所有事件(一个特定的位置+解释的就可以了)，并返回事件的位置，无论是通过图搜索还是FQL。问题:我只得到很少数量的事件返回。
严密的搜索，没有找到我想要的东西。特别是在一个区域内获取所有事件似乎是一个问题。
浏览 0提问于2013-06-01得票数 6
我正在尝试使用seaborn.displot()方法创建一条高斯曲线(没有条形图)。不幸的是，我得到的是y轴上的归一化值，而不是绝对值。如何解决此问题？
下面是我的代码：
height_mu = 165
height_sigma = 15
heights = np.random.normal(height_mu, height_sigma, size=10000)
plt.figure(figsize=(20, 5))
sns.distplot(heights, hist=False)
plt.axvline(165, color='red', label='Mean
浏览 16提问于2020-06-04得票数 1
回答已采纳
2回答
我有一个假设的问题来理解这个概念..
假设我捕获了一个采样率为8000 in的单声道语音片段，即4096字节的数据。通过大小为256的FFT输入前512个字节(16位编码)，将返回128个值，我将其转换为振幅。所以我的输出频率是
FFT BIN #1
0: 0*8000/256
1: 1*8000/256
.
.
127: 127*8000/256
到目前为止还好吧？所以现在还剩下3584个字节的未处理数据。因此，我对512字节的数据执行了另一个256大小的fft。并得到相同数量的结果..因此，对于这个，我再次有以下频率：
FFT BIN #2:
Example1:
0: 0*8000/256
浏览 3提问于2012-07-19得票数 2
回答已采纳
1回答
我读过蜂巢图和直方图图。但我不明白其中的区别。在直方图中，如果回收箱数目增加，那么它的表示法就会改变。那么，蜂温图如何才能帮助理解这些数据呢？
浏览 0提问于2020-04-04得票数 1
回答已采纳
1回答
我正在尝试使用R包ggplot2绘制和比较0和1之间的两组十进制数。当我在qplot中使用geom=" density“绘图时，我注意到密度曲线超过了1.0。我希望有一个不超过集合的值范围的数据的密度图，即所有面积保持在0到1之间。
有没有可能在0和1之间绘制密度图，而不超过1或0？如果是这样的话，我将如何实现这一点？我需要两个图的面积在0和1之间相等，这是数据的范围。这是我用来生成图的代码。
右图：qplot(precision,data = compare, fill=factor(dataset),binwidth = .05,geom="density", al
浏览 2提问于2011-06-20得票数 2
回答已采纳
2回答
从数据中获取FFT峰值
、、、
我正在使用Octave从头开始开发一个语音识别系统。我试图通过检测频率的差异来检测音素。目前，我已经读入了wav文件，将值组织成块，并将fft应用于整个数据。之后，我用plot(abs(real(fft(q))))绘制新数据，它创建了这个图：
我怎样才能得到频率值(图的峰值)？
浏览 1提问于2014-11-23得票数 6
2回答
白噪声季节性
、、
📷
上面的图是白噪音吗？我被某些地方的尖峰弄糊涂了。在对时间序列进行一阶差分后，得到了上述图。如何证明数据在差异后是否有季节性存在/缺失？
浏览 0提问于2019-04-21得票数 0
1回答
我使用一个巨大的4DNIFTI文件(x -y-z- subject；MRI数据)，由于大小我不能转换为csv文件并在R中打开。我想得到一系列覆盖密度图()，每个对象一个，想法只是可视化，在整个样本中密度分布没有太大变化。
但是，我可以提取每个主题的汇总统计信息(感兴趣变量的平均值、中值、SD、极差等)，并使用这些统计信息来创建密度图(至少对于正态分布的变量)。像这样的东西会很棒，但我不确定如何为密度图做这件事。
您的帮助我们将不胜感激。
浏览 2提问于2018-05-10得票数 2
2回答
解释R中的密度图
、、
我有一个以天为单位的年龄列表，我希望在密度图上显示它们。
我采用了两种方法-将x轴上的标签更改为年，并将数据除以365。这些方法给出了不同的密度估计：
df <- data.frame(id = 1:80000, age = rnorm(80000, 46, 5) * 365)
使用以下命令生成第一个图：
breaks <- seq(from = min(df$age), to = max(df$age), by = 10*365)
ggplot(data = df, aes(x = age)) +
geom_density(aes(y = ..density..)) +
浏览 0提问于2012-11-20得票数 2
回答已采纳
2回答
我有一个或多或少重复的信号(从一个到下一个并不完全相同，请参见左侧的图)。如果我使用自相关，我会得到一些最大值(右图)，但它不会告诉我哪里(哪个样本数)相关性高。它给了我滞后，但我丢失了有关位置的信息，即原始数据中出现信号的样本数。例如，在自相关中。图中，样本500处的第二峰值应对应于数据图中约样本750处的信号。我可以通过使用一个小窗口来实现这一点，该窗口在数据轨迹上移动并找到最大值，但这需要太多时间。在matlab中有没有更快的方法来做这件事？谢谢。
浏览 0提问于2012-11-28得票数 0
回答已采纳
我已经梳理了这里的几个问题，我似乎不能弄清楚我的密度图发生了什么。
我有一套属于不同文化的放射性碳测年。我需要通过时间显示日期的频率，但要按文化区分日期。堆叠的直方图可以完成这项工作(图1)，但通常不鼓励使用它们，所以这是不可能的，但我想要比频率图更平滑的图(图2)。
图1:直方图
图2:频率图
当我制作一个由培养物着色的密度图(图3)时，培养物在y轴上的相对分布与它们的计数有很大的不同。例如，在密度图中，蓝色的密度曲线远远高于紫色的密度曲线；然而，在直方图中，我们可以看到有更多的日期归因于紫色组。
图3:密度图
我的代码是不是做错了什么(见下文)？或者，也许我需
浏览 2提问于2018-06-26得票数 1
2回答
我试图在Excel中使用NORMDIST函数来创建钟形曲线，但是输出很奇怪。
我的平均值是0,0000583，标准差是0,0100323，所以当我把它插入函数NORMDIST(0,0000583; 0,0000583; 0,0100323; FALSE)时，我期望得到接近0,5的值，因为我使用的值与这个值的平均值应该是50%相同，但是函数给出的39,77输出显然是不正确的。
为什么是这样？
浏览 10提问于2013-07-30得票数 3
回答已采纳
我需要从R中的密度对象中选择单独的观测，这与有关，但我的数据的格式不同。
我从基因表达值"b“的矩阵开始，其中列代表10个观察，行代表1000个基因的表达水平。
b = matrix(rexp(1000),ncol=10)
接下来，我用每条观测作为一条单独的线来绘制密度图。
d = apply(b, 2, density, na.rm = TRUE)
plot(b, type="n", ylim=c(0,1))
for (i in 1:10) { lines( d[[i]] ) }
abline(h=0.5)
如何选择“b”中的列，其密度在'b,1=0.5‘
浏览 5提问于2015-02-09得票数 0
回答已采纳
1回答
我被赋予了一项任务，给出了离散的数据，如下所示
我需要将它分割成5块，由它创建的模板决定。
我不允许猜测模板，因为每个输入看起来都不一样。
我的方法是在数据中找到峰值(高于或低于零)，然后使用这种峰值模式来分割数据。以下是我得到的信息：(以上数据不适用)
上面的图是图中的峰值，因为我知道我有5个部分，15个点，我可以说每个部分都有3个点，然后切片，这是图中的第二个图。
在40个输入中，我成功地完成了其中的5个，因为我的“峰值检测”算法非常基本。
我应该使用什么峰值检测算法，它也可以找到局部极小值，并且有PHP实现/简单的psudo代码？我是这个领域的数据分析的初学者，所以我
浏览 5提问于2016-08-03得票数 0
首先，我不完全确定这是否是发布这篇文章的正确地点，因为也许它应该放在一个更注重统计的论坛上。然而，由于我计划用R来实现这一点，我想最好在这里发布它。如果我错了请道歉。
所以，我想做的是我想要模拟总共250.000次观测的数据，按照经验数据(离散数据)导出的核密度估计值分配一个连续(非整数)值，原始值从-5到+5不等。下面是我想要使用的分布图。
对我来说非常重要的是，我不应该根据离散概率来模拟新的数据，而是要模拟连续的数据，因为一个值可以说是2.89，而不是3或2。因此，新的值将根据图中描述的概率来分配。模拟数据中最常见的值将在+2附近，而在-4和+5附近的值则相当罕见。
关于R中的模拟数
浏览 2提问于2016-09-26得票数 3
回答已采纳
1回答
我有两列数据，x(角)和y(距离)，它们是我的实验结果(图1)。如何利用这些结果用matlab制作一个“热图”？Z轴是二维坐标系中xy点的密度值，如figure2下面。谢谢!
📷
📷
浏览 0提问于2017-11-13得票数 0
回答已采纳
我有一个numpy数组，它正在使用Matplotlib绘制。我的问题是，这些值是非常相似的，所以当它是图形时，可读性是不存在的。
0,0,0,0,0,0,0,0,46.29821447,49.49781571,49.83072758,50.89081787,98.49113721,98.5522082,99.29547499,99.91765345,99.93779431,99.95351796,99.98066963,99.99294867,100
注意其中一些值是如何聚集的，我的问题是，有什么方法可以在numpy数组上迭代并确定那些紧密结合的集群，然后应用一个不包括零值的放大来分隔它们？当
浏览 5提问于2015-04-17得票数 2
1回答
我有三组数据(df$A，df$B，df$C)，我想绘制它们的密度图(每个值的频率)。我知道我可以使用plot(density(df$A))来实现它。 但是，我想将三个密度图放在一个图中，并让df$B和df$C的密度图作为df$A密度图的置信区间带。也就是说，df$B和df$C的密度图最好是较浅的颜色或甚至是虚线。
浏览 22提问于2021-06-18得票数 1
回答已采纳
我试图在R中覆盖三个不同的密度图，以创建一个显示所有三条线的图(重叠)。我已经安装/加载了sm包，但是我尝试将它与我的数据一起使用，但没有结果。我简单地使用density()和绘制值创建了三个单独的数据图。我的代码如下所示：
library(sm)
set.seed(0)
x <- rnorm(100, 0, 1)
y <- rnorm(126, 0.3, 1.2)
z <- rnorm(93, -0.5, 0.7)
dx <- density(x)
dy <- density(y)
dz <- density(z)
plot(dx)
plot(dy)
浏览 3提问于2016-06-22得票数 4
回答已采纳
1回答
我正在尝试创建一个像这样的3d曲面图，链接在这里：
但问题是，我只有有限的数据可用，我只有峰位置和峰高的数据，而其余的数据都丢失了。在示例中，z-data需要25X25值625个数据点才能生成有效的曲面图。
我的数据看起来像这样：
所以我的问题是，有没有可能使用一些带有峰值位置值的多项式函数作为约束，以基于我所拥有的信息来生成Z数据？
欢迎任何讨论。任何形式的建议都是受欢迎的。
浏览 8提问于2021-11-23得票数 1
我的模型基于以下增长函数: x_t+1 = x_t +r* x_t *(1-x_t/K)--增长率r是通过随机抽签(mean=0.2，标准deviation=0.2)定义的。我在看股票的概率密度函数，在70年后运行10000次。当我计算概率密度函数下的面积时，当标准差为0.1时，它大致等于1，但对于0.2或0.3的标准差，则不等于1。我为区域A创建了一个有100个值的直方图。它在A=2上显示出一个很高的峰值，但是面积的值甚至上升到8，为什么它不等于1？
runs=10000
mean=0.2
sd=0.2
K=1
x_0=0.4
v=0.1
t=70
y=c()
x=x_0
for
浏览 4提问于2016-08-01得票数 2
1回答
有以下数据：
#demo data:
set.seed(1234)
library(tidyverse)
library(fs)
n = 100
time= c(15, 30, 60, 90, 120, 180, 240, 300)
treat = factor(c("trt1", "trt2", "tr1+2", "trt1+2+3"))
intensity = c(sample(1:400, n, replace=TRUE))
df <- expand.grid(time= time, tr
浏览 3提问于2021-03-02得票数 0
回答已采纳
我有一个图显示了我的数据中的点(图1) ? 和使用stat_density_2d生成的等高线图(图2) ? 等高线显然不能很好地表示原始数据。我已经使用相同的代码生成了与数据完美匹配的其他等高线图(图3)。 ? 我使用的代码是： SolidReg<-ggplot(RhyShp[,c(13,15)], aes(x=Solidity, y=Reg) ) +
stat_density_2d(aes(fill = ..level..), geom = "polygon") +
labs(x = "Solidity", y = "
浏览 51提问于2020-03-27得票数 0
1回答
我可以在图表上得到峰值，并在上面绘制线，但现在我需要将数据分成两部分(峰值之前)和(峰值之后)。
怎么能做到呢？
这是我用来获取峰值并在上面画线的代码：
library(data.table)
library(ggplot2)
ggplot(DS, aes(DS$loan_amount_ink)) + geom_density()
Ymax <- which.max(density( DS$loan_amount_ink)$y)
Xmax <- density(DS$loan_amount_ink)$x[Ymax]
ggplot(DS, aes(DS$
浏览 1提问于2019-04-28得票数 0
回答已采纳
3回答
我正在绘制R中F(1,49)的密度。当数值接近于零时，模拟的图似乎不符合与的理论图。
set.seed(123)
val <- rf(1000, df1=1, df2=49)
plot(density(val), yaxt="n",ylab="",xlab="Observation",
main=expression(paste("Density plot (",italic(n),"=1000, ",italic(df)[1],"=1, ",italic(df)[2],
浏览 3提问于2022-05-17得票数 0
我想对我的数据是否符合特定的分布函数进行直观的评估。为此，我使用R生成分位数-分位数(Q-Q)图。分布函数非常具体，没有出现在概率分布的标准列表中，因此我编写了自己的R函数来描述它。它在下面的代码中称为'DistFunc‘，由两个伽马函数的比率组成。
简而言之，我在代码中所做的就是从文件'DistributionEstimate.txt‘中读取数据，该文件包含两列。列1是x值，列2是y值。变量'a‘和'b’是我之前在另一个程序中使用该分布函数对数据的最小二乘拟合确定的最佳拟合参数。然后，我定义了DistFunc，并尝试使用qqmath函数绘制Q-Q图。
问题出现
浏览 5提问于2014-05-24得票数 2
我正在使用astropy中的Lomb-Scargle包。我尝试了人工数据，这是一个振幅为1的正弦函数：
from astropy.stats import LombScargle
t_sin2=np.arange(1000)*1.0
a_sin2=np.sin(t_sin2)
frequency=np.arange(0.001,0.5,0.001)
PSD_LS = LombScargle(t_sin2,a_sin2).power(frequency, normalization='psd')
plt.plot(frequency, PSD_LS)
我得到的图是：。PSD的峰值
浏览 5提问于2020-01-11得票数 0
我正在将我的一些R文件翻译成Python，我在使用Python中的Epanechnikov内核设置KDE图的带宽时遇到了一些问题。
我尝试使用函数seaborn.kdeplot()和statsmodels.nonparametric.kde.kdensity()，但始终没有得到正确的结果。
到目前为止，我还没有找到Python中带宽是如何缩放或标准化的模式。此外，如果我使用高斯内核，R和Python中的曲线图是相同的。因此，也许Epanechnikov内核的计算方式不同？
我使用以下代码在Python中创建KDE图：
import seaborn as sns
sns.kdeplot(y,ke
浏览 21提问于2019-03-27得票数 0
1回答
我有一个Python直方图。
我想将直方图的峰值标准化为1，这样就只有相对于条形图的高度是重要的。
我看到了一些方法，包括改变垃圾箱的宽度，但我不想这样做。
我也意识到我可以改变y轴的标签，但我也有另一个地块覆盖，所以yticks必须是实际值。
是否无法访问和更改每个回收站中的直方图“计数”？
谢谢。
浏览 3提问于2013-10-02得票数 4
回答已采纳
2回答
我正在尝试理解一个数值变量的偏度和峰度，以理解数据的形状。
我首先使用skewness命令计算，如下所示：
skewness(data$responsetime)
[1] 26.56731
以及峰度：
kurtosis(data$responsetime)
[1] 3723.961
偏度是正的，所以尾巴应该向右，峰度是>= 3。
现在我想用一个曲线图来确认偏度和峰度。我试着这样做：
plot(density(data$responsetime))
我得到了一个像下面这样的图，很难得出一些结论。我是R的新手，我正在尝试让这个图表更清晰，比如调整x大小或其他东西，但我找不到这样做的命令。奥
浏览 0提问于2016-10-07得票数 0
回答已采纳
我分析了大量正弦形式的时间序列。我已经成功地用R中的peakPick来计算时间序列中的峰值。然而，山谷带来了一个问题。基本上我用：
smoothpeaks = filter(ts2[,2], rep(1/20, 20), sides=2)
形成相对一致的山峰和山谷。然后我使用：
peakhits <- peakpick(matrix(smoothpeaks,ncol =1 ),
neighlim = 10,
deriv.lim = 0.04,
peak.min.s
浏览 1提问于2020-11-17得票数 0
回答已采纳
2回答
R中的快速傅立叶变换
、、、、
我有一个数据集，其中包含一只动物在12个月内每小时访问的次数。我想使用快速傅立叶变换来检查周期模式和周期性。在过去，我对此使用了Statistica；但是，我想使用R来获得谱密度与周期的关系图。在R中有没有简单的方法可以做到这一点？如果可能的话，我想确定12小时和24小时的活动峰值。
浏览 7提问于2012-12-23得票数 29
回答已采纳
我试着画一些派图，但得到了不同的结果：
这是在Android4.3(正确的变体)
在Android4.1 get next(不正确的变体)上：
我的代码如下：
@Override
protected void onDraw(Canvas canvas) {
super.onDraw(canvas);
// int saveFlag = Canvas.MATRIX_SAVE_FLAG | Canvas.CLIP_SAVE_FLAG | Canvas.HAS_ALPHA_LAYER_SAVE_FLAG
// | Canva
浏览 2提问于2014-11-29得票数 1
回答已采纳
我试图绘制逆gamma(0.001, 0.001)的密度，但是这个图只有一个点。我使用的命令是
alpha1 = 0.001;
beta1 = 0.001;
n = 1e+5;
r=1./gamrnd(alpha1,1/beta1,n,1);
[f,xi] = ksdensity(r);
plot(xi,f,'--m');
f的第一个项是实数，其余的都是NaN。xi的第一个项是实数，其余的都是Inf。
你能帮我解决这个问题吗。
非常感谢。
浏览 1提问于2015-07-19得票数 1
回答已采纳
2回答
我有一个图，它的y轴我需要显示密度而不是频率。
这是我使用的代码：
ggplot(stocks_orig, aes(x=Value)) + geom_histogram(aes(y=..density..), colour="black", fill="white", bins=20)+geom_density(aes(y=..density..),alpha=.2, fill="lightblue", size=1)+
geom_vline(aes(xintercept = -0.019), linetype = "dashed&#
浏览 7提问于2022-02-16得票数 1
回答已采纳
1回答
在r中绘制eda分析的hist +密度图时，输出图不适合于框架。
hist(Absenteeism_Data$Absenteeism.time.in.hours, col = "peachpuff", border = "black", prob = TRUE, xlab = "Absenteeism.time.in.hour", main = "Absenteeism_Data")
lines(density(Absenteeism_Data$Absenteeism.time.in.hours, na.rm = TR
浏览 0提问于2019-02-21得票数 0
回答已采纳
2回答
创建一个带有中断和密度的条形图
、、、、
我希望创建一种特定的条形图，而不是在ggplot2或R图形专家。每个条形的高度对应于下面的密度向量。
在中断向量中标识每个条形的x轴中断。(密度和中断向量的长度始终相同。)中断是连续的。我可以保证断点之间的宽度对于所有中断都是相同的(在下面的示例中，每个间隔的宽度都是.2 )。
如何在R中创建此图表，使条形的高度是密度，而条形的宽度和开始/结束由中断向量确定？
# density for each bar
density = c(0, 2.43053372991266e-05, 0, 2.56155325481663e-05, 7.85928661230842e-05,
6.65974683
浏览 0提问于2011-12-22得票数 1
回答已采纳
我在做这个问题，我必须区分给定的两个波形。我绘制了波形图，并确定了最大值/峰值。我的目标是找出第一个波形中没有的峰值，但在第二个波形中存在。对肉眼来说，识别这一点是相当简单的。但请你帮助我，至少在逻辑/概念，如果可能的话，在MATLAB本身的代码。
注：，我希望能够检测出两个波形的不同峰值点，其中我的意图是检测出与给定的2个波形的其他峰值有很大不同的第8个峰值。
--编辑--
花时间对信号进行滤波以去除不必要的峰值。
下图表示参考框架：
从图中可以明显看出，存在三个峰，即-1，2，3，这幅图被用作参考，以后计算的任何其他点，我们都可以忽略这些峰(可能不在完全相同的坐标上
浏览 1提问于2021-08-14得票数 0
2回答
我一直致力于从一个图中提取峰值，它看起来像这样：
但我注意到，对于我一直在处理的一些xcorr图，值的结果并不像预期的那样，它们通常看起来更像这样：
还有这个：
而不是像代码在第一张图中所做的那样尝试选择峰值，我将如何开始尝试选择向下的斜率暂时稳定的值(如图3所示)？
当我尝试在图2和图3所示的数据上以代码的当前状态运行代码时，没有得到任何有用的数据。
我想我需要一个if语句或类似的“查找极值点”部分，但我不确定这是否正确。到目前为止，我的函数.M代码如下所示：
[inputname, pathname] = uigetfile('*.wav', &#
浏览 4提问于2010-09-09得票数 1
回答已采纳
2回答
我找到了这个自相关代码，但我就是不知道什么是"float_vec_t“。它是一个类还是什么？(我正在尝试将其转换为我很熟悉的Action-script或JavaScript )
class LPCAnalysis{
public:
float_vec_t LPCAnalysis::autoCorrelation(const float_vec_t &x);
};
/* Calculate the (UN-normalized) autocorrelation for a frame of a
signal */
float_vec_t LPCAn
浏览 0提问于2011-02-12得票数 2
1回答
我正在尝试根据数据点的密度对数据进行聚类。
我想根据density.Like在这些区域周围绘制轮廓线，这样：
为了达到这一点，我正在尝试修改中的以下代码：
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import gaussian_kde
# Generate fake data
x = np.random.normal(size=1000)
y = x * 3 + np.random.normal(size=1000)
# Calculate the point density
浏览 3提问于2015-10-06得票数 1
1回答
检测R中的季节性
、、、、
问题:在R中使用周期图和FFT来检测每日数据中的周期模式。
问题是如何在R中对周期图进行编码，以检测数据中的月度、季度、半年、annual..etc周期模式。换句话说，我需要检测低频周期模式的存在(例如:1 year=> 2*pi/365，6个月=> 4*pi/365，等等)
可重现的例子：
library(weatherData)
w2009=getWeatherForYear("sfo",2009)
w2010=getWeatherForYear("sfo",2010)
w2011=getWeatherForYear("sfo&#
浏览 3提问于2015-03-03得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
