Man mano che $ N $ cresce, la tua capacità di usare la matematica diventa assolutamente cruciale. la matematica "inefficiente" ti costerà al PC. Il limite superiore dipende dall'equazione che stai risolvendo. Evitare i calcoli inversi o determinanti della matrice è un grande vantaggio.
Un modo per aiutare ad aumentare il limite è usare i teoremi per scomporre una matrice inversa grande da in inversi di matrice più piccoli. Questo spesso può significare la differenza tra fattibile e non fattibile. Ma questo comporta un duro lavoro e manipolazioni matematiche spesso piuttosto complicate! Ma di solito ne vale la pena. Fai i conti o il tempo!
I metodi bayesiani potrebbero essere in grado di fornire un modo alternativo per ottenere il tuo risultato - potrebbero essere più veloci, il che significa che il tuo "limite superiore" aumenterà (se non altro perché ti dà due modi alternativi per calcolare la stessa risposta: il più piccolo tra due sarà sempre più piccolo di uno di essi!).
Se puoi calcolare un coefficiente di regressione senza invertire una matrice, probabilmente salverai un Un sacco di tempo. Ciò può essere particolarmente utile nel caso bayesiano, perché "all'interno" di un normale integrale di marginalizzazione, la matrice $ X ^ {T} X $ non ha bisogno di essere invertita, si calcola solo una somma di quadrati. Inoltre, la matrice determinante farà parte della costante di normalizzazione. Quindi "in teoria" potresti usare tecniche di campionamento per valutare numericamente l'integrale (anche se ha un'espressione analitica) che sarà eoni più veloce rispetto al tentativo di valutare l '"esplosione combinatoria" di inverse di matrice e determinanti. (sarà ancora una "esplosione combinatoria" di integrazioni numeriche, ma potrebbe essere più veloce).
Questo suggerimento sopra è un po 'una mia "bolla di pensiero". Voglio provarlo davvero, vedere se va bene. Penso che sarebbe (5.000 simulazioni + calcolare exp (somma dei quadrati) + calcolare il beta medio ponderato dovrebbe essere più veloce dell'inversione di matrice per una matrice abbastanza grande.)
Il costo è approssimativo piuttosto che una stima esatta. Non c'è niente che ti impedisca di usare lo stesso insieme di numeri pseudo casuali per valutare numericamente l'integrale, il che ti farà risparmiare molto tempo.
Niente ti impedisce anche di usare una combinazione di entrambe le tecniche. Usa l'esatto quando le matrici sono piccole, usa la simulazione quando sono grandi. Questo perché in questa parte dell'analisi. Sono solo tecniche numeriche diverse: scegli la tecnica più veloce!
Ovviamente questo è solo un po 'di argomenti "ondeggianti", non conosco esattamente i migliori pacchetti software da usare - e peggio ancora, cercando di capire quali algoritmi usano effettivamente.