Indexar um vetor no comando "for"
Na rotina abaixo eu preciso gravar cada vetor results para os 16 valores de j. Preciso colocar um índice porque esses 16 resultados formarão uma nova matriz com 16 colunas e 500 linhas. Como indexar os vetores results no comando "for" e montar a nova matriz? for(j in 1:16){ y <- dados[,j] qqnorm(y);qqline(y) mean(y) results <- c() for(i in 1:500){ results[i] <- mean( sample(y,replace=T))} results qqnorm(results);qqline(results) mean (results) } Obrigado Clodoaldo *"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." * Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/> Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 Métodos Numéricos para Engenharia - Matemática Aplicada Instituto Federal Catarinense - Campus Araquari Rodovia BR 280 - km 27 - Cx. Postal 21 CEP 89245-000 - Araquari/SC Fone: (47) 3803-7240
Caro Fernando (e demais colegas) Gostaria de agradecer sua grande ajuda, mas preciso tirar mais outras dúvidas. Abaixo está o código completo e a matriz dados que é base para os cálculos. A saída que inseri é a coluna resultado [,1] Porque não consigo calcular a média dos valores? As saídas serão usadas como vetores (individualmente) em um outro procedimento para uso do pacote "sem" (Modelagem de equações estruturais). Como tira as aspas, ou não tem nada a ver? Porque aparecem valores 4500 e 3500 no meio dos calculados? (Não tem lógica) setwd("C:/Dados") dados <- read.table("Dados_Boot_txt.txt", h=T) dados
dados X4 X5 X15 X6 X11 X12 X13 X16 X17 X18 X19 X20 X25 X26 X27 X29 1 3 4 3 2 3 4 4 2 3 3 3 5 3 4 3 3 2 2 5 5 5 3 5 5 4 5 4 5 4 4 3 5 3 3 4 5 4 3 5 5 5 4 5 6 6 5 4 4 5 4 4 2 6 4 4 2 6 6 4 4 3 6 6 3 4 4 2 5 4 3 1 3 1 4 5 3 3 1 3 1 2 1 2 1 6 5 5 6 5 5 6 6 6 5 5 6 4 5 5 5 4 7 1 7 4 5 7 6 7 5 4 5 5 3 4 4 7 6 8 5 7 6 5 4 6 6 3 7 5 4 3 5 7 5 6 9 5 6 4 3 5 5 5 4 4 4 3 2 2 5 2 5 10 3 6 2 5 5 6 6 6 4 4 5 5 2 2 6 3 11 2 2 3 2 2 4 4 4 3 4 5 4 5 3 4 3 12 2 3 2 2 3 2 5 3 2 4 5 2 4 5 3 2 13 7 7 4 3 3 6 7 6 6 6 6 6 4 4 6 5 14 5 7 6 6 7 6 6 7 7 6 5 5 3 4 6 6 15 3 6 5 6 5 6 6 5 6 6 5 5 4 4 6 6 16 1 2 1 4 1 5 2 2 3 6 5 4 1 2 1 1 17 4 5 2 4 3 4 3 5 5 4 3 3 2 4 3 4 18 5 6 5 5 3 5 5 6 7 6 6 5 4 3 6 6 19 4 3 5 3 1 5 4 2 3 4 4 5 1 4 2 5 20 3 3 3 4 4 3 4 2 2 3 3 3 2 3 1 2
mean(dados[,1])
mean(dados[,1])[1] 3.5
library(sem) library(boot) #Bootstrap (médias das replicações) obs <- function(x){ results <- c() for(i in 1:500){ results[i] <- mean(sample(dados,replace=T)) saida <- paste(results,i,sep="") } return(saida) } resultado<-apply(dados,2,obs) resultado[,1]
resultado[,1] [1] "3.75500" "4.1500" "4.05500" "3.5500" "3.8500" "4.1500" "3.45500" "4.75500" [9] "4.25500" "3.45500" "3.6500" "3.85500" "3.9500" "4.05500" "3.6500" "3.75500" [17] "4.2500" "4.4500" "3.65500" "4.2500" "4500" "4.15500" "5.05500" "3.75500" [25] "3.8500" "4.45500" "4.25500" "3.65500" "3.65500" "4.05500" "4.05500" "3.3500" [33] "3.4500" "4.15500" "3.15500" "3.8500" "3.7500" "4.3500" "4500" "3.15500" [41] "4.55500" "3.95500" "4.1500" "4.05500" "4.2500" "3.6500" "3.95500" "4.35500" [49] "3.65500" "4.75500" "3.7500" "3.65500" "3.6500" "4.85500" "3.8500" "3.05500" [57] "3.85500" "3.8500" "3.8500" "3.35500" "3.7500" "4.25500" "4.25500" "3.8500" [65] "3.95500" "3.25500" "3.9500" "3.75500" "3.7500" "4.45500" "3.65500" "3.4500" [73] "4500" "4500" "3.55500" "3.9500" "3.75500" "4.3500" "3.3500" "4.25500" [81] "4.1500" "4.2500" "4.2500" "4.2500" "3.95500" "4.15500" "4.35500" "4.1500" [89] "4.1500" "4.05500" "3.9500" "3.65500" "3.4500" "3.65500" "4.15500" "3.35500" [97] "3.55500" "3.75500" "3.8500" "3.8500" "4.25500" "4.15500" "4500" "3.95500" [105] "3.95500" "3.65500" "4.4500" "3.6500" "3.75500" "3.55500" "3.9500" "4.3500" [113] "3.9500" "3.55500" "3.35500" "4.1500" "3.65500" "3.75500" "3.6500" "3.05500" [121] "4500" "3.9500" "4.05500" "3.65500" "4.15500" "4.25500" "4500" "4.05500" [129] "4.7500" "3.1500" "3.2500" "3.85500" "3.5500" "3.5500" "4.4500" "3.45500" [137] "3.7500" "3.65500" "4.1500" "4.05500" "3.35500" "4.1500" "3.9500" "3.1500" [145] "4.2500" "4.15500" "4.45500" "3.95500" "3.25500" "3.9500" "3.9500" "3.75500" [153] "3.65500" "3.6500" "3.5500" "3.65500" "3.9500" "4.05500" "4.05500" "3.85500" [161] "4.05500" "3.45500" "3.05500" "3.45500" "4.3500" "4.6500" "3.8500" "4.35500" [169] "4.05500" "4.1500" "4.15500" "4.8500" "3.6500" "3.05500" "4500" "3.95500" [177] "3.6500" "4.4500" "4.25500" "4.35500" "3.6500" "4.2500" "3.9500" "4.55500" [185] "3.9500" "4.3500" "4.15500" "3.35500" "3.6500" "3.9500" "3.85500" "3.45500" [193] "3.65500" "3.45500" "3.7500" "3.45500" "4.3500" "3.7500" "3.4500" "3.55500" [201] "3.75500" "3.85500" "3.8500" "4.65500" "3.8500" "3.6500" "3.7500" "3.45500" [209] "3.65500" "4.2500" "3.5500" "3.75500" "4.4500" "3.8500" "3.45500" "3.5500" [217] "4.5500" "4.1500" "3.85500" "3.6500" "3.15500" "4.55500" "3.5500" "3.65500" [225] "3.8500" "4500" "4.15500" "3.65500" "3.35500" "3.5500" "3.45500" "3.65500" [233] "3.5500" "4.1500" "3.5500" "4.15500" "4.15500" "4.15500" "4.2500" "4.05500" [241] "4.15500" "3.8500" "3.85500" "4.05500" "4.1500" "3.75500" "3.75500" "3.9500" [249] "3.25500" "4.15500" "3.85500" "4.25500" "4500" "3.9500" "4.15500" "4500" [257] "3.9500" "3.95500" "4500" "4.1500" "3.8500" "4.5500" "4.4500" "3.9500" [265] "3.95500" "4.4500" "3.65500" "3.5500" "4500" "3.8500" "3.4500" "4.2500" [273] "3.05500" "3.6500" "3.6500" "4.2500" "3.4500" "3.55500" "4.25500" "3.95500" [281] "3.4500" "3.9500" "4500" "4.05500" "3.75500" "3.85500" "3.95500" "4.3500" [289] "3.1500" "4.35500" "3.75500" "3.25500" "3.95500" "3.65500" "3.75500" "4.45500" [297] "4.2500" "4.2500" "3.8500" "4500" "3.3500" "4.3500" "3.6500" "4.25500" [305] "3.45500" "3.65500" "5.05500" "4.1500" "3.85500" "4.1500" "3.8500" "3.65500" [313] "3.55500" "4.05500" "3.45500" "3.8500" "4.15500" "3.3500" "3.85500" "4.45500" [321] "4.05500" "3.8500" "3.6500" "4.65500" "3.95500" "3.55500" "3.7500" "3.95500" [329] "3.8500" "3.7500" "3.55500" "3.95500" "3.55500" "3.15500" "4500" "3.35500" [337] "4500" "4.2500" "3.85500" "3.6500" "3.8500" "3500" "4.2500" "3.55500" [345] "3.8500" "4.25500" "3.5500" "3.85500" "3.8500" "3.25500" "3.65500" "4.4500" [353] "4500" "3.6500" "4.4500" "3.3500" "4.45500" "4.25500" "4.05500" "3.45500" [361] "4.3500" "3.85500" "4.15500" "3.65500" "3.55500" "4.3500" "3.9500" "4.4500" [369] "4.15500" "4.15500" "3.85500" "4.25500" "3.7500" "3.8500" "4.2500" "3.35500" [377] "3.2500" "3.8500" "4.2500" "3.95500" "3.85500" "3.25500" "4.35500" "3.1500" [385] "3.7500" "3.35500" "4.1500" "4.15500" "3.85500" "4.55500" "3.65500" "3.75500" [393] "3.8500" "3.25500" "3.2500" "3.5500" "3.75500" "3.45500" "4.65500" "4.15500" [401] "3.25500" "4.1500" "4.3500" "4.2500" "3.4500" "3.9500" "4.6500" "4.05500" [409] "4.35500" "3.55500" "3.8500" "3.8500" "3.55500" "3.8500" "3.7500" "4.05500" [417] "3.65500" "3.55500" "3.5500" "4.15500" "3.75500" "3.6500" "3.85500" "3.3500" [425] "3.85500" "2.95500" "3.85500" "4.45500" "3.85500" "3.45500" "3.65500" "4.1500" [433] "3.6500" "3.9500" "3.8500" "3.15500" "3.65500" "3.4500" "3.5500" "4.25500" [441] "4.2500" "4500" "3.6500" "3.6500" "3.25500" "4.25500" "3.95500" "3.1500" [449] "4.2500" "3.15500" "3.65500" "4500" "3.4500" "4.55500" "4.1500" "4.5500" [457] "3.75500" "3.85500" "4500" "3.85500" "4.55500" "3.6500" "4.1500" "4.3500" [465] "4.05500" "3.7500" "3.95500" "3.55500" "3.1500" "3.95500" "3.85500" "3.3500" [473] "3.95500" "4.25500" "3.95500" "4.1500" "4.05500" "4.25500" "3.75500" "4.3500" [481] "4.3500" "4.3500" "4.7500" "4.45500" "2.9500" "3.5500" "3.6500" "3.55500" [489] "4.4500" "4500" "3.4500" "3.5500" "3.45500" "3.35500" "4.65500" "4.3500" [497] "3.8500" "4.3500" "3.95500" "3.65500"
mean(resultado[,1])
mean(resultado[,1])[1] NAWarning message:In mean.default(resultado[, 1]) : argumento não é numérico nem lógico: retornando NA
Muito obrigado *"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." * Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/> Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 Métodos Numéricos para Engenharia - Matemática Aplicada Instituto Federal Catarinense - Campus Araquari Rodovia BR 280 - km 27 - Cx. Postal 21 CEP 89245-000 - Araquari/SC Fone: (47) 3803-7240 Em 17 de maio de 2017 17:42, Fernando Antonio de souza < nandodesouza@gmail.com> escreveu:
Clodoaldo
O código abaixo , ler seus dados e para cada coluna, aplica a função "obs" a saída será uma tabela contendo o número de colunas do seu dado fornecido e cada coluna conterá 500 linhas correspondente as quinhentas amostragens realizadas pela função "sample".
Assim você pode indentificar cada vetor contendo os 500 dados amostrados através do index de coluna por exemplo : resultado[,1]
Espero ter ajudado
A disposição
dados <- matrix(nrow=10,ncol=16,data=rnorm(160))
obs <- function(x){
results <- c() for(i in 1:500){
results[i] <- mean( sample(x,replace=TRUE)) saida <- paste(results,i,sep="")
}
return(saida) }
resultado<-apply(dados,2,obs)
========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/0?redirect=mailto%3Ae-mail%253Anandodesouza%40gmail.com&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy> Lattes: http://lattes.cnpq.br/6519538815038307 <https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/1?redirect=http%3A%2F%2Flattes.cnpq.br%2F6519538815038307&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy> Blog: https://producaoanimalcomr.wordpress.com/ <https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/2?redirect=https%3A%2F%2Fproducaoanimalcomr.wordpress.com%2F&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy> ==========================================
Sent from Nylas Mail <https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/3?redirect=https%3A%2F%2Fnylas.com%3Fref%3Dn1&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy>, the best free email app for work
On Mai 17 2017, at 4:19 pm, Clodoaldo José Figueredo via R-br < r-br@listas.c3sl.ufpr.br> wrote:
Na rotina abaixo eu preciso gravar cada vetor results para os 16 valores de j. Preciso colocar um índice porque esses 16 resultados formarão uma nova matriz com 16 colunas e 500 linhas. Como indexar os vetores results no comando "for" e montar a nova matriz?
for(j in 1:16){ y <- dados[,j] qqnorm(y);qqline(y) mean(y) results <- c() for(i in 1:500){ results[i] <- mean( sample(y,replace=T))} results qqnorm(results);qqline(results) mean (results) }
Obrigado
Clodoaldo
*"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." * Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy>
Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 Métodos Numéricos para Engenharia - Matemática Aplicada Instituto Federal Catarinense - Campus Araquari Rodovia BR 280 - km 27 - Cx. Postal 21 CEP 89245-000 - Araquari/SC Fone: (47) 3803-7240
_______________________________________________ R-br mailing listR-br@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.
Clodoaldo Estou sem acesso a meu computador aqui é por as não posso testar. Tente rodar a função obs se a variável saida. Acredito que ela está transformando a saída em carácter e por isso a saída. Assim q acessar meu PC testo e lhe confirmo. obs <- function(x){ results <- c() for(i in 1:500){ results[i] <- mean(sample(dados,replace=T)) } return(results) } Em 18/05/2017 9:24 AM, "Clodoaldo José Figueredo" < clodoaldo.figueredo@ifc-araquari.edu.br> escreveu: Caro Fernando (e demais colegas) Gostaria de agradecer sua grande ajuda, mas preciso tirar mais outras dúvidas. Abaixo está o código completo e a matriz dados que é base para os cálculos. A saída que inseri é a coluna resultado [,1] Porque não consigo calcular a média dos valores? As saídas serão usadas como vetores (individualmente) em um outro procedimento para uso do pacote "sem" (Modelagem de equações estruturais). Como tira as aspas, ou não tem nada a ver? Porque aparecem valores 4500 e 3500 no meio dos calculados? (Não tem lógica) setwd("C:/Dados") dados <- read.table("Dados_Boot_txt.txt", h=T) dados
dados X4 X5 X15 X6 X11 X12 X13 X16 X17 X18 X19 X20 X25 X26 X27 X29 1 3 4 3 2 3 4 4 2 3 3 3 5 3 4 3 3 2 2 5 5 5 3 5 5 4 5 4 5 4 4 3 5 3 3 4 5 4 3 5 5 5 4 5 6 6 5 4 4 5 4 4 2 6 4 4 2 6 6 4 4 3 6 6 3 4 4 2 5 4 3 1 3 1 4 5 3 3 1 3 1 2 1 2 1 6 5 5 6 5 5 6 6 6 5 5 6 4 5 5 5 4 7 1 7 4 5 7 6 7 5 4 5 5 3 4 4 7 6 8 5 7 6 5 4 6 6 3 7 5 4 3 5 7 5 6 9 5 6 4 3 5 5 5 4 4 4 3 2 2 5 2 5 10 3 6 2 5 5 6 6 6 4 4 5 5 2 2 6 3 11 2 2 3 2 2 4 4 4 3 4 5 4 5 3 4 3 12 2 3 2 2 3 2 5 3 2 4 5 2 4 5 3 2 13 7 7 4 3 3 6 7 6 6 6 6 6 4 4 6 5 14 5 7 6 6 7 6 6 7 7 6 5 5 3 4 6 6 15 3 6 5 6 5 6 6 5 6 6 5 5 4 4 6 6 16 1 2 1 4 1 5 2 2 3 6 5 4 1 2 1 1 17 4 5 2 4 3 4 3 5 5 4 3 3 2 4 3 4 18 5 6 5 5 3 5 5 6 7 6 6 5 4 3 6 6 19 4 3 5 3 1 5 4 2 3 4 4 5 1 4 2 5 20 3 3 3 4 4 3 4 2 2 3 3 3 2 3 1 2
mean(dados[,1])
mean(dados[,1])[1] 3.5
library(sem) library(boot) #Bootstrap (médias das replicações) obs <- function(x){ results <- c() for(i in 1:500){ results[i] <- mean(sample(dados,replace=T)) saida <- paste(results,i,sep="") } return(saida) } resultado<-apply(dados,2,obs) resultado[,1]
resultado[,1] [1] "3.75500" "4.1500" "4.05500" "3.5500" "3.8500" "4.1500" "3.45500" "4.75500" [9] "4.25500" "3.45500" "3.6500" "3.85500" "3.9500" "4.05500" "3.6500" "3.75500" [17] "4.2500" "4.4500" "3.65500" "4.2500" "4500" "4.15500" "5.05500" "3.75500" [25] "3.8500" "4.45500" "4.25500" "3.65500" "3.65500" "4.05500" "4.05500" "3.3500" [33] "3.4500" "4.15500" "3.15500" "3.8500" "3.7500" "4.3500" "4500" "3.15500" [41] "4.55500" "3.95500" "4.1500" "4.05500" "4.2500" "3.6500" "3.95500" "4.35500" [49] "3.65500" "4.75500" "3.7500" "3.65500" "3.6500" "4.85500" "3.8500" "3.05500" [57] "3.85500" "3.8500" "3.8500" "3.35500" "3.7500" "4.25500" "4.25500" "3.8500" [65] "3.95500" "3.25500" "3.9500" "3.75500" "3.7500" "4.45500" "3.65500" "3.4500" [73] "4500" "4500" "3.55500" "3.9500" "3.75500" "4.3500" "3.3500" "4.25500" [81] "4.1500" "4.2500" "4.2500" "4.2500" "3.95500" "4.15500" "4.35500" "4.1500" [89] "4.1500" "4.05500" "3.9500" "3.65500" "3.4500" "3.65500" "4.15500" "3.35500" [97] "3.55500" "3.75500" "3.8500" "3.8500" "4.25500" "4.15500" "4500" "3.95500" [105] "3.95500" "3.65500" "4.4500" "3.6500" "3.75500" "3.55500" "3.9500" "4.3500" [113] "3.9500" "3.55500" "3.35500" "4.1500" "3.65500" "3.75500" "3.6500" "3.05500" [121] "4500" "3.9500" "4.05500" "3.65500" "4.15500" "4.25500" "4500" "4.05500" [129] "4.7500" "3.1500" "3.2500" "3.85500" "3.5500" "3.5500" "4.4500" "3.45500" [137] "3.7500" "3.65500" "4.1500" "4.05500" "3.35500" "4.1500" "3.9500" "3.1500" [145] "4.2500" "4.15500" "4.45500" "3.95500" "3.25500" "3.9500" "3.9500" "3.75500" [153] "3.65500" "3.6500" "3.5500" "3.65500" "3.9500" "4.05500" "4.05500" "3.85500" [161] "4.05500" "3.45500" "3.05500" "3.45500" "4.3500" "4.6500" "3.8500" "4.35500" [169] "4.05500" "4.1500" "4.15500" "4.8500" "3.6500" "3.05500" "4500" "3.95500" [177] "3.6500" "4.4500" "4.25500" "4.35500" "3.6500" "4.2500" "3.9500" "4.55500" [185] "3.9500" "4.3500" "4.15500" "3.35500" "3.6500" "3.9500" "3.85500" "3.45500" [193] "3.65500" "3.45500" "3.7500" "3.45500" "4.3500" "3.7500" "3.4500" "3.55500" [201] "3.75500" "3.85500" "3.8500" "4.65500" "3.8500" "3.6500" "3.7500" "3.45500" [209] "3.65500" "4.2500" "3.5500" "3.75500" "4.4500" "3.8500" "3.45500" "3.5500" [217] "4.5500" "4.1500" "3.85500" "3.6500" "3.15500" "4.55500" "3.5500" "3.65500" [225] "3.8500" "4500" "4.15500" "3.65500" "3.35500" "3.5500" "3.45500" "3.65500" [233] "3.5500" "4.1500" "3.5500" "4.15500" "4.15500" "4.15500" "4.2500" "4.05500" [241] "4.15500" "3.8500" "3.85500" "4.05500" "4.1500" "3.75500" "3.75500" "3.9500" [249] "3.25500" "4.15500" "3.85500" "4.25500" "4500" "3.9500" "4.15500" "4500" [257] "3.9500" "3.95500" "4500" "4.1500" "3.8500" "4.5500" "4.4500" "3.9500" [265] "3.95500" "4.4500" "3.65500" "3.5500" "4500" "3.8500" "3.4500" "4.2500" [273] "3.05500" "3.6500" "3.6500" "4.2500" "3.4500" "3.55500" "4.25500" "3.95500" [281] "3.4500" "3.9500" "4500" "4.05500" "3.75500" "3.85500" "3.95500" "4.3500" [289] "3.1500" "4.35500" "3.75500" "3.25500" "3.95500" "3.65500" "3.75500" "4.45500" [297] "4.2500" "4.2500" "3.8500" "4500" "3.3500" "4.3500" "3.6500" "4.25500" [305] "3.45500" "3.65500" "5.05500" "4.1500" "3.85500" "4.1500" "3.8500" "3.65500" [313] "3.55500" "4.05500" "3.45500" "3.8500" "4.15500" "3.3500" "3.85500" "4.45500" [321] "4.05500" "3.8500" "3.6500" "4.65500" "3.95500" "3.55500" "3.7500" "3.95500" [329] "3.8500" "3.7500" "3.55500" "3.95500" "3.55500" "3.15500" "4500" "3.35500" [337] "4500" "4.2500" "3.85500" "3.6500" "3.8500" "3500" "4.2500" "3.55500" [345] "3.8500" "4.25500" "3.5500" "3.85500" "3.8500" "3.25500" "3.65500" "4.4500" [353] "4500" "3.6500" "4.4500" "3.3500" "4.45500" "4.25500" "4.05500" "3.45500" [361] "4.3500" "3.85500" "4.15500" "3.65500" "3.55500" "4.3500" "3.9500" "4.4500" [369] "4.15500" "4.15500" "3.85500" "4.25500" "3.7500" "3.8500" "4.2500" "3.35500" [377] "3.2500" "3.8500" "4.2500" "3.95500" "3.85500" "3.25500" "4.35500" "3.1500" [385] "3.7500" "3.35500" "4.1500" "4.15500" "3.85500" "4.55500" "3.65500" "3.75500" [393] "3.8500" "3.25500" "3.2500" "3.5500" "3.75500" "3.45500" "4.65500" "4.15500" [401] "3.25500" "4.1500" "4.3500" "4.2500" "3.4500" "3.9500" "4.6500" "4.05500" [409] "4.35500" "3.55500" "3.8500" "3.8500" "3.55500" "3.8500" "3.7500" "4.05500" [417] "3.65500" "3.55500" "3.5500" "4.15500" "3.75500" "3.6500" "3.85500" "3.3500" [425] "3.85500" "2.95500" "3.85500" "4.45500" "3.85500" "3.45500" "3.65500" "4.1500" [433] "3.6500" "3.9500" "3.8500" "3.15500" "3.65500" "3.4500" "3.5500" "4.25500" [441] "4.2500" "4500" "3.6500" "3.6500" "3.25500" "4.25500" "3.95500" "3.1500" [449] "4.2500" "3.15500" "3.65500" "4500" "3.4500" "4.55500" "4.1500" "4.5500" [457] "3.75500" "3.85500" "4500" "3.85500" "4.55500" "3.6500" "4.1500" "4.3500" [465] "4.05500" "3.7500" "3.95500" "3.55500" "3.1500" "3.95500" "3.85500" "3.3500" [473] "3.95500" "4.25500" "3.95500" "4.1500" "4.05500" "4.25500" "3.75500" "4.3500" [481] "4.3500" "4.3500" "4.7500" "4.45500" "2.9500" "3.5500" "3.6500" "3.55500" [489] "4.4500" "4500" "3.4500" "3.5500" "3.45500" "3.35500" "4.65500" "4.3500" [497] "3.8500" "4.3500" "3.95500" "3.65500"
mean(resultado[,1])
mean(resultado[,1])[1] NAWarning message:In mean.default(resultado[, 1]) : argumento não é numérico nem lógico: retornando NA
Muito obrigado *"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." * Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/> Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 Métodos Numéricos para Engenharia - Matemática Aplicada Instituto Federal Catarinense - Campus Araquari Rodovia BR 280 - km 27 - Cx. Postal 21 CEP 89245-000 - Araquari/SC Fone: (47) 3803-7240 Em 17 de maio de 2017 17:42, Fernando Antonio de souza < nandodesouza@gmail.com> escreveu:
Clodoaldo
O código abaixo , ler seus dados e para cada coluna, aplica a função "obs" a saída será uma tabela contendo o número de colunas do seu dado fornecido e cada coluna conterá 500 linhas correspondente as quinhentas amostragens realizadas pela função "sample".
Assim você pode indentificar cada vetor contendo os 500 dados amostrados através do index de coluna por exemplo : resultado[,1]
Espero ter ajudado
A disposição
dados <- matrix(nrow=10,ncol=16,data=rnorm(160))
obs <- function(x){
results <- c() for(i in 1:500){
results[i] <- mean( sample(x,replace=TRUE)) saida <- paste(results,i,sep="")
}
return(saida) }
resultado<-apply(dados,2,obs)
========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/0?redirect=mailto%3Ae-mail%253Anandodesouza%40gmail.com&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy> Lattes: http://lattes.cnpq.br/6519538815038307 <https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/1?redirect=http%3A%2F%2Flattes.cnpq.br%2F6519538815038307&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy> Blog: https://producaoanimalcomr.wordpress.com/ <https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/2?redirect=https%3A%2F%2Fproducaoanimalcomr.wordpress.com%2F&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy> ==========================================
Sent from Nylas Mail <https://n1.nylas.com/link/982181c7e27b6ef6fdf684e5db1087a882a4023a059b93166c5b32c1a9842d8e/3?redirect=https%3A%2F%2Fnylas.com%3Fref%3Dn1&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy>, the best free email app for work
On Mai 17 2017, at 4:19 pm, Clodoaldo José Figueredo via R-br < r-br@listas.c3sl.ufpr.br> wrote:
Na rotina abaixo eu preciso gravar cada vetor results para os 16 valores de j. Preciso colocar um índice porque esses 16 resultados formarão uma nova matriz com 16 colunas e 500 linhas. Como indexar os vetores results no comando "for" e montar a nova matriz?
for(j in 1:16){ y <- dados[,j] qqnorm(y);qqline(y) mean(y) results <- c() for(i in 1:500){ results[i] <- mean( sample(y,replace=T))} results qqnorm(results);qqline(results) mean (results) }
Obrigado
Clodoaldo
*"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." * Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/&r=Y2xvZG9hbGRvLmZpZ3VlcmVkb0BpZmMtYXJhcXVhcmkuZWR1LmJy>
Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 Métodos Numéricos para Engenharia - Matemática Aplicada Instituto Federal Catarinense - Campus Araquari Rodovia BR 280 - km 27 - Cx. Postal 21 CEP 89245-000 - Araquari/SC Fone: (47) 3803-7240
_______________________________________________ R-br mailing listR-br@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.
participantes (2)
-
Clodoaldo José Figueredo -
Fernando Antonio de souza