Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Nothing spectacular, just followed Martin's example from the least squared fit and produced this simple plot after using the problematic function from stats package. It was a bit trial-and-error, and it might be not so elegant-but it works. Regards, Radaovan Edited by user 04 January 2014 16:55:35(UTC)
| Reason: Not specified File Attachment(s): omorr attached the following image(s): |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
I added a plot with the draw package and tried to streamline the regression. The result of the statistics functions is a so called inference_result object, being essentially a list structure with a set of access functions. items_inference(obj) lists the available attributes, take_inference(attr,obj) returns the value of the specified attribute. An advantage of the draw package is that you can plot 2 column matrices just like in SMath. Thus the example works without a single list<>matrix conversion. File Attachment(s): mkraska attached the following image(s): |
|
1 user thanked mkraska for this useful post.
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Thank you Martin very much, It does really work and it is quite nice . I think I understand a bit more how the things are working The "inference_results" thing is a good one as well. I do not now how much time it would take to figure out myself. Thank you. Regards, Radovan |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Here is how the gaps in StatisticalTools by Davide could be filled. The quantile functions, however, are rather slow. I don't know if that is due to slow translation or to the slow algorithm. Also, the chi2-quantile seems to not be robust (sometimes fails to be plotted). The function definitions in the attached sheet are not entirely consistent to Davide's plugin, as they do not produce a 2 column matrix when they receive vector input. Davide, if you do not want to use the Maxima algorithms in your plugin (which might be more efficient), then we could provide the functions as part of the Maxima-plugin. File Attachment(s): mkraska attached the following image(s): |
|
1 user thanked mkraska for this useful post.
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Originally Posted by: mkraska Davide, if you do not want to use the Maxima algorithms in your plugin (which might be more efficient), then we could provide the functions as part of the Maxima-plugin. Thank you Martin for pointing this out once more. As I know, Davide asked for some appropriate algorithms to fill the gap in the StatisticalTools. These quantile functions which missed there are one of the most frequently used. It would be rather useful if they are inside Statistical tools, although we have them via Maxima as Martin presented. Regards, Radovan |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Newbie
Groups: Registered
Joined: 21/03/2014(UTC) Posts: 9
|
Would like to use the plot/draw functions from maxima in the executable file generated by the Smath viewer.
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Originally Posted by: baijo Would like to use the plot/draw functions from maxima in the executable file generated by the Smath viewer. You might file your feature request in the bug tracking system. Please attach an .sm sheet for testing. However, providing Maxima as part of compiled software might violate their license terms. |
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Here is a small program for generation of Ashby-like material property charts. The data has been taken from the online ressources for Buchmayr 2002. Besides the Maxima plugin and it's drawing functions this program also demonstrates the use of varous plugins by Davide (Table, Slider, DataExchange, Checkbutton List) It would be easy to generate the bubbles in the modeller region. But that has to wait until Andrey finds some time to document how to place text annotations. Then one might be able to click the bubbles in order to create and place label strings. EDIT: Due to a checkboxlist bug I had to remove the material chooser. For now you have to switch on and off the materials by setting the corresponding entry in gr.sel to 1 or 0. Edited by user 18 April 2014 03:21:23(UTC)
| Reason: Not specified File Attachment(s): mkraska attached the following image(s): |
|
2 users thanked mkraska for this useful post.
|
on 18/04/2014(UTC), on 18/04/2014(UTC)
|
|
Rank: Advanced Member Groups: Registered, Advanced Member Joined: 13/01/2012(UTC) Posts: 2,650 Location: Italy Was thanked: 1338 time(s) in 878 post(s)
|
Thanks Martin, nice work CheckBoxList is under development (alpha/not public - in this state it's "naturally" buggy) Best regards, Davide edit: for the testers, there's a plugin in the EM to open and save again worksheets containing ComboBoxList regions. Edited by user 18 April 2014 20:29:30(UTC)
| Reason: Not specified |
If you like my plugins consider to support SMath Studio buying a plan; to offer me a coffee: paypal.me/dcprojects |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Improved version of the Ashby chart generator. - bubbles are elliptic - minimal bubble size independent of axes limits - separate checkboxes for selection of material groups File Attachment(s): mkraska attached the following image(s): |
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Here is a modified version of Radovan's distrib-example. Changes were required, because of - The plugin does not convert Maxima lists to SMath vectors any more, - distrib functions seem to require float arguments. Edited by user 28 May 2014 19:32:27(UTC)
| Reason: Not specified File Attachment(s): mkraska attached the following image(s): |
|
1 user thanked mkraska for this useful post.
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Thank you Martin for correcting the example. Here is another one using F distribution from Maxima distrib package. Statistical plugin by Davide is used as well. There is a nice feature in this plugin of returning two columns when using probability distribution functions (missing quantile functions unfortunately). One more interesting thing is to note about distrib package functions. They will not use vectors as arguments returning vector results - therefore we have to use loops (see the second picture). This is contrary to other Maxima examples where we can avoid loops in Smath Regards, Radovan Edited by user 29 May 2014 12:25:26(UTC)
| Reason: Not specified File Attachment(s): omorr attached the following image(s): |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
1 user thanked omorr for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Originally Posted by: omorr One more interesting thing is to note about distrib package functions. They will not use vectors as arguments returning vector results - therefore we have to use loops (see the second picture). This is contrary to other Maxima examples where we can avoid loops in Smath
Much like the SMath plugins, many Maxima features are provided by third party packages (wrt the core development team) and there are as many opinions on what a good function interface is as there are contributors. It is surely possible to re-define the Maxima functions such that they accept matrix arguments and return two-col matrices in that case. If you have a corresponding idea, then you can add the definition to your maxima.xml and try if it works. We can then include this into the plugin distribution. I added a Maxima based plot to the sheet. File Attachment(s): mkraska attached the following image(s): |
|
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Thank you Martin for this one . You read my mind, as I was browsing your Maxima drawing examples reminding the Draw2D options and trying to do a similar thing. Now, you've made me a template. Thank you. About the changing maxima.xml, unfortunately it is beyond my capabilities at the moment. I am far from being an advanced user of these things - just a curious one. Regards, Radovan |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Here is another simple example of using SMath and Maxima plugin for a textbook example. One has to admit that teaching with SMath might be a fun Regards, Radovan Edited by user 30 May 2014 14:48:47(UTC)
| Reason: Not specified File Attachment(s): omorr attached the following image(s): |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Here is an example (LRC oscillator) for solving differential equations from the Maxima mailing list. The combination of SMath Studio and Maxima has much better and more convenient support for units than native Maxima. The first version solves the generic system of differential and linear algebraic equations directly. A second version is to follow in the next post. Edited by user 05 June 2014 20:31:35(UTC)
| Reason: Not specified File Attachment(s): mkraska attached the following image(s): |
|
1 user thanked mkraska for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
And here is a version of the previous example where we reduce the system to a second order ODE using Maxima's ev() function. Due to spurious imaginary parts in the numeric evaluation of the result, we have to use Maxima again for extraction of the real part. File Attachment(s): mkraska attached the following image(s): |
|
1 user thanked mkraska for this useful post.
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Nice example of desolve(), especially pointing out using of units. Thank you Martin. Moreover, I just forgot this trick about setting x-range on the native Smath 2D plot. Regards, Radovan omorr attached the following image(s): |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
|
|
|
Rank: Administration Groups: Registered, Advanced Member Joined: 23/06/2009(UTC) Posts: 1,740 Was thanked: 318 time(s) in 268 post(s)
|
Just testing desolve() based on the example from Maxima by Example. I am not sure that desolve() can solve other than initial conditions for t=0. Regards, Radovan Edited by user 06 June 2014 14:03:34(UTC)
| Reason: Not specified File Attachment(s): omorr attached the following image(s): |
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!" |
1 user thanked omorr for this useful post.
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,988 Was thanked: 1127 time(s) in 723 post(s)
|
Originally Posted by: omorr I am not sure that desolve() can solve other than initial conditions for t=0.
Maxima doc says that initial conditions can only be specified for x=0. One might wonder why then t has to be specified in atvalue(), but that function has more uses than just specification of initial conditions. BTW, you also can use ODE.2, a wrapper for ode2, with preprocessing of the equation switched off. Also, you can specify the initial conditions in Solve using Davide's at(). File Attachment(s): mkraska attached the following image(s): |
|
|
|
|
Forum Jump
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.