# Jeff Miller

**Statistics**

RANK**86**

of 246,934

REPUTATION**1,368**

CONTRIBUTIONS

**14** Questions

**554** Answers

ANSWER ACCEPTANCE **
85.71%
**

VOTES RECEIVED**141**

RANK**5,692** of 17,223

REPUTATION**178**

AVERAGE RATING**4.70**

CONTRIBUTIONS**8** Files

DOWNLOADS **11**

ALL TIME DOWNLOADS**1648**

CONTRIBUTIONS**0 Posts**

CONTRIBUTIONS**0** Public Channels

AVERAGE RATING

CONTRIBUTIONS**0** Highlights

AVERAGE NO. OF LIKES

**Content Feed**

understanding Linear regression model equation after simplifying.

I think the notation y~ 1+ x1*x2 is shorthand for y~ 1*b0+ b1*x1 + b2*x2 + b3*x1*x2. The four b values are the four estimates...

14 days ago | 0

| accepted

Is there a way to normalizing a probability density function / keep the density values between 0 and 1?

Dividing by the max density value would have the unfortunate ramification that the resulting function would no longer be a densi...

14 days ago | 0

| accepted

plot of log likelihood values

I think this might be what you are after: x = [2 4 3 1 6 2 2.5 7 10]; % data ShapeCandidates = 1:10; Scale = 6; loglik...

21 days ago | 0

How can I see specific likelihood values with MLE (or possibly another function?)

Once you have the parameter estimates (e.g., from mle) you can compute the likelihood value by calling your custom density funct...

29 days ago | 0

Facing an error while using 'anovan' function

anovan needs the original data (y1 and y2) to compute an error term, but you are only giving it the averages (ybar). This is why...

29 days ago | 0

| accepted

Generating random variable from inverse gaussian dist.

Do you really need to write the random generator yourself? This one is already available: >> pd = makedist('inversegaussian') ...

1 month ago | 0

How do I get MATLAB to display the d' when my hit rate = 1 using the norminv function?

No, d' cannot be computed when an observed hit_rate is 1 or fa_rate is zero. One of the most common work-arounds for this probl...

1 month ago | 0

How to delete some of anonymous variables randomly?

One way to approach the problem is to write a function f(x) where x is a vector with 47 positions. When you call f, put numeric...

1 month ago | 0

| accepted

Random Laplacian distribution in a specified interval

Here is a quick and dirty way to do this with Cupid: location = 0; % set whatever Laplacian parameters you want. scale = 1; ...

2 months ago | 0

| accepted

kstest for uniform distribution

you can test for the fit of a discrete distribution, including uniform, with chi2gof. One of the examples (about 1/3 of the way...

2 months ago | 1

| accepted

How to remove the outliers in distribution fitter

If you click on the "Exclude" button in the picture you posted, you will get the option to set lower and upper limits to exclude...

2 months ago | 1

Bootstrap multi linear model

The problem is that fitlm produces a complicated output structure and bootstrp doesn't know what to do with it. x = (1:1000)'; ...

2 months ago | 0

| accepted

How Can I simulate Data for Least-Squares Variance Component Estimation (LS-VCE) for Multiple regression?

You could use mvnrnd to generate random values of x1-x4. Is that what you need?

2 months ago | 0

How to group data for MANOVA1 analysis

84*84 is not going to work, sorry. For MANOVA, you need to have more cases in each group than the number of variables. Some pos...

2 months ago | 0

How to test if values are close togather using p-value (ttest)

A small p value means that the data are incompatible with some underlying assumption/model. For example, you could potentially ...

3 months ago | 0

Finding the minimum of a complicated function with variables and parameters

If I understand the problem correctly, then this should be close to what you want: function bestEsts = MyOuterFn(x,y) MyErrFn ...

3 months ago | 0

How to get correct p-values for the Steel-Dwass test?

The problem is with the 'inf' value given as the degrees of freedom in this line p_steelDwass(i, 1) = 1 - tcdf(t_steelDwass...

3 months ago | 0

How to calculate for significant difference between Cohen's Kappa values?

As I understand it, the fundamental question is whether tests A & B agree better than tests A & C, beyond a minor improvement th...

3 months ago | 0

fsolve doesn't find a solution when the value of the objective function is obtained through a Montecarlo simulation

I don't know if it can be done in the context of your simulation, but sometimes it is helpful to pick the "random" variables sys...

3 months ago | 0

Question

speed up integrating the same function over many overlapping intervals?

I need to compute many different integrals of the same function fn, essentially like this: Xint = size(X); for i=1:numel(X) ...

3 months ago | 2 answers | 0

### 2

answersDifferences between functions ezfit and Matlab's gamfit in fitting data to the gamma distribution?

There are two somewhat different parameterizations of the gamma distribution. In one (which I think is used by gamfit), the mean...

3 months ago | 0

How do I calculate the R^2?

Either of these: Rsqr = corr(x(:,2),y)^2; RsqrAlso = 1-(std(resids)/std(y))^2;

3 months ago | 1

I would like to find elements of a matrix that fulfill a condition but after every loop the condition is changing.

I think this will almost do what you want. edges = 0:0.05:30; [~,~,bin_numbers] = histcounts(diameters,edges); The output bin...

3 months ago | 0

| accepted

Generating a categorical array with group names

Is this what you are after? part1 = repmat({'Name1'},1055,1); part2 = repmat({'Name2'},1641,1); part3 = repmat({'Name3'},2055...

3 months ago | 0

| accepted

Forming a CDF by adding different types of PDF and generating random number

It seems like there are several different questions here. If you want to generate random numbers from a custom CDF, one simple ...

3 months ago | 2

| accepted

Calculating pdf and cdf for data extracted as a table

This command pdfwindspeed = pdf('Normal',A) returns the pdf of the A values within the standard normal with mu=0 and sigma=1. ...

3 months ago | 0

How to use one data set to find the mean of another?

NON1 = NON == 1; % this makes a boolean vector that is true for the positions where NON==1 mean(PLCm(NON1)); % this computes ...

3 months ago | 1

Question

Request for help computing convolution of random variables via fft

Inspired by the answers of @David Goodmanson to this earlier question, I have been trying to understand exactly how in general t...

3 months ago | 2 answers | 0

### 2

answersrun glmfit , why get Warning: X is ill conditioned, or the model is overparameterized

The problem is that the x1-x5 values sum to 24 in every row. This means that any one of the x's can be predicted from the rest ...

3 months ago | 0

| accepted

CDF of VOn Mises distribution

The Von Mises distribution is included in Cupid. You could use it like this to calculate CDF values: >> location = 10; >> conc...

3 months ago | 0