# additive_CNLS.R

library(readr)
library(Benchmarking)


df <- read_csv("https://raw.githubusercontent.com/ds2010/pyStoNED/master/pystoned/data/electricityFirms.csv")

y <- df$Energy
x <- cbind(df$OPEX, df$CAPEX)

sol <- stoned(x, y, RTS = "vrs", COST = 0, MULT = 0, METHOD = "MM")

sol$solutionNorm

sol$residuals