Welcome Guest! To enable all features please Login. New Registrations are disabled.

Notification

Icon
Error

Login


2 Pages12>
Options
Go to last post Go to first unread
Offline - Norm -  
#1 Posted : 07 September 2018 20:02:40(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
I have tried to convert the attached MC15 worksheet to SMATH. I cannot work around the solve block (Maximize) successfully. I welcome the collabs review and recommendation(s).
Sincerely,
nas
Chain Failure Weibull Analysis.sm (64kb) downloaded 15 time(s). chain failure weibull analysis(1).xmcd (247kb) downloaded 9 time(s).

Wanna join the discussion?! Login to your SMath Studio Forum forum account. New Registrations are disabled.

Offline - Norm -  
#2 Posted : 07 September 2018 20:19:42(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
Solve block.jpg

The solve block is in the "Weibull Distribution" area (collapsed).
nas
Offline Jean Giraud  
#3 Posted : 07 September 2018 22:45:21(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: nasch Go to Quoted Post
I have tried to convert the attached MC15 worksheet to SMATH.


Smath converts zero, absolute zero MCD.
You have to create the project from square one .
At least, you can plot 3 graphs
Offline - Norm -  
#4 Posted : 08 September 2018 00:26:00(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
If I force the shape and scale factors all seems to work. I just need help with obtaining these. Chain Failure Weibull Analysis Forced.sm (67kb) downloaded 7 time(s).
Offline - Norm -  
#5 Posted : 08 September 2018 04:06:47(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
Chain Failure Weibull Analysis Forced.sm (72kb) downloaded 17 time(s).

I noticed that my previous SMATH files lacked some narrative. Updated now.

Thank you,
Offline - Norm -  
#6 Posted : 09 September 2018 20:35:14(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
Chain Failure Weibull Analysis - 2.sm (53kb) downloaded 11 time(s).

Taking jmg's advice; started going thru the steps and found a function that didn't come over from MC that I used the wrong "Sum" function in SMATH. Worksheet function ll.W(....) now seems to work plotting varying lambda and beta (sheet attached). I still cannot determine/maximize the function, finding the optimim lambda and beta.

I welcome any suggestions.

nas
Offline Razonar  
#7 Posted : 10 September 2018 02:57:40(UTC)
Razonar


Rank: Advanced Member

Groups: Registered
Joined: 28/08/2014(UTC)
Posts: 1,356
Uruguay

Was thanked: 815 time(s) in 516 post(s)
Hi. Try this. Chain Failure Weibull Analysis - 3.sm (66kb) downloaded 16 time(s).

Alvaro.
Offline Jean Giraud  
#8 Posted : 10 September 2018 04:31:47(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: nasch Go to Quoted Post
I still cannot determine/maximize the function, finding the optimim lambda and beta.


Is that what you mean: that the n*F.W(x,λ,β) is the best fit to data {TY] ?
Very rarely, a model function can be declared exact. That is the case of all LS based methods.
I guess the data set [TY] can be fitted by many models, surely giving λ,β closely to Alvaro.
Easy to plot [TY] over the n*F.W(x,λ,β) ... module added .
Interesting but not of immediate understanding for the novice.

Cheers ... Jean
Opened in 6179, re-posted from same original.

Chain Failure Weibull Analysis - 3.sm (72kb) downloaded 21 time(s).
Offline mkraska  
#9 Posted : 10 September 2018 12:20:39(UTC)
mkraska


Rank: Advanced Member

Groups: Registered
Joined: 15/04/2012(UTC)
Posts: 1,986
Germany

Was thanked: 1124 time(s) in 721 post(s)
Fit using Maxima's lsquares package. MSE() is the mean square error, which is minimized by Fit()

2018-09-10 11_11_54-SMath Studio 0.99.6671 - [Chain%20Failure%20Weibull%20Analysis%20-%203.sm_].png
Chain%20Failure%20Weibull%20Analysis%20-%203_Kr.sm (79kb) downloaded 26 time(s).
Martin Kraska

Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx
thanks 1 user thanked mkraska for this useful post.
on 10/09/2018(UTC)
Offline - Norm -  
#10 Posted : 10 September 2018 15:44:02(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
Alvaro, jmg and Martin, many thanks!

Martin's requires a newer version of SMATH installed, I'll have to request the upgrade to be installed.
Offline mkraska  
#11 Posted : 10 September 2018 16:10:51(UTC)
mkraska


Rank: Advanced Member

Groups: Registered
Joined: 15/04/2012(UTC)
Posts: 1,986
Germany

