code-statistics

0熱度

1回答

想從三角形分佈生成數字，它只有三個參數：在我的情況下，a，b，c和c被定義爲（a + b）/ 2。比方說，我有一個向量x： x <- c(1,-1,2,-2,3,-3,4,-4,5,-5,11,-11,12,-12,13,-13) 而且我想因爲有在向量x負數產生許多新的值。因此，我可以用三角分佈生成的數字代替負值。 library(triangle) c = abs(x[x<0])/2

0熱度

1回答

錯誤使用nlme;畸形因子

我正在使用nlme軟件包來學習多層次模型，並遵循教科書「Discovering Statistics Using R」的示例。 Mixed Models Code 的數據集是蜜月Period.dat，在他們的同伴網站也可下載。 Data Set - Multilevel Models require(nlme) require(reshape2) satisfactionData = read

0熱度

1回答

如果有幾個條件適用，如何合併SAS中的行？

我需要檢查我的數據集中的三列中的相同值，然後添加滿足該條件的兩行。輸出應該看起來像之前的行，但有一個單元格具有總計值。例子： Security date flag volume CUSIP A 20150707 1 100 xxx B 20150708 0 100 xyz B 20150708 0 200 xyz B 20150709 1 150 xyz C 2015070

0熱度

1回答

SAS中的DO循環和數據輸入

我有SAS的大學版，即它在瀏覽器窗口中打開。我有來自治療組A，B和C的數據。我嘗試使用DO循環分別處理這些組以進行比較。當數據長度相同時，我可以在一個嵌套的DO循環中執行此操作。但這些團體的觀察數量不同，我遇到了麻煩。這裏是我的代碼： data AirPoll1 (keep = Group Ozone); label Group = "Treatment Group"; label Ozon

1熱度

1回答

sql - 針對列樣本中的每個唯一值2來自另一列的不同值

我困在一個困難的sql聚合問題上。考慮下面的表/視圖： Column1 Column2 1 2564 2 6550 1 3578 2 6548 2 4789 1 9876 我想設計一個查詢，以做到以下幾點：對於每一個不同的列1值，樣品2分的記錄。採樣策略可能是某種引導/重新採樣，因爲可能沒有太多的數據點。因此該表將成爲： Column1 Column2 1

0熱度

2回答

我可以推定一個條件變量嗎？

我想要推測關於是否有人在英國出生的wave 1到wave 2的數據。我懷疑egen函數可以工作，但我不確定代碼的樣子是什麼？正如你所看到的，我需要分配在波1人ID爲1的相同born in the uk響應波2 我知道我可以通過重塑數據集寬做格式，但你知道是否有其他方式嗎？

0熱度

1回答

預測與華冠通訊（）平穩序列

我在預測鍛鍊工作。首選模型是ARIMA（0,0,1）（0,1,1）4，其中有三個外生變量（Forestalling.1，Forestalling.2，Break）。我的因變量是Pmean，平均房價，外生變量是指示立法和財產危機變化的虛擬變量（這些變量由以下值0,1，-1組成）。我最初的方法是者區分原始的和適合ARIMA（）模型;然而，這使我試圖預測系列作爲預測是在平穩序列完成後麻煩 - DIFF

0熱度

1回答

包含求和的非線性方程R

我正在嘗試在R中實現這個方程到非線性求解器中最困難的時候。我正在嘗試nleqslv和BB程序包，但到目前爲止只能得到錯誤。我已經搜索並閱讀文檔，直到我的眼睛流血，但我無法將我的大腦包裹在它周圍。公式本身是這樣的： The Equation s2 * sum(price^(2*x+2)) - s2.bar * sum(price^(2*x)) = 0 凡s2，s2.bar，和price是已知的相等長

0熱度

1回答

R矩陣中的代碼錯誤：'dimnames'必須是列表

我目前在此網站上使用蒙特卡洛方法http://www.quantpsy.org/medmc/medmc.htm。雖然代碼（有一些小的調整）與我的2x2或3x3矩陣的工作，我不斷收到以下錯誤代碼爲我的4x4矩陣：錯誤矩陣（C（0.0461705，0，0，0，0， 0.0028639，0，0，0，0，0.0740766，： 'dimnames' 必須是一個列表什麼我做錯了，我該如何解決這一錯誤消息

0熱度

1回答

在數值優化似然函數與R，達到最小，但海森矩陣不是正半定

最近，我已經構建了一個統計模型負對數似然性被最小化。有9個參數需要估算（實際上我想再增加兩個）。在R中已經使用了幾種優化方法，包括optim,GenSA, DEoptim,Solnp。然後我得到了最低的滿意。在下一過程中，以計算t值，有必要計算SE： sqrt(diag(solve(hessian))) 然而，發生由於Hessian矩陣誤差沒有半正定的，在主存在負數對角元素。我試過optim