R assegnare più elementi di lista allo stesso oggetto

Question

Attualmente ho un ciclo - beh, in realtà un loop in loop, in un modello di simulazione che rallenta con un numero maggiore di individui. Ne ho vettorializzato la maggior parte e ne ho reso molto più veloce. Ma c'è una parte in cui ho assegnare più elementi di una lista come la stessa cosa, semplificando un grande anello al solo compito che voglio raggiungere:

new.matrices[[length(new.matrices)+1]]<-old.matrix 

Con ogni iterazione del ciclo della linea di cui sopra si chiama, e lo stesso oggetto matrice viene assegnato al prossimo nuovo elemento di una lista.

Sto provando a vettorizzare questo - se possibile, o renderlo più veloce di un ciclo o applicare una dichiarazione.

Finora ho provato cose lungo le linee di:

indices <- seq(from = length(new.matrices) + 1, to = length(new.matrices) + reps) 
new.matrices[indices] <- old.matrix 

Tuttavia questo si traduce nel messaggio:

Warning message: 
In new.effectors[effectorlength] <- matrix : 
    number of items to replace is not a multiple of replacement length 

Si cerca anche di assegnare un valore del old.matrix ad un elemento di new.matrices in questo modo:

[[1]] 
[1] 8687 

[[2]] 
[1] 1 

[[3]] 
[1] 5486 

[[4]] 
[1] 0 

Quando il risultato desiderato è una lista eleme nt = una matrice intera, una copia di old.matrix

C'è un modo per vettorizzare una matrice negli elementi di lista senza eseguire il ciclo? Con loops come è attualmente implementato stiamo parlando di molte migliaia di ripetizioni che rallentano notevolmente le cose, da qui il mio desiderio di vettorializzare questo se possibile.

Answer 1

Probabilmente hai già risolto il problema, in ogni caso, il problema nel codice

new.matrices[indices] <- old.matrix 

è stato causato da cercando di sostituire alcuni oggetti (le NULL elementi nel proprio elenco new.matrices) con qualcosa di diverso, un matrix. Quindi, R importa old.matrix in un vettore e prova ad attaccare ogni singolo valore a un altro elemento di elenco, (questo è il motivo per cui hai ottenuto questo risultato e quando, per esempio, reps è 4 o 8 e old.matrix NON è una matrice 2 x 2, anche tu prendi il warning). Facendo

new.matrices[indices] <- list(old.matrix) 

funzionerà e R replicherà l'elenco singolo elemento list(old.matrix) "reps" tempi automaticamente.