Was thanked: 1124 time(s) in 721 post(s)
Originally Posted by: nasch Go to Quoted Post
Alvaro, jmg and Martin, many thanks!

Martin's requires a newer version of SMATH installed, I'll have to request the upgrade to be installed.


Should work with SMath build 6671 as well (latest stable with updated plug-ins), just ignore the complaints.
Martin Kraska

Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx
Offline Jean Giraud  
#12 Posted : 10 September 2018 18:40:07(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: nasch Go to Quoted Post
Martin's requires a newer version of SMATH installed,


Maybe NOT ... Maybe I'm too ignorant about Weibull.
What I understand: over 'T', you have observed cumulative failure, is it ?
From there you can approximate the probability of failures from the logistic Verhulst.
I would be curious you find simpler and more realistic and more immediate insight.
Done Smath 6971, should work in more recent SS.

Cheers ... Jean

Page0Weibull.sm (41kb) downloaded 20 time(s).

Offline Jean Giraud  
#13 Posted : 11 September 2018 04:10:19(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
What I understand: over 'T', you have observed cumulative failure, is it ?


Not like said. The 'Y' data are the days between failure.
Refreshed and added the cumulative integration of the Verhulst fit.

Genfit Weibull.sm (46kb) downloaded 24 time(s).
Offline - Norm -  
#14 Posted : 11 September 2018 21:41:46(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
Thank you Jean, I have a lot to look over. I noticed that the MC file calculated a beta > than 1, which is a "wear-out" characteristic vs. the SMATH beta <1; "infant mortality". I found an earlier MC file that used a different method to find beta, and that too beta was < 1. Using excel, that was also in agreement; beta < 1. The values vary, the old MC file = 0.99, 0.98 for the SMATH file, and .99 for excel.

I use the files as a reference, being a novice using tools for stats. I did get a lot of new ideas and understand more using SMATH, which I am thankful to this collab.

Sincerely,

Norm
Offline - Norm -  
#15 Posted : 11 September 2018 22:03:52(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
Hello Jean, Beta will not calculate, I have the same SMATH version 0.99.6671. Maybe I'm missing a plugin?
Offline Jean Giraud  
#16 Posted : 11 September 2018 23:26:15(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: nasch Go to Quoted Post
"infant mortality"

... because there 3 kinds of lies:
1. Kid's white lie [Mamy, Mamy ... the cat did suck your last candy]
2. True lies [Politician's lie]
3. Statistics

What happens in the β "infant mortality" is consequence of
reflexivity in the function wrt parameters. Many models
can't be concluded by NO maths ... just manual from St. Patience !

If you stick with Weibull failure analysis, good luck.
In about ½ hr I found nothing from Gooling.

Offline - Norm -  
#17 Posted : 12 September 2018 00:24:22(UTC)
- Norm -


Rank: Advanced Member

Groups: Registered
Joined: 29/11/2017(UTC)
Posts: 82
United States
Location: Midwest

Was thanked: 9 time(s) in 7 post(s)
Sorry Jean, your genfit file won't find beta. It works fine when I download, but fails when I open on my pc.

nas
Offline Jean Giraud  
#18 Posted : 12 September 2018 06:27:03(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: nasch Go to Quoted Post
Sorry Jean, your genfit file won't find beta. It works fine when I download, but fails when I open on my pc.


1. The work sheet attached finds the β's of the Verhulst logistic function,
which corresponds to your collection. For the two parameters Weibull, we have
to scale the Y data to [0..1] and the X data to something else.
Your project up until now is just an undefined collection.
I will attempt a pure Weibull fit.

2. Can't understand it works fine from download but fails on your PC.
Like saying: moon is flat because she is flat as well as spherical
because she is spherical. Just circular like Shrodinger cat ...
hit the cat he is dead, hit again back live.

Cheers ... Jean

Stat Treasury_7 Weibull.sm (10kb) downloaded 10 time(s).
Offline Jean Giraud  
#19 Posted : 12 September 2018 06:55:35(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
I will attempt a pure Weibull fit.


... here is a typical two parameters Weibull

Stat Treasury_7 Weibull.sm (19kb) downloaded 13 time(s).
Offline Jean Giraud  
#20 Posted : 12 September 2018 07:17:41(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
... here is a typical two parameters Weibull


... too reactive for :Genfit Minimize.
Thus: simply manual fit by the look.

Stat Treasury_7 Weibull.sm (21kb) downloaded 21 time(s).
Users browsing this topic
Guest
2 Pages12>
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.