Answered

log likelihood output from the distribution fitter app

Just be a bit careful here... I think those log likelihood values are likelihoods under the maximum likelihood parameter estima...

log likelihood output from the distribution fitter app

Just be a bit careful here... I think those log likelihood values are likelihoods under the maximum likelihood parameter estima...

6 days ago | 0

Answered

Finding P-Values with the Null Hypothesis Being a Coefficient Equal to 0

The easiest way is with the corr function, because the p value for the slope is the same as the p value for the correlation coef...

Finding P-Values with the Null Hypothesis Being a Coefficient Equal to 0

The easiest way is with the corr function, because the p value for the slope is the same as the p value for the correlation coef...

6 days ago | 0

Answered

How do I set a seed to generate different random initial numbers and storing them

Maybe the FileExchange contribution RNGMgr will be helpful.

How do I set a seed to generate different random initial numbers and storing them

Maybe the FileExchange contribution RNGMgr will be helpful.

10 days ago | 0

| accepted

Answered

How to accurately process an Anderson Darling test?

Looks to me like the same problem as in your previous question. MATLAB's 2-parameter Weibull starts from a minimum of 0, but Mi...

How to accurately process an Anderson Darling test?

Looks to me like the same problem as in your previous question. MATLAB's 2-parameter Weibull starts from a minimum of 0, but Mi...

13 days ago | 2

| accepted

Answered

Where is the coefficient for the reference condition when using 'fitlm' to perform ANOVA with intercept omitted?

Why not just use this? for iGrp=1:5 stderr = std(y(x==iGrp)) / sqrt( sum(x==iGrp) ) end

Where is the coefficient for the reference condition when using 'fitlm' to perform ANOVA with intercept omitted?

Why not just use this? for iGrp=1:5 stderr = std(y(x==iGrp)) / sqrt( sum(x==iGrp) ) end

16 days ago | 0

| accepted

Answered

How to solve for 3 parameters for a Weibull distribution

