Regresiones no lineales en Libreoffice

Muchos me dijeron que no se podía y que no había nada mejor que el Solver del Excel. Me puse manos a la obra y logré resolver una curva de selectividad sigmoide  utilizando solamente Libreoffice Calc.

La verdad que me pareció de lo mas sencillo y eficiente, y me parecía una mezquindad no compartir esto, así que les dejo un pequeño y sencillo tutorial.

Supongamos que tenemos los siguientes datos extraidos de un muestreo con sobrecopo para eficiencia de un determinado elemento de selectividad en redes de arrastre:

Talla (cm) Copo Sobrecopo Total Fracción Retenida SL Calculado
(Fracción Retenida Calculada)
Error^2
31 1 88 89 0,0112359551 1,16286056408873E-078 0,0001262467
32 2 128 130 0,0153846154 5,30778667818561E-070 0,0002366864
33 7 136 143 0,048951049 2,42269798212669E-061 0,0023962052
34 10 160 170 0,0588235294 1,10582166700926E-052 0,0034602076
35 9 114 123 0,0731707317 5,04743706499355E-044 0,005353956
36 21 95 116 0,1810344828 2,30386342437776E-035 0,0327734839
37 25 94 119 0,2100840336 1,05158055659529E-026 0,0441353012
38 22 87 109 0,2018348624 4,79985773161847E-018 0,0407373117
39 40 90 130 0,3076923077 2,19085775681281E-009 0,0946745549
40 34 85 119 0,2857142857 0,5 0,0459183673
41 42 59 101 0,4158415842 0,9999999978 0,3412410522
42 42 44 86 0,488372093 1 0,2617631152
43 37 22 59 0,6271186441 1 0,1390405056
44 48 13 61 0,7868852459 1 0,0454178984
45 65 12 77 0,8441558442 1 0,0242874009
46 51 8 59 0,8644067797 1 0,0183855214
47 72 0 72 1 1 0
48 64 4 68 0,9411764706 1 0,0034602076
49 53 4 57 0,9298245614 1 0,0049245922
50 86 0 86 1 1 0
51 80 2 82 0,9756097561 1 0,000594884
52 69 1 70 0,9857142857 1 0,0002040816
53 75 0 75 1 1 0
54 66 1 67 0,9850746269 1 0,0002227668

 

NOTA: Los valores de SL calculado, están calculados mediante la fórmula(1) y usando valores arbitrarios de c y L50

Y que el gráfico de todas estas muestras es similar al siguiente:

Dispersion

Y que necesitamos ajustar una curva sigmoide cuya fórmula(1) es:
SigmoideDonde L = Talla(cm), -c=Constante a calcular, L50=Cte a Calcular

Necesitamos calcular c y L50. Si c y L50 están ajustados lo lógico es que el error cuadrático se minimice. Entonces, en dos celas libres distintas, declaremos valores arbitrarios para c y L50 (digamos celdas J3 y J4 respectivamente) y en otra celda (J5) se coloca la suma de los errores cuadráticos (esto no es arbitrario y es la suma de la primer tabla). A mi se me ocurrió darle estos valores a C y L50, pero le pueden poner lo que quieran.

Variables C 20
L50 40
Error Total 1,1093543468

Como ven, el error cuadrático es bastante grande y nuestra idea es minimizarlo. Para ello, vamos a Herramientas/Solucionador. En la siguiente ventana ponemos como “Celda Objetivo” el error cuadrático $J$5, En “optimizar Resultados”, tildamos “minimizar” y en “Cambiando celdas” colocamos las variables a estimar, en este caso $J$3:$J$4. Le damos solucionar, esperamos un ratito y ya tenemos automágicamente  ¡La curva ajustada!

Curva Ajustada

Honorable Consejo Federal Pesquero fortifica la pesca del calamar

Primer Ministro Chino
Por segunda vez consecutiva, el CFP aprobó la presentación de proyectos pesqueros para incorporar 20 barcos poteros, nacionales o extranjeros, hasta por 30 años. Se argumenta que con 62 buques no se alcanzan los objetivos de aprovechamiento integral del recurso. La experiencia en 2012 fue un fracaso. Sigue leyendo

FELICITACIONES PARA ALEJANDRO FÁBREGA

Estimados compañeros y seguidores, tengo el grandísimo placer y agrado de comunicarles que nuestro queridísimo y gran amigo Alejandro Fábrega, se recibió ayer de Técnico Universitario Pesquero. En lo personal quiero hacer mención de que el Sr. Alejandro es un amigo único y de los mejores que podrían tener, una gran persona. Espero que esté sea un gran momento de felicidad para todos los que rodeamos a nuestro amigo.

IMG_7195

Cabe destacar que Alejandro se ha desarrollado luchando codo a codo con sus compañeros-amigos de la carrera y que ha logrado grandes cosas. Esperamos que este paso solo sea un aliciente para que sigas creciendo profesionalmente y que llegues muy lejos. Un abrazo gigante, y de nuevo, te felicito amigo.

¡Necesitamos colaboradores!

“Se buscan hombres para peligroso viaje. Salario reducido. Frío penetrante. Largos meses de completa oscuridad. Constante peligro. Dudoso regreso sano y salvo. Honor y reconocimiento en caso de éxito.” (Ernest Shackleton)

Estimados, docentes, estudiantes y amigos: Necesitamos gente con buena onda y ganas de aprender para escribir en el blog de la carrera. No hace falte ser erúditos en informática, ni tampoco en ninguna temática en general. Con buena predisposición y pilas nos sobra.

Si estás interesado, dejá un mensaje y a la brevedad nos comunicaremos con vos. Por favor difundir

Volvemos después de una larga pausa…

No queremos universidades para formar charlatanes y generalizadores. No queremos escuelas para formar hombres que les digan a los demás cómo hay que hacer las cosas sino hombres que sepan hacer por sí las cosas (…) y para esto hay que tener manos de trabajador y vivir con olor a aceite de las máquinas”.

(Juan D. Perón en la Inauguración de la Universidad Obrera Nacional)

El blog, este humilde espacio que alguna vez nos puso en la mira nacional e internacional, que generaba debates sobre temas de los que poco se discutía a nivel provincial, que trató de trasponer las barreras de las aulas para cumplir sus roles sociales de una forma mas efectiva, hoy vuelve con bríos renovados. Somos los estudiantes, profesores, porteros y el resto de la comunidad que no se rindieron cuando nos dijeron que Ingeniería Pesquera no iba mas en Ushuaia, cuando remábamos contra el anacronismo de un modelo que no contemplaba los desarrollos tecnológicos vinculados con la sociedad. Con nuestros logros y fracasos, seguimos adelante, y desde aquel inolvidable 12 de Octubre del 2009 hemos demostrado que se puede construir otra realidad a golpe de esfuerzos, sacrificios y trabajo.

Nuestro deseo mas grande es que este espacio no se personalice y que seamos todos los que lo construyamos. La Universidad somos todos, y todos de alguna u otra forma podemos colaborar para que esto crezca. El proyecto debe ser dinámico y adaptable a las realidades de la Universidad. Si se necesitan foros para cada materia, se crearán foros, si se hace imprescindible una herramienta de E-learning, la levantamos y y la dejamos funcionando. Las ganas y la fuerza las tenemos, solo esperamos que se sumen.

¡Un abrazo!