Mrthfx

How do you funnel food off a cutting board?

What makes the Forgotten Realms "forgotten"?

En Passant For Beginners

Do authors have to be politically correct in article-writing?

Avoiding morning and evening handshakes

Recrystallisation of dibenzylideneacetone

It took me a lot of time to make this, pls like. (YouTube Comments #1)

Does my logo design convey the right feelings for a University Student's Council?

just upgraded iMac late 2015 ram from 12 GB to 28 GB and it became too slow to use

Can we use the stored gravitational potential energy of a building to produce power?

Why would the Pakistan airspace closure cancel flights not headed to Pakistan itself?

Program that converts a number to a letter of the alphabet

What incentive is there for baking on the test chain during voting period?

Tikzing a circled star

If I delete my router's history can my ISP still provide it to my parents?

What to do if authors don't respond to my serious concerns about their paper?

What's a good word to describe a public place that looks like it wouldn't be rough?

Using loops to create tables

Finding radius of circle

Does Windows 10's telemetry include sending *.doc files if Word crashed?

How to acknowledge an embarrassing job interview, now that I work directly with the interviewer?

How to generate a matrix with certain conditions

How did the original light saber work?

Why is working on the same position for more than 15 years not a red flag?

How to generate a matrix with certain conditions

Standardizing a coset table via matrix manipulationEigenvalues and Determinant of a large matrixHow to generate a random matrix with specific parameters?Generate Non-Singular Matrix of $ntimes n$ dimensionHow to create a matrix with some conditions?Problems with RandomChoicePicking first row of a matrix with a positive first elementSimplest way to construct a matrix its elements are defined by known functionsmatrix with chosen elements distributed in a random positionUsing Table to build a matrix

1

$begingroup$

I want to generate a $n times n$ matrix.

1. I want the diagonal entries to be all 0


2. I want a random choice of matrix elements with 0 or 1.


3. The probability of having a 1 as a matrix element is $1/m$ and the probability of having a 0 as a matrix element is $1-1/m$.

I used the following command but it is wrong.

A[n_, m_] :=Table[If[i == j, 0,RandomVariate[BernoulliDistribution[m],{n,n}]]]

And I tried to test this command with n=4, m=0.4 but it didn't work.

Could anyone kindly tell me how to do this please?
Thank you!

matrix probability-or-statistics random sampling

share|improve this question

edited 1 hour ago

J. M. is computer-less♦

97k10303463

asked 6 hours ago

tiffanytiffany

61

New contributor

tiffany is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  For a start, your code seems to set the diagonal elements to 1 rather than zero.
  $endgroup$
  – MarcoB
  6 hours ago
  

  
  
  
  

add a comment | 

1

$begingroup$

I want to generate a $n times n$ matrix.

1. I want the diagonal entries to be all 0


2. I want a random choice of matrix elements with 0 or 1.


3. The probability of having a 1 as a matrix element is $1/m$ and the probability of having a 0 as a matrix element is $1-1/m$.

I used the following command but it is wrong.

A[n_, m_] :=Table[If[i == j, 0,RandomVariate[BernoulliDistribution[m],{n,n}]]]

And I tried to test this command with n=4, m=0.4 but it didn't work.

Could anyone kindly tell me how to do this please?
Thank you!

matrix probability-or-statistics random sampling

share|improve this question

edited 1 hour ago

J. M. is computer-less♦

97k10303463

asked 6 hours ago

tiffanytiffany

61

New contributor

tiffany is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  For a start, your code seems to set the diagonal elements to 1 rather than zero.
  $endgroup$
  – MarcoB
  6 hours ago
  

  
  
  
  

add a comment | 

$begingroup$

I want to generate a $n times n$ matrix.

1. I want the diagonal entries to be all 0


2. I want a random choice of matrix elements with 0 or 1.


3. The probability of having a 1 as a matrix element is $1/m$ and the probability of having a 0 as a matrix element is $1-1/m$.

I used the following command but it is wrong.

A[n_, m_] :=Table[If[i == j, 0,RandomVariate[BernoulliDistribution[m],{n,n}]]]

And I tried to test this command with n=4, m=0.4 but it didn't work.

Could anyone kindly tell me how to do this please?
Thank you!

matrix probability-or-statistics random sampling

share|improve this question

edited 1 hour ago

J. M. is computer-less♦

97k10303463

asked 6 hours ago

tiffanytiffany

61

New contributor

tiffany is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

$endgroup$

A[n_, m_] :=Table[If[i == j, 0,RandomVariate[BernoulliDistribution[m],{n,n}]]]

share|improve this question

edited 1 hour ago

J. M. is computer-less♦

97k10303463

asked 6 hours ago

tiffanytiffany

61

New contributor

tiffany is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

share|improve this question

edited 1 hour ago

J. M. is computer-less♦

97k10303463

asked 6 hours ago

tiffanytiffany

61

New contributor

tiffany is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

edited 1 hour ago

J. M. is computer-less♦

97k10303463

edited 1 hour ago

J. M. is computer-less♦

97k10303463

asked 6 hours ago

tiffanytiffany

61

New contributor

tiffany is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

asked 6 hours ago

tiffanytiffany

61

3 Answers
3

active

oldest

votes

3

$begingroup$

Binary random variables are often modeled using the BernoulliDistribution. You can use the function RandomVariate to get a matrix of such variables.

mat = RandomVariate[BernoulliDistribution[0.9], {5, 5}]; 

mat - DiagonalMatrix[Diagonal[mat]]

%//MatrixForm

You can change the 0.9 to any value (this is your m). The second line sets all the diagonal elements to zero.

share|improve this answer

answered 6 hours ago

bill sbill s

53.6k376153

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Thank you bill s How will the command be if I do not restrict what n and m be? For example, I would like to generate a set of commands which I can replace m and n easily by any number.
  $endgroup$
  – tiffany
  5 hours ago
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  @tiffany, just do With[{n = 8, m = 3}, (# - DiagonalMatrix[Diagonal[#]]) &[RandomVariate[BernoulliDistribution[1/m], {n, n}]]].
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

3

$begingroup$

You can use weight option in RandomChoice

n = 5;

m = 2;

mat = RandomChoice[{1/m, 1 - 1/m} -> {1, 0}, {n, n}];

(mat - DiagonalMatrix[Diagonal@mat]) // MatrixForm

$left(
begin{array}{ccccc}
0 & 1 & 1 & 0 & 1 \
1 & 0 & 1 & 0 & 1 \
1 & 0 & 0 & 0 & 1 \
0 & 0 & 1 & 0 & 0 \
0 & 0 & 0 & 1 & 0 \
end{array}
right)$

share|improve this answer

answered 4 hours ago

Okkes DulgerciOkkes Dulgerci

5,0891917

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Altho using the BernoulliDistribution[] is best, this is likely to be more easy to read for someone who is not accustomed to discrete probability distributions.
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

0

$begingroup$

m = 2;

n = 5;

rnd[x_] := If[RandomReal[{0, 1}] < 1/m, 1, 0];

t = Table[rnd[x_]*(1 - KroneckerDelta[i, j]), {i, 1, n}, {j, 1, n}];

t // MatrixForm

share|improve this answer

answered 6 hours ago

Vsevolod A.Vsevolod A.

478211

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Thank you Vsevolod A. How will the command be if I do not restrict what n and m be?
  $endgroup$
  – tiffany
  6 hours ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @tiffany you set m and n in the first two lines...
  $endgroup$
  – Vsevolod A.
  4 hours ago
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Just rnd := If[RandomReal[{0, 1}] < 1/m, 1, 0]; will do, since the function never uses its argument.
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

Your Answer

StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "387"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});

tiffany is a new contributor. Be nice, and check out our Code of Conduct.

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f192468%2fhow-to-generate-a-matrix-with-certain-conditions%23new-answer', 'question_page');
}
);

