2015-12-08 39 views
2

我想創建一個QQ圖來測試,如果我的數據可以由Weibull分佈使用命令R:試圖創造威布爾分佈QQ圖,得到錯誤「NAS通過Coersion推出」

建模使用數據

x =c(3.367, 0.769,0.8,1,1.2) 

qqplot(x,'weibull') 

我不斷收到帶有與錯誤

"In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion" 

,並找不出原因。這是否意味着我無法將Weibull分佈符合我的數據?如果有人能幫我指出爲什麼這不起作用,我將非常感激。

+2

是什麼讓你覺得qqplot可以這樣使用?你是否在使用'stats :: qqplot'以外的東西? – rawr

回答

0

製作qqplot需要一些分佈來比較您的數據。您需要先建立一些威布爾分佈,然後創建該圖。例如:

x =sort(c(3.367, 0.769,0.8,1,1.2)) 
dist = rweibull(5, 2, 1) 
qqplot(dist, x) 

您是否具有威布爾分佈的參數?有關更多信息,請參閱?rweibull,但似乎需要n,形狀參數和比例參數。

+0

謝謝你的幫助!我能夠計算分配的一些參數! –

0

您可以使用qualityTools包和qqPlot()函數。 如果您需要估計參數,您可以使用fitdistrplus包。

fit.weibull <- fitdist(x, "weibull") 
    qqplot(qweibull(ppoints(length(x), shape = fit.weibull$estimate[1], 
      scale = fit.weibull$estimate[2]), x)