Great Deal! Get Instant $10 FREE in Account on First Order + 10% Cashback on Every Order Order Now

indnum,year,y,k,l,int,ind 1,1947, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,11516,Farms 1,1948, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12515,Farms 1,1949, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,11448,Farms 1,1950,...

1 answer below »
indnum,year,y,k,l,int,ind
1,1947, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,11516,Farms
1,1948, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12515,Farms
1,1949, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,11448,Farms
1,1950, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12495,Farms
1,1951, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14750,Farms
1,1952, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14952,Farms
1,1953, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,13543,Farms
1,1954, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,13972,Farms
1,1955, XXXXXXXXXX,8090.13, XXXXXXXXXX,14342,Farms
1,1956, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14526,Farms
1,1957, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14893,Farms
1,1958, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,16652,Farms
1,1959, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,17729,Farms
1,1960, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,17546,Farms
1,1961, XXXXXXXXXX,7842.09, XXXXXXXXXX,18399,Farms
1,1962, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,19702,Farms
1,1963, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,20715,Farms
1,1964, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,20320,Farms
1,1965, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,21647,Farms
1,1966, XXXXXXXXXX,9312.11, XXXXXXXXXX,23803,Farms
1,1967, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,24577,Farms
1,1968, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,25021,Farms
1,1969, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,26673,Farms
1,1970, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,28097,Farms
1,1971, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,30067,Farms
1,1972, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,33724,Farms
1,1973, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,44539,Farms
1,1974, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,47744,Farms
1,1975, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,48760,Farms
1,1976, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,53102,Farms
1,1977, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,55957,Farms
1,1978, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,66865,Farms
1,1979, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,80358,Farms
1,1980, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,86426,Farms
1,1981, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,87655,Farms
1,1982, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,86565,Farms
1,1983, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,85089,Farms
1,1984, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,90727,Farms
1,1985, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,83965,Farms
1,1986, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,79296,Farms
1,1987, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,85012,Farms
1,1988, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,92295,Farms
1,1989, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,97979,Farms
1,1990, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,102106,Farms
1,1991, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,104894,Farms
1,1992, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,102719,Farms
1,1993, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,110884,Farms
1,1994, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,113902,Farms
1,1995, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,123494,Farms
1,1996, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,123174,Farms
1,1997, XXXXXXXXXX,48275.5, XXXXXXXXXX,131244,Farms
1,1998, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,128792,Farms
1,1999, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,129015,Farms
1,2000, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,128272,Farms
1,2001, XXXXXXXXXX, XXXXXXXXXX,43177,134241,Farms
1,2002, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,127864,Farms
1,2003, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,133736,Farms
1,2004, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,142914,Farms
1,2005, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,148506,Farms
1,2006, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,156408,Farms
1,2007, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,188952,Farms
1,2008, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,209751,Farms
1,2009, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,191372,Farms
1,2010, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,198522,Farms
1,2011, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,220731,Farms
1,2012, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,249436,Farms
1,2013, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,245982,Farms
1,2014, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,271154,Farms
2,1947, XXXXXXXXXX,531.245,328.398,717,"Forestry fishing and related activities"
2,1948, XXXXXXXXXX,670.891,370.417,843,"Forestry fishing and related activities"
2,1949, XXXXXXXXXX,602.628,381.694,701,"Forestry fishing and related activities"
2,1950,2161.18,746.585,417.594,997,"Forestry fishing and related activities"
2,1951, XXXXXXXXXX,853.577,465.849,1145,"Forestry fishing and related activities"
2,1952, XXXXXXXXXX,807.22,513.167,1145,"Forestry fishing and related activities"
2,1953, XXXXXXXXXX,774.451,569.046,1225,"Forestry fishing and related activities"
2,1954, XXXXXXXXXX,738.561,617.489,1035,"Forestry fishing and related activities"
2,1955, XXXXXXXXXX,874.389,623.955,1303,"Forestry fishing and related activities"
2,1956, XXXXXXXXXX,853.488,672.666,1342,"Forestry fishing and related activities"
2,1957, XXXXXXXXXX,778.904,684.279,1177,"Forestry fishing and related activities"
2,1958, XXXXXXXXXX,778.016,742.421,1077,"Forestry fishing and related activities"
2,1959, XXXXXXXXXX,900.162,776.411,1500,"Forestry fishing and related activities"
2,1960, XXXXXXXXXX,890.883,780.349,1593,"Forestry fishing and related activities"
2,1961, XXXXXXXXXX,892.698,808.25,1681,"Forestry fishing and related activities"
2,1962, XXXXXXXXXX, XXXXXXXXXX,847.269,1860,"Forestry fishing and related activities"
2,1963, XXXXXXXXXX, XXXXXXXXXX,867.701,2231,"Forestry fishing and related activities"
2,1964, XXXXXXXXXX,1325.42,926.063,2413,"Forestry fishing and related activities"
2,1965, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,2594,"Forestry fishing and related activities"
2,1966, XXXXXXXXXX,1462.64, XXXXXXXXXX,2933,"Forestry fishing and related activities"
2,1967, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,3125,"Forestry fishing and related activities"
2,1968, XXXXXXXXXX,1817.25, XXXXXXXXXX,3692,"Forestry fishing and related activities"
2,1969, XXXXXXXXXX,2021.88, XXXXXXXXXX,4217,"Forestry fishing and related activities"
2,1970, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,4146,"Forestry fishing and related activities"
2,1971, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,4882,"Forestry fishing and related activities"
2,1972, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,6779,"Forestry fishing and related activities"
2,1973, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,8023,"Forestry fishing and related activities"
2,1974, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,8764,"Forestry fishing and related activities"
2,1975, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,8002,"Forestry fishing and related activities"
2,1976, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,9874,"Forestry fishing and related activities"
2,1977, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12379,"Forestry fishing and related activities"
2,1978, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14014,"Forestry fishing and related activities"
2,1979, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14942,"Forestry fishing and related activities"
2,1980, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14138,"Forestry fishing and related activities"
2,1981, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,15417,"Forestry fishing and related activities"
2,1982, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12897,"Forestry fishing and related activities"
2,1983, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,13854,"Forestry fishing and related activities"
2,1984, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14168,"Forestry fishing and related activities"
2,1985, XXXXXXXXXX,4930.83, XXXXXXXXXX,11561,"Forestry fishing and related activities"
2,1986, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12270,"Forestry fishing and related activities"
2,1987, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,10731,"Forestry fishing and related activities"
2,1988, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12208,"Forestry fishing and related activities"
2,1989, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12239,"Forestry fishing and related activities"
2,1990, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12935,"Forestry fishing and related activities"
2,1991, XXXXXXXXXX,5719.68, XXXXXXXXXX,12520,"Forestry fishing and related activities"
2,1992, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,15621,"Forestry fishing and related activities"
2,1993, XXXXXXXXXX,5902.41, XXXXXXXXXX,17190,"Forestry fishing and related activities"
2,1994, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,18464,"Forestry fishing and related activities"
2,1995, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,18502,"Forestry fishing and related activities"
2,1996, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,18845,"Forestry fishing and related activities"
2,1997,38171, XXXXXXXXXX, XXXXXXXXXX,17752,"Forestry fishing and related activities"
2,1998, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,20973,"Forestry fishing and related activities"
2,1999, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,23426,"Forestry fishing and related activities"
2,2000, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,16744,"Forestry fishing and related activities"
2,2001, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,16965,"Forestry fishing and related activities"
2,2002, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,17673,"Forestry fishing and related activities"
2,2003, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,18342,"Forestry fishing and related activities"
2,2004, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,21545,"Forestry fishing and related activities"
2,2005, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,22556,"Forestry fishing and related activities"
2,2006, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,18745,"Forestry fishing and related activities"
2,2007, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,15991,"Forestry fishing and related activities"
2,2008, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,16264,"Forestry fishing and related activities"
2,2009, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12898,"Forestry fishing and related activities"
2,2010, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,14013,"Forestry fishing and related activities"
2,2011, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,16689,"Forestry fishing and related activities"
2,2012, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,12066,"Forestry fishing and related activities"
2,2013, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,13348,"Forestry fishing and related activities"
2,2014, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,17670,"Forestry fishing and related activities"
3,1947, XXXXXXXXXX, XXXXXXXXXX,577.479,1708,"Oil and gas extraction"
3,1948,5309.86, XXXXXXXXXX,714.396,2278,"Oil and gas extraction"
3,1949, XXXXXXXXXX, XXXXXXXXXX,726.232,1800,"Oil and gas extraction"
3,1950,5149.55, XXXXXXXXXX,438.208,1864,"Oil and gas extraction"
3,1951, XXXXXXXXXX, XXXXXXXXXX,905.155,2266,"Oil and gas extraction"
3,1952, XXXXXXXXXX, XXXXXXXXXX,743.226,2251,"Oil and gas extraction"
3,1953, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,2549,"Oil and gas extraction"
3,1954, XXXXXXXXXX, XXXXXXXXXX,938.448,2396,"Oil and gas extraction"
3,1955, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,2420,"Oil and gas extraction"
3,1956, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,2741,"Oil and gas extraction"
3,1957, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,3341,"Oil and gas extraction"
3,1958, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,3003,"Oil and gas extraction"
3,1959,8553.13, XXXXXXXXXX, XXXXXXXXXX,3252,"Oil and gas extraction"
3,1960, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,3387,"Oil and gas extraction"
3,1961, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,3480,"Oil and gas extraction"
3,1962,9224.36, XXXXXXXXXX,1134.67,3665,"Oil and gas extraction"
3,1963, XXXXXXXXXX, XXXXXXXXXX, XXXXXXXXXX,3889,"Oil and gas extraction"
3,1964, XXXXXXXXXX, XXXXXXXXXX
Answered 4 days After Nov 04, 2021

