Tengo varios datos los cuales algunos los modelo como una regresión no lineal (nlsLM) sigmoidal y otros como una regresión lineal (lm) cuadrática y necesito representarlas en el mismo gráfico. Para esto, necesito poner los datos en porcentaje absoluto (que cada columna sea expresada en porcentaje siendo el 100% el valor máximo de cada columna), por tener escalas distintas. Sé que esto se puede realizar fácilmente en excel, pero también me gustaría aprender a hacerlo en R. Para los datos de M, N, P y K yo los ajusto de la siguiente manera:
plotPoints(N ~ Dia, data = data)
NT1 <- nlsLM(N ~ a + ((b - a)/(1 + exp(-c * (Dia - d)))),
data = data,
start = list(a = min(data$N),
b = max(data$N),
c = 1, d = median(data$Dia)),
trace = TRUE, algorithm = "port")
pcrGOF(NT1, PRESS = FALSE)
overview(NT1)
plotfit(NT1, smooth = TRUE, xlab="Tiempo (dds)", ylab="KgN/ha", col.fit = "blue", lwd = 3)
Y para el parametro I lo realizo de la siguiente manera:
data %>%
ggplot(aes(x=Dia, y=I)) +
geom_point() +
geom_smooth(formula = y ~ x + I(x^2)-1,
method = "lm",
se=FALSE) +
labs(title="",
x="Tiempo (dds)",
y="I(m2)",
caption="") +
theme_minimal()
El dataframe es lo dejo aquí:
https://www.dropbox.com/s/g9ic9e0vcau9cvu/data.xlsx?dl=0
Se observa también la dificultad que varios de los datos son N/A, para evitar esto use la función predict
para predecir los valores en blanco, estos están en las columnas variable_predict_vals.
Muchas gracias!
plotPoints
, el errorError in .plotMaps(polys = data, xlim = xlim, ylim = ylim, projection = projection, : Invalid PolySet
. Hay algún proceso que hagas sobre los datos?library(mosaic) library(ggplot2) library(nlstools) library(minpack.lm) library(qpcR) library(broom) library(purrr) library(tidyverse) library(scales)
Puede que el gráfico de error porque cuando los gráfico individualmente borro las columnas que contienen los N/A, estos están por problemas en el laboratorio. Este problema lo solucione recién con la función predict.dev.off(profe_de_metodologia)