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

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

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

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

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

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

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

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

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

8 mesi 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') ...

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

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

9 mesi 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; ...

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

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

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

9 mesi 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?

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

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

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

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

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

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

10 mesi 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) ...

10 mesi ago | 2 answers | 0

### 2

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

10 mesi 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;

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

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

10 mesi ago | 0

| accepted