Community Profile Jeff Miller

Last seen: Today Active since 2013

Statistics

All
• • • • • • • • • • • Content Feed

View by

How to turn a function handle with 3 inputs (1 variable and 2 parameters) and then assign the inputs parameters and get a function handle with one input? ?
Not sure I understand the situation, but this sounds to me like a case where an OO approach could be very helpful. Consider: cl...

11 days ago | 0

Repeated measures ANOVA Matlab (Statistic)
% Note the transpose operator: single-column vectors are more convenient. A = [10 11 13 16 15 18 17 16 11 12]'; B = [20 22 21 ...

13 days ago | 0

| accepted

How to fit an uniform distribution to a histogram?
histogram(data, numBars); avgHeight = numel(data)/numBars; hold on; plot([min(data), max(data)],[avgHeight, avgHeight])

14 days ago | 1

Spearman correlation fails in ranking strings
The first two inputs to the corr command have to be the numerical values that are to be correlated, so you will want something l...

15 days ago | 1

| accepted

I want to compute mean median and mode for the distribution from its c.d.f.
The median is the value at p=0.50, which you already seem to have in your table. You can't really compute the mean and mode fr...

21 days ago | 0

How to calculate StdErr for intercepts and slopes in a multivariate ANCOVA (function ANOVAN() with continuous independent variable)?
If you set up dummy variables to code your categorical independent variables, you can then do analysis of covariance with fitlm....

22 days ago | 2

| accepted

Is it appropriate to use ranksum to compare statistical means for time series data?
The ranksum test ignores the sequential order of scores, so for example ranksum(x,y) will give you the same answer no matter how...

1 month ago | 0

pca with more than one grouping variable
gscatter will accept multiple grouping variables, in your case organism & dead/alive. There is an example in the gscatter docu...

1 month ago | 0

| accepted

Distribute a sum using binomial distribution
If I understand correctly, it is fine to use binornd with n = row total. It sounds like that is what was done in the paper. j...

2 months ago | 0

Why is fitlm affected by variable scale?
Your understanding is correct for linear regression but your model is nonlinear because of the interaction terms. Consider: zX ...

2 months ago | 0

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...

2 months 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...

2 months 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...

3 months 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...

3 months 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...

3 months 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') ...

3 months 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...

3 months 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...

3 months 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; ...

3 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...

4 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...

4 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)'; ...

4 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?

4 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...

4 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 ...

4 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 ...

4 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...

4 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...

4 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...

5 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) ...

5 months ago | 2 answers | 0