Riccardo, thanks for replying. Let me provide more information , as I reckon I did a bad job the first time ;)
Thanks for the welcoming! The tricky bit of analysis will handle angles, distances and deltas of those variables, since everything needs to be georeferenced, every observation takes a coordinate pair (X,Y or N,E) but also Z (X,Y,Z). The data will take the form of a panel model, with a benchmark value for the base variables (the average of a few dozen observations ) and the recurring measurements shall be compared against the benchmark, and once enough data is collected, we'll try to understand various relationships amongst variables and their observed values. My emphasis on the (XYZ) is that though rain, temp, exogenous impacts, etc play a role in the *seasonal* observations therefore recurring, the coordinates are observed down to a high precision (mm) and changes in them may flag changes in the observed values which are not related to any recurring cycle but a deterioration in the condition of a specific location, and many tests and hypothesis will be checked for each observation. There geo variable should be so interrelated that I expected y= a + b*sin(x) is a generalization for an accumulation of a geographical position, such as Y = Yi + d*sine(alpha) for example. When I run an OLS the results seem to explain very little of what I expected to be a straightforward and linear expression, specially because it's a panel and I am comparing the same data which is expected to change very little between observations/measurements. So first of all I need to script for checking 1mm change in coordinates (above or below each variable X, Y, Z) then check for season impacts etc. I am attaching an example file with test measurements. Important detail (in order to make sense of data) is that I am creating a BENCHMARK value for XYZ (average value) and creating variables DeltaX/Y/Z which are the ones I will compare against 0.001m (1mm) change by doing ex. Delta_X(i) = X(i) - mean(X) etc. I also should say that the dataset used a less precsa equipment so it is very hard to check for 1mm changes, as per OLS , the results will tell. My idea is to use tests to check if the changes breach a critical value and if they relate to seasonal factors (no big deal) or they are indeed red flags. Any ideas how you'd run this? Cheers from Brazil! Em sex., 8 de dez. de 2023 às 11:51, Riccardo (Jack) Lucchetti < p002...@staff.univpm.it> escreveu: > On 08/12/2023 14:52, Meridiana GeoTopo wrote: > > Hello Gretl Enthusiasts! > > > > I hope this message finds you well. I am reaching out to this > > knowledgeable community as I embark on a unique project that involves > > utilizing Gretl for the estimation of environmental and engineering > > parameters. While I've already dived into Gretl for this purpose, the > > project's diversity from traditional econometrics poses a few challenges. > > Welcome to the community! > > > Specifically, I am looking to run regression models where I can create > > custom formulas and expressions. To give you a glimpse, my requirements > > involve working with a Cartesian coordinate system, something like y = a > > + d(sine(b)), as I heavily rely on spatial data in my study and also > > structural health. > > Some clarification is needed here. If (as I believe form the context) a > and b are parameters and d is an observable series, then all you need is > OLS plus the application of the delta method in order to recover the > estimate of b and its standard error. Otherwise, if b is an observable > series, then sin(b) is an observable series too, and a distinction must > be made according to what the "d()" expression means. If it means "a > parameter called d multiplied by something", then OLS suffices, again. > Otherwise, if it means "the function d() applied to something", where > this function continuous and differentiable, then what you need is NLS. > > All these methods are provided by gretl, so feel free to ask for guidance. > > > I could also use TLS and GLS models - can I run them in gretl and if so > > where can I find examples of implementations? > > The TLS model as such is not implemented in gretl natively, since it's > very seldom used (if ever) in econometrics. However, we do have a > function for performing the SVD and it shouldn't be difficult to write a > script to perfom TLS estimation. As for GLS, it really depends on what > you need (GLS is a very generic term that could be applied to many > different things). If what you need is the special case known in > econometrics as "weighted least squares", then yes, we do have a "wls" > command. > > > Later we will attach soil, air, temperature, rain etc to the model but > > these will follow a linear model approach easily. > > > > As a newcomer to Gretl, I am seeking honest assistance, guidance, and > > insights from the experienced members of this group. If anyone could > > spare a moment to share your expertise, provide script help, or guide me > > on using the GUI to achieve these objectives, I would greatly appreciate > it. > > > > Your support will not only contribute significantly to the success of my > > project but will also help me enhance my understanding of Gretl in a > > practical context. I am eager to learn and absorb insights from the > > Gretl community. > > > > Thank you in advance for your time and assistance. I look forward to the > > valuable input from this esteemed group. > > Feel free to ask! > > ------------------------------------------------------- > Riccardo (Jack) Lucchetti > Dipartimento di Scienze Economiche e Sociali (DiSES) > > Università Politecnica delle Marche > (formerly known as Università di Ancona) > > r.lucche...@univpm.it > http://www2.econ.univpm.it/servizi/hpp/lucchetti > ------------------------------------------------------- > _______________________________________________ > Gretl-users mailing list -- gretl-users@gretlml.univpm.it > To unsubscribe send an email to gretl-users-le...@gretlml.univpm.it > Website: > https://gretlml.univpm.it/postorius/lists/gretl-users.gretlml.univpm.it/ >
ORIENT EQUIP STN CAMPANHA PT_CTR PONTO DESCR AH AV DH X Y Z Irrad ARC5 STN1 09/nov/23 A A1 C 13,80916667 87,68638889 149,0954 1.035,59 5.144,79 17,7138 IrradInv ARC5 STN1 09/nov/23 A A2 C 13,81638889 87,69083333 149,0988 1.035,61 5.144,78 17,7023 Irrad ARC5 STN1 09/nov/23 A A3 C 13,81 87,68555556 149,0963 1.035,59 5.144,79 17,716 IrradInv ARC5 STN1 09/nov/23 A A4 C 13,81527778 87,69138889 149,0989 1.035,60 5.144,79 17,7009 Irrad ARC5 STN1 09/nov/23 A A5 C 13,81055556 87,685 149,0952 1.035,59 5.144,78 17,7174 IrradInv ARC5 STN1 09/nov/23 A A6 C 13,81611111 87,69055556 149,0988 1.035,61 5.144,78 17,703 Irrad ARC5 STN1 09/nov/23 A A7 C 13,81 87,68472222 149,0972 1.035,59 5.144,79 17,7182 IrradInv ARC5 STN1 09/nov/23 A A8 C 13,81444444 87,69138889 149,0999 1.035,60 5.144,79 17,7009 Irrad ARC5 STN1 09/nov/23 A A9 C 13,80888889 87,68583333 149,0953 1.035,59 5.144,79 17,7152 IrradInv ARC5 STN1 09/nov/23 A A10 C 13,81583333 87,69055556 149,0998 1.035,61 5.144,79 17,7031 Irrad ARC5 STN1 09/nov/23 B B1 C 79,50916667 100,3155556 7,8559 1.007,72 5.001,43 10,2601 IrradInv ARC5 STN1 09/nov/23 B B2 C 79,51444444 100,3213889 7,8558 1.007,72 5.001,43 10,2593 Irrad ARC5 STN1 09/nov/23 B B3 C 79,51 100,3163889 7,8559 1.007,72 5.001,43 10,26 IrradInv ARC5 STN1 09/nov/23 B B4 C 79,51333333 100,3219444 7,8548 1.007,72 5.001,43 10,2594 Irrad ARC5 STN1 09/nov/23 B B5 C 79,51166667 100,3097222 7,8551 1.007,72 5.001,43 10,2611 IrradInv ARC5 STN1 09/nov/23 B B6 C 79,51555556 100,3125 7,855 1.007,72 5.001,43 10,2607 Irrad ARC5 STN1 09/nov/23 B B7 C 79,51055556 100,3072222 7,8552 1.007,72 5.001,43 10,2615 IrradInv ARC5 STN1 09/nov/23 B B8 C 79,51555556 100,3133333 7,854 1.007,72 5.001,43 10,2608 Irrad ARC5 STN1 09/nov/23 B B9 C 79,50888889 100,3094444 7,8561 1.007,72 5.001,43 10,261 Irrad ARC5 STN1 09/nov/23 C C1 C 154,0019444 91,05888889 22,7361 1.009,97 4.979,56 11,2698 IrradInv ARC5 STN1 09/nov/23 C C2 C 154,0088889 91,06305556 22,7351 1.009,96 4.979,56 11,2681 Irrad ARC5 STN1 09/nov/23 C C3 C 154,0038889 91,05777778 22,7351 1.009,97 4.979,57 11,2702 IrradInv ARC5 STN1 09/nov/23 C C4 C 154,0083333 91,06305556 22,7351 1.009,96 4.979,56 11,2681 Irrad ARC5 STN1 09/nov/23 C C5 C 154,0058333 91,05861111 22,7351 1.009,96 4.979,56 11,2699 IrradInv ARC5 STN1 09/nov/23 C C6 C 154,0097222 91,06277778 22,7351 1.009,96 4.979,56 11,2682 Irrad ARC5 STN1 09/nov/23 C C7 C 154,0030556 91,05861111 22,7351 1.009,97 4.979,57 11,2699 IrradInv ARC5 STN1 09/nov/23 C C8 C 154,0086111 91,06333333 22,7351 1.009,96 4.979,56 11,268 Irrad ARC5 STN1 09/nov/23 C C9 C 154,0036111 91,05805556 22,7361 1.009,97 4.979,56 11,2701 IrradInv ARC5 STN1 09/nov/23 C C10 C 154,01 91,06444444 22,7351 1.009,96 4.979,56 11,2676 Irrad ARC5 STN1 10/nov/23 A A50 C 13,80694444 87,68555556 149,0943 1.035,58 5.144,79 17,7159 IrradInv ARC5 STN1 10/nov/23 A A51 C 13,81138889 87,6925 149,097 1.035,59 5.144,79 17,6979 Irrad ARC5 STN1 10/nov/23 A A52 C 13,80722222 87,68555556 149,0953 1.035,58 5.144,79 17,7159 IrradInv ARC5 STN1 10/nov/23 A A53 C 13,81138889 87,69111111 149,0979 1.035,59 5.144,79 17,7016 Irrad ARC5 STN1 10/nov/23 A A54 C 13,80694444 87,68527778 149,0952 1.035,58 5.144,79 17,7167 IrradInv ARC5 STN1 10/nov/23 A A55 C 13,81222222 87,69166667 149,0999 1.035,60 5.144,79 17,7002 Irrad ARC5 STN1 10/nov/23 A A56 C 13,8075 87,68527778 149,0952 1.035,58 5.144,79 17,7167 IrradInv ARC5 STN1 10/nov/23 A A57 C 13,81277778 87,69222222 149,097 1.035,60 5.144,79 17,6986 Irrad ARC5 STN1 10/nov/23 A A58 C 13,80805556 87,68527778 149,0962 1.035,58 5.144,79 17,7167 IrradInv ARC5 STN1 10/nov/23 A A59 C 13,8125 87,69194444 149,0979 1.035,60 5.144,79 17,6994 Irrad ARC5 STN1 10/nov/23 A A60 C 13,80666667 87,685 149,0942 1.035,58 5.144,79 17,7173 IrradInv ARC5 STN1 10/nov/23 A A61 C 13,81222222 87,69111111 149,0989 1.035,60 5.144,79 17,7016 Irrad ARC5 STN1 10/nov/23 B B50 C 79,5125 100,3097222 7,8551 1.007,72 5.001,43 10,2611 IrradInv ARC5 STN1 10/nov/23 B B51 C 79,51888889 100,3169444 7,8539 1.007,72 5.001,43 10,2603 Irrad ARC5 STN1 10/nov/23 B B52 C 79,51305556 100,3088889 7,8551 1.007,72 5.001,43 10,2612 IrradInv ARC5 STN1 10/nov/23 B B53 C 79,52 100,3161111 7,854 1.007,72 5.001,43 10,2604 Irrad ARC5 STN1 10/nov/23 B B54 C 79,51277778 100,3108333 7,8561 1.007,72 5.001,43 10,2608 IrradInv ARC5 STN1 10/nov/23 B B55 C 79,51833333 100,3161111 7,854 1.007,72 5.001,43 10,2604 Irrad ARC5 STN1 10/nov/23 B B56 C 79,51305556 100,31 7,8561 1.007,72 5.001,43 10,2609 IrradInv ARC5 STN1 10/nov/23 B B57 C 79,52027778 100,3161111 7,854 1.007,72 5.001,43 10,2604 Irrad ARC5 STN1 10/nov/23 B B58 C 79,5125 100,3113889 7,857 1.007,73 5.001,43 10,2605 IrradInv ARC5 STN1 10/nov/23 B B59 C 79,51805556 100,3166667 7,8539 1.007,72 5.001,43 10,2603 Irrad ARC5 STN1 10/nov/23 C C50 C 154,0027778 91,0525 22,7352 1.009,97 4.979,57 11,2723 IrradInv ARC5 STN1 10/nov/23 C C51 C 154,0094444 91,05666667 22,7351 1.009,96 4.979,56 11,2707 Irrad ARC5 STN1 10/nov/23 C C52 C 154,0033333 91,05194444 22,7352 1.009,97 4.979,57 11,2725 IrradInv ARC5 STN1 10/nov/23 C C53 C 154,0080556 91,05583333 22,7351 1.009,96 4.979,56 11,271 Irrad ARC5 STN1 10/nov/23 C C54 C 154,0041667 91,05138889 22,7352 1.009,97 4.979,57 11,2728 IrradInv ARC5 STN1 10/nov/23 C C55 C 154,0088889 91,05722222 22,7351 1.009,96 4.979,56 11,2704 Irrad ARC5 STN1 10/nov/23 C C56 C 154,005 91,05111111 22,7342 1.009,96 4.979,57 11,2729 IrradInv ARC5 STN1 10/nov/23 C C57 C 154,0091667 91,05722222 22,7341 1.009,96 4.979,57 11,2705 Irrad ARC5 STN1 10/nov/23 C C58 C 154,0047222 91,05111111 22,7342 1.009,96 4.979,57 11,2729 IrradInv ARC5 STN1 10/nov/23 C C59 C 154,0091667 91,05666667 22,7341 1.009,96 4.979,57 11,2707
_______________________________________________ Gretl-users mailing list -- gretl-users@gretlml.univpm.it To unsubscribe send an email to gretl-users-le...@gretlml.univpm.it Website: https://gretlml.univpm.it/postorius/lists/gretl-users.gretlml.univpm.it/
