Gerando variáveis aleatórias que satisfazem restrições

9

Preciso gerar uma lista de variáveis aleatórias sujeitas a restrições que podem ser expressas no formato que é uma matriz se tem entradas. Em todos os casos com os quais estou lidando, $\bf{x}$ $\bf{E}x=b$ $\bf{E}$ $m \times n$ $\bf{x}$ $n$ $n >> m$ , por exemplo, será em torno de 14.000 e será 50. Não tenho certeza de qual método utilizarei para amostragem aleatória, normal ou uniforme, não é claro qual é o melhor para o problema que estou tentando resolver, mas preciso que todas as variáveis sejam amostradas de distribuições com a mesma média e faixa / variação. $n$ $m$

O que eu tenho feito para resolver isso é reduzir para a forma escalonar linhas, definir todas as variáveis correspondentes às colunas à direita do último pivô para valores aleatórios e resolver a igualdade de matriz quadrada restante. $\bf{E}$

Existe um problema, no entanto, para resolver a igualdade da matriz quadrada, subtraio os valores já definidos do lado direito. Infelizmente, as variações também aumentam, então meus últimos 50 valores tendem a variar bastante, o que infelizmente é inaceitável neste problema.

Existe uma maneira melhor de fazer isso? Não consigo pensar em uma maneira de corrigir o método atual que estou usando. Eu uso R.

random-generation random-matrix Mike Flynn
fonte

2

x_{1} = 0

$x_1 = 0$

x_{2} + x_{3} = 1

$x_2+x_3 = 1$

x_{1}

$x_1$

x_{2}

$x_2$

x_{3}

$x_3$

x_{2}

$x_2$

x_{3} = 0

$x_3 = 0$

jbowman

1

$E$

http://www.vliz.be/imisdocs/publications/149403.pdf

e aqui está o pacote:

http://cran.r-project.org/web/packages/limSolve/index.html

Mike Flynn
fonte

11

vliz.be/imisdocs/publications/149403.pdf parece não estar disponível, atualize o link

eee

0

Pode parecer trivial (e não muito eficiente em termos de máquina), mas considere repetir o processo até obter uma resposta adequada? De preferência, apenas modificando um subconjunto menor a cada vez.

Você pode criar uma medida de "distância" para o quão longe você está da sua resposta ideal? Isso pode ajudá-lo a "otimizar"?

Sean
fonte

Eu posso tentar isso. Um problema é que eu ainda precisaria garantir que o formato da distribuição fosse o mesmo. Também restrições de tempo serão proibitivas.

Mike Flynn

Gerando variáveis ​​aleatórias que satisfazem restrições

Respostas:

Gerando variáveis aleatórias que satisfazem restrições