corspec(seewave)
corspec()所属R语言包:seewave
Cross-correlation between two frequency spectra
两个频谱之间的互相关
译者:生物统计家园网 机器人LoveR
描述----------Description----------
This function tests the similarity between two frequency spectra by returning their maximal correlation and the frequency shift related to it.
此功能测试频谱通过返回其最大的相关性和相关的频移之间的相似性。
用法----------Usage----------
corspec(spec1, spec2, f = NULL, plot = TRUE, plotval = TRUE,
method = "spearman", col = "black", colval = "red",
cexval = 1, fontval = 1, xlab = "Frequency (kHz)",
ylab = "Coefficient of correlation (r)", type="l",...)
参数----------Arguments----------
参数:spec1
a first data set resulting of a spectral analysis obtained with spec or meanspec (not in dB). This can be either a two-column matrix (col1 = frequency, col2 = amplitude) or a vector (amplitude).
第一spec或meanspec(不以dB为单位)的频谱分析,得到的数据。这可以是一个两列的矩阵(col1的=频率,col2上=振幅),或一个矢量(振幅)。
参数:spec2
a first data set resulting of a spectral analysis obtained with spec or meanspec (not in dB). This can be either a two-column matrix (col1 = frequency, col2 = amplitude) or a vector (amplitude).
第一spec或meanspec(不以dB为单位)的频谱分析,得到的数据。这可以是一个两列的矩阵(col1的=频率,col2上=振幅),或一个矢量(振幅)。
参数:f
sampling frequency of waves used to obtain spec1 and spec2 (in Hz). Not necessary if spec1 and/or spec2 is a two columns matrix obtained with spec or meanspec.
用于获取spec1和spec2(赫兹)的波的采样频率。没有必要的,如果spec1和/或spec2是一个两列的矩阵得到spec或meanspec。
参数:plot
logical, if TRUE plots r values against frequency shift (by default TRUE).
逻辑,如果TRUE图R值对频移(默认情况下TRUE)。
参数:plotval
logical, if TRUE adds to the plot maximum r value and frequency offset (by default TRUE).
逻辑,TRUE如果增加的图最大的r值和频率偏移(默认情况下,TRUE)。
参数:method
a character string indicating which correlation coefficient is to be computed ("pearson", "spearman", or "kendall") (see cor).
一个字符串,表示其中的相关性系数的计算(“皮尔森”,“斯皮尔曼”,或“肯德尔”)(见cor)。
参数:col
colour of r values.
r值的颜色。
参数:colval
colour of r max and frequency offset values.
色的r max和频率偏移值。
参数:cexval
character size of r max and frequency offset values.
字符尺寸的r max和频率偏移值。
参数:fontval
font of r max and frequency offset values.
字体的r max和频率偏移值。
参数:xlab
title of the frequency axis.
在频率轴上的标题。
参数:ylab
title of the r axis.
R轴的标题。
参数:type
if plot is TRUE, type of plot that should be drawn. See plot for details (by default "l" for lines).
如果plot是TRUE,应该得出的图。见plot的详细信息(默认情况下,“L”线)。
参数:...
other plot graphical parameters.
其他plot图形参数。
Details
详细信息----------Details----------
It is important not to have data in dB.<br> Successive correlations between spec1 and spec2 are computed when regularly shifting spec2 towards lower or higher frequencies.<br> The maximal correlation is obtained at a particular shift (frequency offset). This shift may be positive or negative.<br> The corresponding p value, obtained with cor.test, is plotted.<br> Inverting spec1 and spec2 may give slight different results, see examples.
重要的是不要以dB为单位的数据。<BR>连续之间的相关性spec1和spec2计算时经常转移spec2向更高或更低频率的最大相关性。参考在一个特定的移位(频率偏移)得到。这种转变可能是积极的还是消极的。<BR>相应的P值,得到的cor.test,绘制。<BR>的反相spec1和spec2可能会略有不同的结果,看到的例子。
值----------Value----------
If plot is FALSE, corspec returns a list containing four components:
如果plotFALSE,corspec返回一个列表,其中包含四个组成部分:
参数:r
a two-column matrix, the first colum corresponding to the frequency shift (frequency x-axis) and the second column corresponding to the successive r correlation values between spec1 and spec2 (correlation y-axis).
一个两列的矩阵中,第一柱的频移(频率x轴)和第二列对应于spec1和spec2(相关y-轴)之间的值的逐次r关联对应。
参数:rmax
the maximum correlation value between spec1 and spec2.
最大相关值之间的spec1和spec2。
参数:p
the p value corresponding to rmax.
p值对应的rmax。
参数:f
the frequency offset corresponding to rmax.
的频率偏移对应的rmax。
(作者)----------Author(s)----------
Jerome Sueur <a href="mailto:sueur@mnhn.fr">sueur@mnhn.fr</a>
参考文献----------References----------
参见----------See Also----------
spec, meanspec, corspec,
spec,meanspec,corspec,
实例----------Examples----------
data(tico)
# compare the two first notes spectra[比较这两个第一个音符谱]
a<-spec(tico,f=22050,wl=512,at=0.2,plot=FALSE)
c<-spec(tico,f=22050,wl=512,at=1.1,plot=FALSE)
op<-par(mfrow=c(2,1), mar=c(4.5,4,3,1))
spec(tico,f=22050,at=0.2,col="blue")
par(new=TRUE)
spec(tico,f=22050,at=1.1,col="green")
legend(x=8,y=0.5,c("Note A", "Note C"),lty=1,col=c("blue","green"),bty="o")
par(mar=c(5,4,2,1))
corspec(a,c, ylim=c(-0.25,0.8),xaxs="i",yaxs="i",las=1)
par(op)
# different correlation methods give different results...[不同的相关方法得到不同的结果...]
op<-par(mfrow=c(3,1))
corspec(a,c,xaxs="i",las=1, ylim=c(-0.25,0.8))
title("spearmann correlation (by default)")
corspec(a,c,xaxs="i",las=1,ylim=c(0,1),method="pearson")
title("pearson correlation")
corspec(a,c,xaxs="i",las=1,ylim=c(-0.23,0.5),method="kendall")
title("kendall correlation")
par(op)
# inverting x and y does not give exactly similar results[反相x和y不给完全相似的结果]
op<-par(mfrow=c(2,1),mar=c(2,4,3,1))
corspec(a,c)
corspec(c,a)
par(op)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|