현재 Henningsen/Henningsen (2011)의 micEconCES 패키지를 사용하여 일부 추정을하려고합니다. 내 문제는 내가 R에 익숙하지 않고 패키지로 견적을 얻기 위해 자체 데이터 집합을 구현하려고한다는 것입니다. 논문의 저자는이 데이터 세트를 사용하여 견적을 작성했습니다.CES 생산 기능 micEconCES를 사용한 추정
R> set.seed(123)
R> cesData <- data.frame(x1 = rchisq(200, 10), x2 = rchisq(200, 10), x3 = rchisq(200, 10), x4 = rchisq(200, 10))
R> cesData$y2 <- cesCalc(xNames = c("x1", "x2"), data = cesData, + coef = c(gamma = 1, delta = 0.6, rho = 0.5, nu = 1.1))
R> cesData$y2 <- cesData$y2 + 2.5 * rnorm(200)
R> cesData$y3 <- cesCalc(xNames = c("x1", "x2", "x3"), data = cesData, coef = c(gamma = 1, delta_1 = 0.7, delta = 0.6, rho_1 = 0.3, rho = 0.5, + nu = 1.1), nested = TRUE)
R> cesData$y3 <- cesData$y3 + 1.5 * rnorm(200)
R> cesData$y4 <- cesCalc(xNames = c("x1", "x2", "x3", "x4"), data = cesData, coef = c(gamma = 1, delta_1 = 0.7, delta_2 = 0.6, delta = 0.5, rho_1 = 0.3, rho_2 = 0.4, rho = 0.5, nu = 1.1), nested = TRUE)
R> cesData$y4 <- cesData$y4 + 1.5 * rnorm(200)
첫 번째 줄에서는 난수 생성기의 "시드"를 설정하여 이러한 예제를 정확히 같은 데이터 집합으로 복제 할 수 있습니다. 두 번째 줄은 각각 200 개의 관측치가 있고 10 자유도가있는 임의의 χ2 분포에서 생성되는 4 개의 입력 변수 (x1, x2, x3 및 x4라고 함)가있는 데이터 세트를 만듭니다. 세 번째, 네 번째 및 일곱 번째 명령은 micEconCES 패키지에 포함 된 cesCalc 함수를 사용하여 2, 3 및 4 개의 입력 (각각 y2, y3 및 y4라고 함)이있는 CES 함수의 결정 론적 출력 변수를 계산합니다. CES 생산 기능이 주어진다. 이제 저는 제 논문에서 자본과 노동의 두 가지 투입 사례에 대해 미국의 CES 기능을 집계 수준에서 평가하려고합니다. 그래서 제가 한 것은 1990-2015 년 세계 은행 데이터베이스의 데이터를 수집 한 것입니다. 자본과 총 노동력을 위해 총 고정 자본 형성을 사용했습니다.
저자는 f.e. 비선형 추정 다음과 같은 방법
R> cesNls <- nls(y2 ~ gamma * (delta * x1^(-rho) + (1 - delta) * x2^(-rho))^(-phi/rho), + data = cesData, start = c(gamma = 0.5, delta = 0.5, rho = 0.25, phi = 1)) R> print(cesNls)
는 지금은 Data_Extract_From_World_Development_Indicators라고 내 자신의 데이터 세트에 대한 정확한 같은 일을합니다. 그래서 내가 무슨 짓을했는지 먼저입니다
ceslan <- cesCalc(xNames = c("GrossFixedCapitalFormation", "LaborForce"), data = Data_Extract_From_World_Development_Indicators, coef = c(gamma = 1, delta = 0.6, rho = 0.5, nu = 1.1))
R> cesData$y2 <- cesCalc(xNames = c("x1", "x2"), data = cesData, coef = c(gamma = 1, delta = 0.6, rho = 0.5, nu = 1.1))
모든
복제 및 자본과 노동에 대한 내 두 변수 X1과 X2를 대체했다. 는 이후 나는 내 문제이고 이것이 지금 비선형 추정cesulan <- nls(y2 ~ gamma * (delta * GrossFixedCapitalFormation^(-rho) + (1-delta)*LaborForce^(-rho))^(-phi/rho), data = Data_Extract_From_World_Development_Indicators, start = c(gamma = 0.5, delta = 0.5, rho = 0.25, phi = 1))
을하려고 노력 : 내 데이터 세트에서 y2 라인 수 무엇을 의미하는지 변수를 잘 모릅니다. 나는 수식에서 y2 ~ gamma * ...라는 용어를 볼 수 있습니다. 나머지 용어들에 대해 그려져 있지만, 어떤 종류의 값을 입력해야하는지 알지 못합니다. 누구라도 조언을하나요? Hennigsen & Hennisgen (2011)에 미리