Here is what I get using Cupid: Build_Array_Final{1}=[214.42,214.99,215.35,215.98,221.08,222.54,214.78,215,215.32,215.75,217.06...

How to solve for 3 parameters for a Weibull distribution

Here is what I get using Cupid: Build_Array_Final{1}=[214.42,214.99,215.35,215.98,221.08,222.54,214.78,215,215.32,215.75,217.06...

17 days ago | 0

| accepted

Answered

How to find the minimum of the function?

fminsearch is a good place to start

How to find the minimum of the function?

fminsearch is a good place to start

17 days ago | 0

Answered

Distribution that works with positive and negative skewness

The skew normal is in Cupid. The last parameter controls the amount and direction of skew. E.g., posSkew = SkewNor(10,1,2); p...

Distribution that works with positive and negative skewness

The skew normal is in Cupid. The last parameter controls the amount and direction of skew. E.g., posSkew = SkewNor(10,1,2); p...

25 days ago | 0

| accepted

Answered

how can I do a double summation with dependent indices?

Sorry, I am having trouble relating your code to your equation, and I would have no idea how to do this with syms. But maybe go...

how can I do a double summation with dependent indices?

Sorry, I am having trouble relating your code to your equation, and I would have no idea how to do this with syms. But maybe go...

1 month ago | 0

| accepted

Answered

Random samples from combination of poisson distributions

One easy way to generate a random score from a poisson mixture distribution is to randomly select the particular poisson first a...

Random samples from combination of poisson distributions

One easy way to generate a random score from a poisson mixture distribution is to randomly select the particular poisson first a...

1 month ago | 0

| accepted

Answered

How can I calculate the 95% confidence interval for spearman correlation values?

I don't think there are formulas for that but you can estimate it with bootstrapping: mycorr = @(x1,x2) corr(x1,x2,'Type','Spea...

How can I calculate the 95% confidence interval for spearman correlation values?

I don't think there are formulas for that but you can estimate it with bootstrapping: mycorr = @(x1,x2) corr(x1,x2,'Type','Spea...

1 month ago | 0

| accepted

Answered

Finding the t values of the Student's t test (t1-α,ν)

critt = tinv(1-alpha,v)

Finding the t values of the Student's t test (t1-α,ν)

critt = tinv(1-alpha,v)

1 month ago | 0

| accepted

Answered

Draw PDF and CDF of random values of distance from a 10X10 matrix

I'm guessing you want the empirical PDF and CDF functions (that will change slightly every time you rerun the code, due to rando...

Draw PDF and CDF of random values of distance from a 10X10 matrix

I'm guessing you want the empirical PDF and CDF functions (that will change slightly every time you rerun the code, due to rando...

1 month ago | 0

Answered

How to find 2d spearman and kendall correlations using corr2?

corr2 is doing the same computation as corr, except that it allows 2d inputs rather than 1d. So, you can get your Spearman or K...

How to find 2d spearman and kendall correlations using corr2?

corr2 is doing the same computation as corr, except that it allows 2d inputs rather than 1d. So, you can get your Spearman or K...

1 month ago | 0

| accepted

Answered

calculate empirical distribution function and interpolation

You want to use the regress function, something like this: X = [ones(size(temp) temp)]; % temp is a column vector of temperatu...

calculate empirical distribution function and interpolation

You want to use the regress function, something like this: X = [ones(size(temp) temp)]; % temp is a column vector of temperatu...

1 month ago | 0

Answered

How to find P value in chi2inv

Not sure what you are asking, but maybe p = chi2cdf(3.7920,df,'upper') where df is the degrees of freedom for your computed ch...

How to find P value in chi2inv

Not sure what you are asking, but maybe p = chi2cdf(3.7920,df,'upper') where df is the degrees of freedom for your computed ch...

1 month ago | 0

Answered

Updating multiple variables randomly with a constraint in a while cycle

You might be able to do what you want with nested while loops, something like this. Pretty tedious though. x = zeros(6,1); st...

Updating multiple variables randomly with a constraint in a while cycle

You might be able to do what you want with nested while loops, something like this. Pretty tedious though. x = zeros(6,1); st...

1 month ago | 0

Answered

How to test for a main effect of the between subjects-variable with fitrm and ranova?

% Define the within-Ss factor and give it a name. The 1 & 2 refer to Awake / Asleep WithinStructure = table([1 2]','VariableNa...

How to test for a main effect of the between subjects-variable with fitrm and ranova?

% Define the within-Ss factor and give it a name. The 1 & 2 refer to Awake / Asleep WithinStructure = table([1 2]','VariableNa...

1 month ago | 0

Answered

fitgmdist for loop; nested for loop, combine array of different sizes

I don't really follow what you are trying to do, but maybe some of these comments will help a bit: I. Use a cell array to outpu...

fitgmdist for loop; nested for loop, combine array of different sizes

I don't really follow what you are trying to do, but maybe some of these comments will help a bit: I. Use a cell array to outpu...

2 months ago | 1

| accepted

Answered

How do I fit multiple Gaussians to my data?

From the new graph you've posted, it looks like the values in AbsY are not the type of data expected by fitgmdist, and I suspect...

How do I fit multiple Gaussians to my data?

From the new graph you've posted, it looks like the values in AbsY are not the type of data expected by fitgmdist, and I suspect...

2 months ago | 1

Answered

How to compare appropriateness of fitted Weibull and lognormal models via the likelihood ratio test

Unfortunately the likelihood ratio test cannot be used this way. It is only for comparing nested models (i.e., one model has to...

How to compare appropriateness of fitted Weibull and lognormal models via the likelihood ratio test

Unfortunately the likelihood ratio test cannot be used this way. It is only for comparing nested models (i.e., one model has to...

2 months ago | 0

| accepted

Answered

how can i check the geometric distribution

Look at the Wikipedia page on this distribution and it will tell you how to estimate the parameter from your observations (it de...

how can i check the geometric distribution

Look at the Wikipedia page on this distribution and it will tell you how to estimate the parameter from your observations (it de...

2 months ago | 1

Answered

Use 'trandn' function to draw random numbers from a lognormal distribution

An easy way to get random numbers from a lognormal distribution is to get random numbers from a normal distribution and then exp...

Use 'trandn' function to draw random numbers from a lognormal distribution

An easy way to get random numbers from a lognormal distribution is to get random numbers from a normal distribution and then exp...

2 months ago | 0

Answered

Ploting the result of a fitgmdist (of 1D data)

Say you want to plot the distribution from -15 to 15 in steps of 0.05. This would be the code: z = (-15:0.05:15)'; % Note th...

Ploting the result of a fitgmdist (of 1D data)

Say you want to plot the distribution from -15 to 15 in steps of 0.05. This would be the code: z = (-15:0.05:15)'; % Note th...

2 months ago | 0

| accepted

Answered

Ranking Variables by Value

You might find the second output of the sort command useful. E.g. fruitVals = [p_a p_b p_o p_p]; [~, idxOrder] = sort(fruitVa...

Ranking Variables by Value

You might find the second output of the sort command useful. E.g. fruitVals = [p_a p_b p_o p_p]; [~, idxOrder] = sort(fruitVa...

2 months ago | 0

| accepted

Answered

help random numbers generator

Good. Here it is as an 'official' answer: pd = makedist('Normal'); t = truncate(pd,-0.5,0.5); r = random(t,10,1);

help random numbers generator

Good. Here it is as an 'official' answer: pd = makedist('Normal'); t = truncate(pd,-0.5,0.5); r = random(t,10,1);

2 months ago | 0

| accepted

Answered

How to generate random numbers from a scaled beta PDF?

Generate regular beta random numbers and then rescale them. I guess you want to generate randomly from the regular beta(3,3) di...

How to generate random numbers from a scaled beta PDF?

Generate regular beta random numbers and then rescale them. I guess you want to generate randomly from the regular beta(3,3) di...

2 months ago | 0

| accepted

Answered

Why does stepwiselm function output NaN with model interactions?

Multivariate techniques like stepwise regression require more cases (rows) than variables--in practice, you need at least 5-10 t...

Why does stepwiselm function output NaN with model interactions?

Multivariate techniques like stepwise regression require more cases (rows) than variables--in practice, you need at least 5-10 t...

2 months ago | 0

Answered

Concatenate string to string array in efficient way

strMsgArray = strcat(strStart,strMsgArray,strEnd)

Concatenate string to string array in efficient way

strMsgArray = strcat(strStart,strMsgArray,strEnd)

2 months ago | 0

| accepted

Answered

question related to Statistics in matlab

r=poissrnd(9,200,1); will give you 200 random numbers from a poisson distribution with lambda=9. Lots of those numbers will be...

question related to Statistics in matlab

r=poissrnd(9,200,1); will give you 200 random numbers from a poisson distribution with lambda=9. Lots of those numbers will be...

2 months ago | 0