Mínimos cuadrados restringidos y ajuste exponencial.Estoy buscando encajar Y ~ 1-exp {- [C1 * (x / Xmax) ^ 1 + ....... CN * (x / Xmax) ^ N]}. El orden podría estar entre 6 y 12 (valor máximo de N). Pero antes de generar el ajuste, necesito estimar el coeficiente usando LSQ modificado. Miré el ajuste exponencial (analizar ajuste Y -> X) y traté de ingresar fórmulas, pero no quiero ajustar directamente la curva exponencial, sino que quiero estimar el mejor coeficiente de ajuste para algunas restricciones de límite superior debido a la naturaleza estrictamente monóticamente creciente de Y. Esencialmente lo más cercano que pude encontrar que representa mi modelo es https://www.mathworks.com/help/optim/ug/lsqlin.html . Parece un proceso de 2 pasos con la LSQ st. restricciones lineales, seguidas por el ajuste exponencial real. No estoy seguro si hay un enfoque de 2 pasos o de recta para el ajuste para un límite superior dado en el C1..CN.
Cualquier puntero sería muy útil.
Gracias,
Re: Mínimos cuadrados restringidos y ajuste exponencial.¿Puedes linealizar la relación con los registros: -Log (1 - Y) = c1 * (x / Xmax) + c2 * (x / Xmax) ^ 2 + ... + cN * (x / Xmax) ^ N?
(Esta forma podría no ser válida, según cómo contribuyan los errores a la expresión original).
Re: Mínimos cuadrados restringidos y ajuste exponencial.Lo intenté, pero los coeficientes estimados están desactivados en comparación con algunas de las referencias que tengo de la literatura. Tal vez me falten algunos límites o las restricciones no se han introducido correctamente. Lo intentaré de nuevo.
También intenté agregar el modelo personalizado a la biblioteca como se describe aquí https://www.jmp.com/support/help/14-2/create-a-formula-using-the-model-library.shtml
pero obteniendo error de E / S. Tengo derechos de administrador no estoy seguro de por qué sucede esto. ¿Hay algún modo de "edición" o "entrada de usuario" que me permita agregar mi modelo no lineal a la biblioteca?
Re: Mínimos cuadrados restringidos y ajuste exponencial.Las estimaciones estarán en el espacio de registro, no en las dimensiones originales. ¿Esa diferencia explicaría los resultados inesperados? Además, estos resultados no pretendían ser las mejores estimaciones, sino que solo sirven como valores iniciales para los parámetros antes de ajustar los datos con la plataforma No lineal. Debe volver a transformar los valores iniciales con la función Exp ().
No sé sobre el error de E / S. ¿Es posible que el archivo de script esté bloqueado?
Sugiero hacer su modelo personalizado de forma manual / interactiva primero. Es decir, definir el modelo primero como una fórmula de columna. De esta manera demostraremos que tal modelo es factible. La Biblioteca de modelos es solo una forma de ajustar un modelo personalizado. No es requerido. De esta manera, también proporcionará la expresión correcta si decide incluirla más adelante en la Biblioteca de modelos.
