Community Profile

# Jeff Miller

##### Last seen: Today
528 total contributions since 2013

View details...

Contributions in
View by

Problems using `fitdist` to Rician distribution
I'm not sure why you are getting that error message, but Rician(6.5538e-05,0.11714) looks quite good: In case you are curious...

6 days ago | 0

OOP: Reference (subclass) enumeration names from static superclass method. Possible?
If doWork really just needs the names, maybe pass allnames as an argument of doWork, something like this: function resu...

7 days ago | 2

Two-Way repeated Measures Anova: Pairwise comparison of multiple factor levels using Multcompar
Go back to the original data and compute the means for each group of interest, so you can compare which group mean is higher.

11 days ago | 0

Fitting Gaussian keeping max amplitude same
One option is to fit a smoothed curve and then multiply the (smoothed, predicted) height by whatever constant is needed to make ...

12 days ago | 1

Develop uncertainty estimates using Monte-Carlo simulation, bootstrap-resampling
You seem to be asking how to select a random subset of the years for each iteration of the simulation. If that is right, then y...

14 days ago | 0

| accepted

How to interpret Repeated Measures Anova Example Table
(Intercept):Measurements is the "main effect" of measurement type--a factor with 4 levels in this example (i.e. length and width...

14 days ago | 1

| accepted

PDF and CDF of rician distribution simulation source
OK, here's a more complete answer: pd = makedist('Rician','s',0,'sigma',2); % Change to your required parameters N = 1000; %...

16 days ago | 0

PDF and CDF of rician distribution simulation source
Why can't you use this: pd = makedist('Rician')

19 days ago | 0

how to refer variable from string name
This is pretty clumsy, but I think something like this will do what you want: varnames = {'m1_n', 'm3_n'}; % names of the vari...

27 days ago | 0

| accepted

generating correlated random variables
Generate C1-C6 separately (independently) from whatever lognormal you want. Then generate C7-C8 as a pair with the desired corr...

1 month ago | 0

| accepted

Stop signal task - pre error speeding
I think it would be cleanest to start by making a vector to indicate which trials precede unsuccessful stops, which can be done ...

1 month ago | 0

| accepted

How to make Anova table for data from a factorial experiment
x = readtable('mack_data.txt'); x.Properties.VariableNames = {'o1','o2','o3','o4','o5','m1','m2','m3','m4','m5'}; withinDesign...

2 months ago | 1

| accepted

How to create a cumulative frequency array
This is inelegant, but it works: cumFreq = zeros(size(V)); for i=1:numel(V) cumFreq(i) = sum(Z<V(i)); end

2 months ago | 1

| accepted

Define a density function and draw N samples of it
I don't think that density function is right--seems like you at least need to include the mixture probabilities in there--e.g., ...

2 months ago | 0

Function to create struct with parameters
Maybe you would like something like this? month = 1; year = 1951; day = 11; inp = CreateStruct(month,day,year) function o...

2 months ago | 0

| accepted

How do I declare correlation between a random variable with beta distribution and a random variable with normal distribution?
Something like this would work, though you would have to adjust rho by trial and error to get the final corr(x1,x2) value that y...

2 months ago | 0

| accepted

Extract columns into new variables efficently
A= [1 1 2 2 3 3 1 1 2 2 3 3; 1 1 2 2 3 3 1 1 2 2 3 3]; B = [10 20 30 40 50 60 70 80 90 100 110 120; 1 2 3 4 5 6 7 8 9 10 11 12]...

2 months ago | 0

| accepted

Random numbers with specific sum and varying maximum
How about first picking the maximum "thismax" randomly from whatever distribution you want, and then using randfixedsum to gener...

2 months ago | 0

Using ranova and multcompare, ranova doen't show significance (p=0.33) but multcompare shows significance between groups (p<0.001)
Note that your anova table doesn't have enough lines--there is no main effect of IndVar, just the IndVar*Time interaction (p=0.3...

2 months ago | 0

| accepted

How may I show a mixture distribution?
Here is an example of plotting a mixture that produces a graph similar to yours example: pd1 = makedist('normal',1,1); % proba...

2 months ago | 1

| accepted

How to design a Class to apply operations for a collection of data?
This is perhaps a bit subjective, but personally I would strongly prefer the solution with classes that process only a single au...

2 months ago | 0

| accepted

histogram normalization not working
The histogram is not wrong, because pdf values can be more than 1. pdf values are normalized such that the total area under the...

2 months ago | 0

Equal Bin Sizes in Histogram Probability Distribution Function (PDF)
The histogram function has a parameter 'BinEdges' that will allow you to specify any bins you like. So, make whatever bin edges...

2 months ago | 1

Minimization of Chi-Square Statistic using the Simplex Algorithm
One way to proceed is to write a function which accepts a vector of values for your three parameters, and returns the chi-square...

2 months ago | 0

How to generate correlated samples from a single distribution? Each sample should be correlated with the previous.
If you don't like the sorting solution, then I think you must consider the joint and conditional distribution of two successive ...

3 months ago | 0

| accepted

Inverse Cumulative distribution function
Suppose you have a function called thisCDF which computes the CDF for your distribution. Then define function x = inverseF(p) ...

3 months ago | 0

| accepted

How to generate correlated samples from a single distribution? Each sample should be correlated with the previous.
If you don't care about the order of the data points, one option that seems to fulfill your stated requirements is the following...

3 months ago | 0

I want this for loop to dispaly the smallest calculated value of PE. What would be the best way of going about this?
% something like this: minPE = realmax; for t1 = 0:1:90,t2 = 0:1:90; t1; t2; h1 = -b1*cosd(t1); ...

3 months ago | 0

Normally distributed Random numbers generator issue
As the other posters have said, you aren't necessarily going to get an exactly normal-looking distribution if you generate the n...

3 months ago | 0