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

Bancamos a los 3 porteros ilegítimamente expulsados de la Universidad Obrera Ushuaia

Estimadxs: Pasamos a comunicarles que tres históricos porteros (Juan Carlos, Carlos y Enrique) de nuestra Universidad Obrera, fueron suspendidos de sus puestos laborales por pedir simplemente que se respeten derechos adquiridos: el pago de horas extras y licencias por fallecimiento de familiares.
Como estudiantes nos sentimos tremendamente heridos y por ello, es que con el único objetivo de apoyar a aquellos que son compañeros de ruta de nuestra vida académica, queremos invitarlos a apoyar su lucha.

Creemos que los integrantes de la Universidad Obrera debemos tomar partido de alguna forma y los invitamos a todos a una reunión para sacar una solicitada en conjunto con UTEDyC.

DATOS
———–
DIA: 18/08/13
HORARIO: 17.00

LUGAR: Sede UTEDyC -Gobernador Paz 104 (esquina Antártida Argentina)-

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

Sistema Vito de Navegación Marítima

Presentación: MIERCOLES 4 de setiembre  14:45 hs – CADIC
El sistema de Navegación “VITO”, desarrollado por estudiantes de Ingeniería Pesquera de la UTN Ushuaia, es el resultado de la iniciativa y trabajo en conjunto de diversas personas a lo largo de los últimos diez años.
En un primer momento, Oficiales y Tripulantes del Buque “Ushuaia”, de la Empresa Antarpply (única empresa argentina de expedición antártica), empezaron a realizar levantamientos hidrográficos expeditivos en sectores de la Antártida, en que los peligros a la navegación se encontraban deficientemente indicados en las cartas náuticas vigentes.
Dicha iniciativa, orientada a proveer a la seguridad en la navegación del propio Buque, se plasmaba en cartas “papel” que se dibujaban a mano:
En un paso posterior, las cartas comenzaron a ser elaboradas a bordo en formato digital, mediante la búsqueda y empleo de softwares libres que permitieran, a través del acceso al código fuente, el dibujo de contornos de costa e inclusión de profundidades e indicaciones necesarias para la derrota segura:
A esta altura, la digitalización se cristalizaba en la elaboración en sí de la carta, pero su empleo seguía siendo en versión papel. Surgieron entonces dos requerimientos nuevos y complementarios:1) Poder “correr” las cartas de aquellos lugares difíciles, en un soft de navegación.

2) Vincular los puntos notables de las áreas más críticas, es decir aquellas más cercanas a bajofondos, a las imágenes visual y radar que pudiera llegar a tener el navegante que se aproximara al lugar.  De esta manera, se buscaba minimizar la incertidumbre que pudiera originar la modificación de líneas de costa por presencia de hielo, propendiendo a su vez que fuera posible situar con seguridad al buque con multiplicidad de sistemas de posicionamiento. Es decir: no solo con GPS, sino también y muy principalmente con radar y en visual. La imagen que sigue describe perfectamente el sentido de este nuevo “requisito” a la cartografía antártica:

En este punto se recurrió al “GRUPO DE TECNOLOGIA MARITIMA UTN Ushuaia”, agrupación conformada por estudiantes de Ingeniería Pesquera y marinos mercantes del Centro de Patrones Fluviales, de Pesca y Cabotaje Marítimo.

El grupo venía trabajando en diversos proyectos vinculados a campañas oceanográficas, particularmente para la Comisión Nacional del Límite Exterior de la Plataforma Continental. En especial, venía haciendo experiencia en  cuestiones asociadas a las maniobras, equipos de batimetría y de muestreo de fondo del Buque ARA “PUERTO DESEADO” del CONICET.

Los estudiantes propusieron entonces desarrollar un sistema integrado de navegación, pero a partir de un entorno de software libre.
Luego de un año de trabajo, en mayo de 2012 se expuso el proyecto en el marco del “2do Taller de Operatoria Marítima Antártica”. A esa fecha se había avanzado, entre otros aspectos, en la elección y adaptación del software libre que serviría de base para correr las cartas y, particularmente, en la transformación de éstas, pasándolas del soporte digital en que eran elaboradas, a un sistema en el que pudieran ser trabajadas con los diferentes sensores del Buque.
La integración del soft de navegación con los sensores se pensó través del sistema “VITO”, desarrollado por los propios estudiantes y denominado así en recuerdo del gran navegante argentino Vito DUMAS, primer hombre en circunnavegar en solitario la Antártida por la “Ruta Imposible”. El “VITO”, de este modo, no solo permitiría navegar en modo digital con la cartografía levantada por el Buque antártico “Ushuaia”, sino que a su vez constituiría el primer sistema de navegación integral elaborado específicamente, en la Argentina, para satisfacer requerimientos de la navegación en zonas polares.
Con el levantamiento de nuevas cartas en zonas críticas durante la temporada 2013, y con vistas a su puesta a punto, calibración y pruebas, la empresa Antarpply autorizó la instalación del VITO a bordo del “Ushuaia”, en cuanto el sistema se encontrara en un estado de desarrollo que permitiera su interconexión con los sensores del Buque.
A tal efecto, el grupo de Tecnología Marítima se vinculó con Gustavo Biancalana, Técnico Electrónico especializado en ecosondas, radares y sistemas de navegación, con gran experiencia en proveer soluciones tecnológicas a embarcaciones.
El resultado ha sido un nuevo avance hacia el desarrollo final del sistema, el cual es esperable que se encuentre en etapa efectiva de prueba durante el presente año, en aguas del Canal Beagle, y con muchas probabilidades de ser incorporado a bordo del Ushuaia en la temporada 2014/15 (si no llega a concretarse durante la que se inicia próximamente, 2013/14).
El sistema, finalmente, se presenta el día miércoles 4 de setiembre, a 14:45 hs, en el CADIC, en el marco de la SEMANA DEL EMPRENDEDOR TECNOLOGICO (jornada cuyas actividades comienzan a 14:00):
Desde el Centro de Patrones invitamos a participar de la exposición que brindarán los estudiantes del “Grupo de Tecnología Marítima de la UTN Ushuaia”, quienes, como se ha expresado, en trabajo en conjunto con marinos mercantes, empresa antártica argentina y especialistas en equipamiento electrónico, están aportando y apostando con su propio esfuerzo al desarrollo local, materializando un producto que contribuirá a la seguridad en la navegación, y la continuidad de la operatoria antártica de turismo con base en Ushuaia.———————————————————————————————————————————-
El Centro de Patrones y Oficiales Fluviales, de Pesca y Cabotaje Marítimo, junto a los estudiantes de la UTN Ushuaia, agradecen al Ministerio de Industria e Innovación Productiva de Tierra del Fuego, el interés puesto de manifiesto por la labor del Grupo de Tecnología Marítima.

Agradecemos a su vez, tanto a dicho Ministerio como a la UNTDF, la invitación a participar de la SEMANA DEL EMPRENDEDOR TECNOLOGICO.

Matriz para el cálculo de redes de pesca de arrastre

Estoy trabajando en el desarrollo de Software Libre de cálculo de redes de pesca de arrastre escrito en Python.

Después de mucho pensar un método de cálculo sistematizado, llegué a la conclusión de que, para determinar parámetros básicos, tales como Superficie Sólida y Resistencia de la red, se puede hacer uso de una matriz como la siguiente:

Matriz de red 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.