目录
?d/p/q/r+概率分布
?QQplot
d/p/q/r+概率分布
概率密度函数生成某概率分布的相关命令:d/p/q/r+概率分布
d:概率密度函数(density)
p:分布函数(probability)
q:分布函数的反函数(quantile)
r:产生相同分布的随机数(random)
#x=1时的概率密度值
> dnorm(1, mean = 0, sd = 1, log = FALSE)
[1] 0.2419707
#x=2时分布函数值
> pnorm(2, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)
[1] 0.9772499
#分布函数值为0.9772499时对应的x值
> qnorm(0.9772499, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)
[1] 2.000001
#随机生成两个符合标准正态分布的数据
> rnorm(2, mean = 0, sd = 1)
[1] 0.9217442 -1.4989003
#随机生成10个符合均匀分布的数据,在0-5范围内
> runif(10,0,5)
[1] 3.882177 1.583489 4.170966 4.642111 4.359600 1.614131 4.638566 4.383169 4.065598
[10] 1.279148
?QQplot
QQplot:判断两组数据是否来自同一分布
相关命令
- qqnorm(y, ylim, main = "Normal Q-Q Plot",xlab = "Theoretical Quantiles", ylab = "Sample Quantiles",plot.it = TRUE, datax = FALSE, ...)
只需要传入一组数据作为y,系统根据输入的y自动生成相同数量的符合标准正态分布的数据作为x。
- qqplot(x, y, plot.it = TRUE, xlab = deparse(substitute(x)),ylab = deparse(substitute(y)), ...)
需要输入两组数据,依次作为x与y,若分布相同应该只会看到一条直线。输入的两组数据数量不要求相同,如果不同系统会自动使用插值法进行填充。
- qqline(y, datax = FALSE, distribution = qnorm, probs = c(0.25, 0.75), qtype = 7, ...)
qqplot绘制的是散点图,而qqline绘制的是近似曲线
|