Rank: 1739 based on 87 downloads (last 30 days) and 3 files submitted
photo

Willem-Jan de Goeij

E-mail

Personal Profile:

n

Professional Interests:

 

Watch this Author's files

 

Files Posted by Willem-Jan View all
Updated   File Tags Downloads
(last 30 days)
Comments Rating
04 Jun 2011 Screenshot PARFOR Progress Monitor v2 This version is an update of the progress monitor made by Edric Ellis. Author: Willem-Jan de Goeij parallel, parfor, progress, monitor 53 11
  • 4.5
4.5 | 8 ratings
08 Dec 2009 Beta regression Estimation of a beta regression model (logit link function). The parameters are estimated with ML. Author: Willem-Jan de Goeij beta, regression, model, generalized, linear, maximum 15 0
16 Aug 2009 Student's t regression Estimates a Student's t regression model Author: Willem-Jan de Goeij studentt, student, t, regression, model, generalized 19 1
  • 5.0
5.0 | 1 rating
Comments and Ratings by Willem-Jan View all
Updated File Comments Rating
22 Oct 2012 Shapiro-Wilk and Shapiro-Francia normality tests. Shapiro-Wilk & Shapiro-Francia parametric hypothesis test of composite normality. Author: Ahmed Ben Saïda

Can you explain the following?
x is normally distributed.
If I perform a 2 tailed test, your function rejects the null hypothesis.

x = norminv((1:9)/10);
[h,p,w]=swtest(x,0.05,0)
h =
1
p =
0.0028
w =
0.9925

04 Feb 2012 inpaint_nans Interpolates (& extrapolates) NaN elements in a 2d array. Author: John D'Errico

Actually, I made a mistake. I've mixed two examples I was looking at.

The average of 2 and 0 is of course 1.
So I would pick 1 as a solution.

What I see in your code is that the problem above is solved by the following set of equations:

2 x = 1+3 = 4
x = 10+6+4 - 4*5 = 0

The least squares solution for this set of equations, is indeed x=1.6.
But why is the first equation multiplied by 2?

If you keep the factor for x in both equations equal to one, you get

x = (1+3)/2 = 2
x = 10+6+4-4*5 = 0

This gives x=1 as least squares solution.

03 Feb 2012 inpaint_nans Interpolates (& extrapolates) NaN elements in a 2d array. Author: John D'Errico

A =
1 4 7
NaN 5 10
3 6 11

Lets call the NaN value x.

Your function gives x=1.6.
I think x=1.5 is more logical, how about you?

If we look at the first column, we have x=0.5*(1+3). ==> x=2
If we look at the 5, we have
5=0.25*(x+10+4+6). ==> x=0

If we average these values, we have x=1.5.

08 Aug 2011 Shuffle Random permutation of array elements, C-Mex: much faster than RANDPERM Author: Jan Simon

Dear Jan,

I've tested calling your function from inside a parfor loop. The results show that repetitions occur.

>> matlabpool local
Starting matlabpool using the 'local' configuration ... connected to 3 labs.
>> shuffleparfor
output =

8 1 3 7 5 2 6 4
2 3 8 7 1 5 6 4
8 1 3 7 5 2 6 4
2 3 8 7 1 5 6 4
8 1 3 7 5 2 6 4
2 3 8 7 1 5 6 4
4 8 7 3 6 2 5 1
4 8 7 3 6 2 5 1
4 8 7 3 6 2 5 1

output = zeros(9,8);
Shuffle(1234567890, 'seed');

parfor ii=1:9
output(ii, :) = Shuffle(8, 'index');
end

03 Aug 2011 Shuffle Random permutation of array elements, C-Mex: much faster than RANDPERM Author: Jan Simon

Hi Jan,

In reply to your first comment from 31/07:
Setting the seed to randi([0, 2^32-1], 1, 4) in every parfor loop would be advisible, I guess?
On my computer 2^32-1+d becomes zero when casted to uint32_T, with d >= 0.999999762.

I'm setting seeds because I also want to be able to reproduce the results.

Comments and Ratings on Willem-Jan's Files View all
Updated File Comment by Comments Rating
23 Oct 2013 PARFOR Progress Monitor v2 This version is an update of the progress monitor made by Edric Ellis. Author: Willem-Jan de Goeij Jerome

Do you have any idea why this clear all global variables?
Without this waitbar my globals stay fine.

29 Aug 2013 PARFOR Progress Monitor v2 This version is an update of the progress monitor made by Edric Ellis. Author: Willem-Jan de Goeij Arwel

For very large arrays the crash still happens - i.e for an array of 2.5e10 elements (as an array of singles), this still throws "java.lang.IllegalArgumantException: invalid range properties", irrespective of the value of progressStepSize.

03 Dec 2012 PARFOR Progress Monitor v2 This version is an update of the progress monitor made by Edric Ellis. Author: Willem-Jan de Goeij Montorsi, Francesco

very useful in my parallel simulations! by the way, whenever you run parallel simulations where you don't know in advance the number of cycles of the parfor, it would be very useful to have a function to update the total number of cycles given at construction time...

24 Jan 2012 PARFOR Progress Monitor v2 This version is an update of the progress monitor made by Edric Ellis. Author: Willem-Jan de Goeij Edmondo

it could be great to have also an estimated remaining time

20 Jan 2012 PARFOR Progress Monitor v2 This version is an update of the progress monitor made by Edric Ellis. Author: Willem-Jan de Goeij Thomas, Adam

Works perfectly in my application, which is a deployed parallel function.

When deploying a function using this progress monitor note that "pctRunOnAll javaaddpath java" will have no effect, you need to embed the folder containing the java classes into the CTF archive when you compile... then it runs like a charm :)

Thanks

Contact us