Solution

Mohd answered on Nov 09 2021
123 Votes
-
-
-
11/8/2021
li
ary(readr)
li
ary(magrittr)
li
ary(dplyr)
li
ary(ggplot2)
li
ary(rmarkdown)
li
ary(MASS)
li
ary(skimr)
li
ary(ggeffects)
li
ary(readr)
klema3 <- read_csv("data/klema3.csv")
View(klema3)
skim(klema3)
Data summary
    Name
    klema3
    Number of rows
    4420
    Number of columns
    7
    _______________________
    
    Column type frequency:
    
    characte
    1
    numeric
    6
    ________________________
    
    Group variables
    None
Variable type: characte
    skim_variable
    n_missing
    complete_rate
    min
    max
    empty
    n_unique
    whitespace
    ind
    0
    1
    5
    66
    0
    65
    0
Variable type: numeric
    skim_variable
    n_missing
    complete_rate
    mean
    sd
    p0
    p25
    p50
    p75
    p100
    hist
    indnum
    0
    1
    33.00
    18.76
    1.00
    17.00
    33.00
    49.00
    65.0
    ▇▇▇▇▇
    yea
    0
    1
    1980.50
    19.63
    1947.00
    1963.75
    1980.50
    1997.25
    2014.0
    ▇▇▇▇▇
    y
    0
    1
    140818.40
    257726.32
    230.56
    12308.21
    45605.68
    147478.41
    2820728.0
    ▇▁▁▁▁
    k
    0
    1
    33430.09
    109737.20
    20.90
    2063.38
    6984.37
    24506.90
    1939544.1
    ▇▁▁▁▁
    l
    0
    1
    44697.50
    90442.06
    14.52
    3194.63
    13133.90
    40102.09
    790195.5
    ▇▁▁▁▁
    int
    0
    1
    62690.81
    103330.22
    15.65
    5241.25
    20151.50
    71683.25
    788034.0
    ▇▁▁▁▁