Post as a guest

Name

Email

Required, but never shown

3

$begingroup$

Binary random variables are often modeled using the BernoulliDistribution. You can use the function RandomVariate to get a matrix of such variables.

mat = RandomVariate[BernoulliDistribution[0.9], {5, 5}]; 

mat - DiagonalMatrix[Diagonal[mat]]

%//MatrixForm

You can change the 0.9 to any value (this is your m). The second line sets all the diagonal elements to zero.

share|improve this answer

answered 6 hours ago

bill sbill s

53.6k376153

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Thank you bill s How will the command be if I do not restrict what n and m be? For example, I would like to generate a set of commands which I can replace m and n easily by any number.
  $endgroup$
  – tiffany
  5 hours ago
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  @tiffany, just do With[{n = 8, m = 3}, (# - DiagonalMatrix[Diagonal[#]]) &[RandomVariate[BernoulliDistribution[1/m], {n, n}]]].
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

3

$begingroup$

Binary random variables are often modeled using the BernoulliDistribution. You can use the function RandomVariate to get a matrix of such variables.

mat = RandomVariate[BernoulliDistribution[0.9], {5, 5}]; 

mat - DiagonalMatrix[Diagonal[mat]]

%//MatrixForm

You can change the 0.9 to any value (this is your m). The second line sets all the diagonal elements to zero.

share|improve this answer

answered 6 hours ago

bill sbill s

53.6k376153

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Thank you bill s How will the command be if I do not restrict what n and m be? For example, I would like to generate a set of commands which I can replace m and n easily by any number.
  $endgroup$
  – tiffany
  5 hours ago
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  @tiffany, just do With[{n = 8, m = 3}, (# - DiagonalMatrix[Diagonal[#]]) &[RandomVariate[BernoulliDistribution[1/m], {n, n}]]].
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

$begingroup$

Binary random variables are often modeled using the BernoulliDistribution. You can use the function RandomVariate to get a matrix of such variables.

mat = RandomVariate[BernoulliDistribution[0.9], {5, 5}]; 

mat - DiagonalMatrix[Diagonal[mat]]

%//MatrixForm

You can change the 0.9 to any value (this is your m). The second line sets all the diagonal elements to zero.

share|improve this answer

answered 6 hours ago

bill sbill s

53.6k376153

$endgroup$

mat = RandomVariate[BernoulliDistribution[0.9], {5, 5}]; 

mat - DiagonalMatrix[Diagonal[mat]]

%//MatrixForm

share|improve this answer

answered 6 hours ago

bill sbill s

53.6k376153

answered 6 hours ago

bill sbill s

53.6k376153

answered 6 hours ago

bill sbill s

53.6k376153

3

$begingroup$

You can use weight option in RandomChoice

n = 5;

m = 2;

mat = RandomChoice[{1/m, 1 - 1/m} -> {1, 0}, {n, n}];

(mat - DiagonalMatrix[Diagonal@mat]) // MatrixForm

$left(
begin{array}{ccccc}
0 & 1 & 1 & 0 & 1 \
1 & 0 & 1 & 0 & 1 \
1 & 0 & 0 & 0 & 1 \
0 & 0 & 1 & 0 & 0 \
0 & 0 & 0 & 1 & 0 \
end{array}
right)$

share|improve this answer

answered 4 hours ago

Okkes DulgerciOkkes Dulgerci

5,0891917

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Altho using the BernoulliDistribution[] is best, this is likely to be more easy to read for someone who is not accustomed to discrete probability distributions.
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

3

$begingroup$

You can use weight option in RandomChoice

n = 5;

m = 2;

mat = RandomChoice[{1/m, 1 - 1/m} -> {1, 0}, {n, n}];

(mat - DiagonalMatrix[Diagonal@mat]) // MatrixForm

$left(
begin{array}{ccccc}
0 & 1 & 1 & 0 & 1 \
1 & 0 & 1 & 0 & 1 \
1 & 0 & 0 & 0 & 1 \
0 & 0 & 1 & 0 & 0 \
0 & 0 & 0 & 1 & 0 \
end{array}
right)$

share|improve this answer

answered 4 hours ago

Okkes DulgerciOkkes Dulgerci

5,0891917

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Altho using the BernoulliDistribution[] is best, this is likely to be more easy to read for someone who is not accustomed to discrete probability distributions.
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

$begingroup$

You can use weight option in RandomChoice

n = 5;

m = 2;

mat = RandomChoice[{1/m, 1 - 1/m} -> {1, 0}, {n, n}];

(mat - DiagonalMatrix[Diagonal@mat]) // MatrixForm

$left(
begin{array}{ccccc}
0 & 1 & 1 & 0 & 1 \
1 & 0 & 1 & 0 & 1 \
1 & 0 & 0 & 0 & 1 \
0 & 0 & 1 & 0 & 0 \
0 & 0 & 0 & 1 & 0 \
end{array}
right)$

share|improve this answer

answered 4 hours ago

Okkes DulgerciOkkes Dulgerci

5,0891917

$endgroup$

n = 5;

m = 2;

mat = RandomChoice[{1/m, 1 - 1/m} -> {1, 0}, {n, n}];

(mat - DiagonalMatrix[Diagonal@mat]) // MatrixForm

share|improve this answer

answered 4 hours ago

Okkes DulgerciOkkes Dulgerci

5,0891917

answered 4 hours ago

Okkes DulgerciOkkes Dulgerci

5,0891917

answered 4 hours ago

Okkes DulgerciOkkes Dulgerci

5,0891917

0

$begingroup$

m = 2;

n = 5;

rnd[x_] := If[RandomReal[{0, 1}] < 1/m, 1, 0];

t = Table[rnd[x_]*(1 - KroneckerDelta[i, j]), {i, 1, n}, {j, 1, n}];

t // MatrixForm

share|improve this answer

answered 6 hours ago

Vsevolod A.Vsevolod A.

478211

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Thank you Vsevolod A. How will the command be if I do not restrict what n and m be?
  $endgroup$
  – tiffany
  6 hours ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @tiffany you set m and n in the first two lines...
  $endgroup$
  – Vsevolod A.
  4 hours ago
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Just rnd := If[RandomReal[{0, 1}] < 1/m, 1, 0]; will do, since the function never uses its argument.
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

0

$begingroup$

m = 2;

n = 5;

rnd[x_] := If[RandomReal[{0, 1}] < 1/m, 1, 0];

t = Table[rnd[x_]*(1 - KroneckerDelta[i, j]), {i, 1, n}, {j, 1, n}];

t // MatrixForm

share|improve this answer

answered 6 hours ago

Vsevolod A.Vsevolod A.

478211

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  Thank you Vsevolod A. How will the command be if I do not restrict what n and m be?
  $endgroup$
  – tiffany
  6 hours ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @tiffany you set m and n in the first two lines...
  $endgroup$
  – Vsevolod A.
  4 hours ago
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Just rnd := If[RandomReal[{0, 1}] < 1/m, 1, 0]; will do, since the function never uses its argument.
  $endgroup$
  – J. M. is computer-less♦
  1 hour ago
  

  
  
  
  

add a comment | 

$begingroup$

m = 2;

n = 5;

rnd[x_] := If[RandomReal[{0, 1}] < 1/m, 1, 0];

t = Table[rnd[x_]*(1 - KroneckerDelta[i, j]), {i, 1, n}, {j, 1, n}];

t // MatrixForm

share|improve this answer

answered 6 hours ago

Vsevolod A.Vsevolod A.

478211

$endgroup$

m = 2;

n = 5;

rnd[x_] := If[RandomReal[{0, 1}] < 1/m, 1, 0];

t = Table[rnd[x_]*(1 - KroneckerDelta[i, j]), {i, 1, n}, {j, 1, n}];

t // MatrixForm

share|improve this answer

answered 6 hours ago

Vsevolod A.Vsevolod A.

478211

answered 6 hours ago

Vsevolod A.Vsevolod A.

478211

answered 6 hours ago

Vsevolod A.Vsevolod A.

478211