试图fitdistr()在R.中使用chi_square分布.这里的文档(对我来说不是很有用):https://stat.ethz.ch/R-manual/R-devel/library/MASS/html/fitdistr. HTML
问题1:chi_df下面有以下输出: 3.85546875 (0.07695236).第二个数字是多少？方差或标准差？
问题2:fitdistr生成由Chi-SQ分布定义的'k'.如何拟合数据,使得缩放常数为"A"？我愚蠢地使用下面的第14-17行.显然不好.
问题3:Chi-SQ分布仅定义为某个x范围吗？(方差定义为2K,而均值= k.这必须要求一些约束的x范围......统计问题不是编程......)
nnn = 1000;
## Generating a chi-sq distribution
chii <- rchisq(nnn,4, ncp = 0);
## Plotting Histogram
chi_hist <- hist(chii);
## Fitting. Gives probability density which must be scaled.
chi_df <- fitdistr(chii,"chi-squared",start=list(df=3));
chi_k <- chi_df[[1]][1];
## Plotting a fitted line:
## Spanning x-length of chi-sq data
x_chi_fit <- 1:nnn*((max(chi_hist[[1]][])-min(chi_hist[[1]][]))/nnn);
## Y data using eqn for probability function
y_chi_fit <- (1/(2^(chi_k/2)*gamma(chi_k/2)) * x_chi_fit^(chi_k/2-1) …我在我的代码的一个实例中获得了以下查询的'integer(0)'结果,但是它运行正常:
data.dat文件:
xx,线性,平方,高斯,rando,水果,颜色,类型,xxx,yyy
1,1,1,1,1.1,苹果,蓝色,金色,1,1
2,3,4,1,2.5,苹果,红,金,2,1
3,2,9,2,4.4,橙色,蓝色,银色,1,1
4,4,16,3,5.9,橙色,蓝色,金色,1,1
5,5,25,5,5.5,桃,蓝,金,1,1
6,6,32,12,6.9,桃,蓝,金,1,2
7,7,48,24,7.2,苹果,蓝色,银色,1,1
8,9,66,30,7.4,苹果,蓝色,金色,1,2
9,8,84,31,7.6,梨,红,金,1,1
10,10,102,30,1.5,橙色,红色,金色,1,1
data2 <- read.csv(file="data.dat",head=TRUE,sep=",");
which(data2$color=="red" , arr.ind=TRUE)
这也不起作用:
which(as.character(data2$color)=="red" , arr.ind=TRUE)
我觉得我正在失去理智,我已经使用了这个功能数百次没有问题......
我希望限制 optim() 的输出参数。对于我的正弦函数（其中“designL”是我的自变量，“ratio”是我的因变量数据，dfm 是我的数据帧），它不必要地收敛了数千个异相相位：
lo_0 = 2e-6
kc_0 = 80000
min.RSS <- function(data, par) {
with(data, sum( (sin(par[2] *(par[1] + designL))^2 - ratio)^2) )
}
resultt <- optim(par = c(lo_0, kc_0), min.RSS, data = dfm)
我想从 0:2e-5 限制 lo_0（相移）。我找到了一些关于此的文档，但它没有详细说明如何实现：https : //ubuntuforums.org/showthread.php?t=1420061