klema3$log_y<-log(klema3$y)
klema3$log_l<-log(klema3$l)
klema3$log_k<-log(klema3$k)
k_mod<-lm(log_y~I(log_k)+I((log_k)^2)+I(log_l)+I((log_l)^2)+I(log_k*log_l),data=klema3)
summary(k_mod)
##
## Call:
## lm(formula = log_y ~ I(log_k) + I((log_k)^2) + I(log_l) + I((log_l)^2) +
## I(log_k * log_l), data = klema3)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.78160 -0.22818 -0.05519 0.18501 1.90069
##
## Coefficients:
## Estimate Std. E
or t value Pr(>|t|)
## (Intercept) 2.057833 0.108113 19.03 <2e-16 ***
## I(log_k) 0.402945 0.029025 13.88 <2e-16 ***
## I((log_k)^2) 0.058295 0.002982 19.55 <2e-16 ***
## I(log_l) 0.491054 0.027102 18.12 <2e-16 ***
## I((log_l)^2) 0.058983 0.003497 16.86 <2e-16 ***
## I(log_k * log_l) -0.113367 0.006044 -18.76 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard e
or: 0.3494 on 4414 degrees of freedom
## Multiple R-squared: 0.959, Adjusted R-squared: 0.959
## F-statistic: 2.067e+04 on 5 and 4414 DF, p-value: < 2.2e-16
emmeans::emmeans(k_mod,specs = c("log_l","log_k"),at=c("l","k"),type="response")
## log_l log_k emmean SE df lower.CL upper.CL
## 9.324 8.869 10.55 0.007373 4414 10.53 10.56
##
## Confidence level used: 0.95
plot(k_mod)
a2_1<-ggpredict(k_mod,c("log_l","log_k"))
a2_1
## # Predicted values of log_y
##
## # log_k = 7.07
##
## log_l | Predicted | 95% CI
## ----------------------------------
## 2.68 | 7.41 | [ 7.23, 7.59]
## 7.47 | 8.79 | [ 8.78, 8.81]
## 8.60 | 9.51 | [ 9.49, 9.53]
## 9.48 | 10.18 | [10.15, 10.21]
## 10.23 | 10.82 | [10.77, 10.87]
## 13.58 | 14.48 | [14.24, 14.72]
##
## # log_k = 8.87
##
## log_l | Predicted | 95% CI
## ----------------------------------
## 2.68 | 9.26 | [ 8.97, 9.55]
## 7.47 | 9.67 | [ 9.64, 9.69]
## 8.60 | 10.16 | [10.14, 10.17]
## 9.48 | 10.64 | [10.63, 10.66]
## 10.23 | 11.13 | [11.11, 11.15]
## 13.58 | 14.11 | [13.98, 14.24]
##
## # log_k = 10.67
##
## log_l | Predicted | 95% CI
## ----------------------------------
## 2.68 | 11.49 | [11.06, 11.93]
## 7.47 | 10.92 | [10.84, 10.99]
## 8.60 | 11.18 | [11.14, 11.21]
## 9.48 | 11.48 | [11.46, 11.51]
## 10.23 | 11.82 | [11.80, 11.84]
## 13.58 | 14.11 | [14.05, 14.18]
plot(a2_1)
lfs21 <- readRDS("~/data/lfs21.rds")
lfs<-lfs21%>%
filter(age_12!="70 and over")
lfs$age=lfs$age_12
lfs$wage=lfs$hrlyearn
lfs$age=as.numeric(lfs$age)
lfs$wage=as.numeric(lfs$wage)
lfs<-lfs%>%
...
SOLUTION.PDF

Answer To This Question Is Available To Download

Related Questions & Answers

More Questions »

Submit New Assignment

Copy and Paste Your Assignment Here