Confesión: desarrollamos software en R y Excel
Ejemplos del ecosistema basado en R que tenemos en el trabajo
07/11/2024
Contexto
- Empresa A quiere cuantificar el impacto de su publicidad en su negocio.
- Empresa A contrata a Ebiquity.
- Los consultores en Ebiquity usan estadística para medir ese impacto.
Más contexto
Los consultores de Ebiquity:
- Saben Excel
- Saben de estadística
- Saben de publicidad
Pero no tienen por qué saber programar
Paula y yo los ayudamos cuando programar es imprescindible
Por ejemplo:
- Datos en formatos complicados o datos no pequeños.
- Automatización de tareas con herramientas internas.
Sí, usamos Excel
Nuestra herramienta de modelos econométricos se llama XSec. 1
Y también VBA para llamar a R
Function RLocation() As String
Dim RLoc As String
RLoc = gwsRSpec.Range("location_RProg").Value & "\bin\x64\Rscript.exe"
If Not FileExists(RLoc) Then
RLoc = gwsRSpec.Range("location_RProg").Value & "\bin\Rscript.exe"
End If
RLocation = RLoc
End Function
La combinación VBA y R hace que ajustar un modelo sea lo más parecido a pulsar un botón
¿Qué os vamos a contar?
- El ecosistema R que sustenta todo esto
- El papel que juega Shiny
- Qué hacemos mal (sabiéndolo)
- Qué más nos gustaría hacer
Algunas de las transformaciones de nuestros modelos son muy típicas en el sector
Realmente lo hacemos en Shiny
¿Has usado Shiny alguna vez?
Un Shiny sirve para explorar y diagnosticar residuos
Otro Shiny sirve para explorar los datos de inversión en medios publicitarios
Otro lanza predicciones
¿Sabías que puedes hacer Shinys como si fueran librerías?
Terminar antes es bueno a veces. Rcpp te ayuda
//' Apply a dim ret transformation to a variable.
//'
//' @param x A numeric vector, the variable to be transformed.
//' @param type A string, the type of transformation. `"adbug", "power"`.
//' @param a The diminishing return parameter.
//' @param e The exponent of the transformation.
//' @export
// [[Rcpp::export]]
NumericVector transform_variable(NumericVector x,
Rcpp::String type,
double a,
double e)
{
if (type == "adbug") {
return pow(x, e) / (pow(x, e)+a);
} else if (type == "power") {
return pow(x,e);
} else if (type == "log2") {
for (int i = 0; i < x.size(); i++) {
// If the element is zero or less than zero, set it to one
if (x[i] <= 0) {
x[i] = 1;
}
}
return log(x);
} else if (type == "log") {
// zero and negative values will be replaced by a log average.
// The logic is the same as the one in ebimodel. The code is from ChatGPT
double sumLog = 0;
int countPos = 0;
for (int i = 0; i < x.size(); i++) {
if (x[i] > 0) {
sumLog += x[i];
countPos++;
}
}
double avgLog = countPos > 0 ? log(sumLog / countPos) : 0; // Calculate the average of the logarithms of positive numbers
NumericVector transformed(x.size());
for (int i = 0; i < x.size(); i++) {
if (x[i] > 0) {
transformed[i] = log(x[i]);
} else {
transformed[i] = avgLog;
}
}
return transformed;
} else if ((type == "pnegativo") | (type == "pneg")) {
int nTimeSteps;
nTimeSteps = x.size();
NumericVector out(nTimeSteps);
// Loop through each element in adst
for (int i = 0; i < nTimeSteps; i++) {
// Check if the element is non-negative
if (x[i] >= 0) {
// Raise the element to the given exponent
out[i] = pow(x[i], e);
}
else {
// Raise the absolute value of the element to the given exponent
// and negate it
out[i] = -1 * pow(-1 * x[i], e);
}
}
return out;
} else if ((type == "anegativo") | (type == "aneg")) {
int nTimeSteps;
nTimeSteps = x.size();
NumericVector out(nTimeSteps);
for (int i=0; i<nTimeSteps; i++){
if (x[i] > 0)
out[i] = pow(x[i], e) / (pow(x[i], e)+a);
else if (x[i] < 0){
out[i] = pow(-1.0 * x[i], e) / (pow(-1.0 * x[i], e)+a);
out[i] = -1.0 * out[i];
} else {
out[i] = x[i];
}
}
return out;
} else {
return x;
}
}
bslib hace que tus Shinys no sean tan feos como son por defecto
Cosas que (sabemos que) hacemos mal
testthat- Shiny Modules
importFromen lugar deimport
No deberíamos importar librerías enteras, sino solo las funciones que necesitamos
Para que no pasen cosas como esta y se nos llene la consola de warnings:
Quarto
A disfrutar de Sevilla
