[R-br] Ajuda em uma função
Sérgio Henrique almeida da silva ju
sergio.edfisica em gmail.com
Terça Junho 19 18:33:02 BRT 2012
Oi Benilton
Funcionou direitinho!
Valeu pela ajuda. Eu e meu estatístico local agradecemos
Em 19 de junho de 2012 18:21, Benilton Carvalho
<beniltoncarvalho em gmail.com>escreveu:
> a correcao da minha funcao eh
>
> f = function(dados){
> res = vector('numeric', ncol(dados)/2)
> cols = seq(1, ncol(dados), 2)
> for (i in 1:length(cols))
> res[i] = kappa(dados[, c(cols[i], cols[i]+1)])
> res
> }
>
> deixo para vc (e, possivelmente, seu estatistico local) entender(em) o
> motivo.
>
> se vc quer uma lista que contenha um elemento chamado Kappa que seja
> uma matrix de 1 coluna...
>
> z = list(Kappa=matrix(f(dados)))
>
> ou altere sua funcao para:
>
> f = function(dados){
> res = vector('numeric', ncol(dados)/2)
> cols = seq(1, ncol(dados), 2)
> for (i in 1:length(cols))
> res[i] = kappa(dados[, c(cols[i], cols[i]+1)])
> list(Kappa=matrix(res))
> }
>
> b
>
> 2012/6/19 Sérgio Henrique almeida da silva ju <sergio.edfisica em gmail.com>:
> > Oi Benilton
> >
> > Aumentei meu conjunto de dados:
> >
> > a = c(1,3,5,4,3,2,1,4,5,3,4,5,2,1,1,3,4,5,5,3,2)
> > a1 = c(4,5,5,4,4,3,2,1,4,3,4,5,3,2,3,4,5,5,4,4,1)
> > d = c(3,3,4,4,5,4,3,2,4,5,5,3,4,4,5,4,3,2,3,4,5)
> > d1 = c(2,2,3,3,4,3,2,3,4,5,5,5,4,3,2,3,4,5,4,3,2)
> > e = c(1,3,5,4,3,2,1,4,5,3,4,5,2,1,1,3,4,5,5,3,2)
> > e1 = c(4,5,5,4,4,3,2,1,4,3,4,5,3,2,3,4,5,5,4,4,1)
> > f = c(1,3,5,4,3,2,1,4,5,3,4,5,2,1,1,3,4,5,5,3,2)
> > f1= c(3,3,4,4,5,4,3,2,4,5,5,3,4,4,5,4,3,2,3,4,5)
> >
> > dados = data.frame(a,a1,d,d1,e,e1,f,f1)
> >
> > Porque quando rodo a função dá esse resultado?
> >
> > f = function(dados){
> > res = vector('numeric', ncol(dados)/2)
> > for (i in seq(1, ncol(dados), 2))
> > res[i] = kappa(dados[, c(i, i+1)])
> > res
> > }
> >
> >> f(dados)
> > [1] 6.462970 0.000000 5.618469 0.000000 6.462970 NA 4.358520
> >
> >
> >
> > Outra duvida
> >
> > Como colar esses resultados numa lista?
> >
> > inclui os comandos
> >
> > f = function(x){
> > res = vector('numeric', ncol(x)/2)
> > for (i in seq(1, ncol(x), 2))
> > res[i] = kappa(x[, c(i, i+1)])
> > res
> > z=cbind(res)
> > rownames(z) <- ????? (sequência que vai do número 1 até n resultados)
> > resultado = list(Kappa = z)
> > resultado
> > }
> >
> >
> >
> > --
> > Sérgio Henrique Almeida da Silva Junior
> > Doutorando em Epidemiologia em Saúde Pública
> > Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ
> > http://lattes.cnpq.br/1611345552843383
> > Tel: (21) 94429486/78101651 id: 123*20942
> >
> >
> >
> > _______________________________________________
> > R-br mailing list
> > R-br em listas.c3sl.ufpr.br
> > https://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.
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://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.
>
--
Sérgio Henrique Almeida da Silva Junior
Doutorando em Epidemiologia em Saúde Pública
Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ
http://lattes.cnpq.br/1611345552843383
Tel: (21) 94429486/78101651 id: 123*20942
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120619/186c4c77/attachment.html>
Mais detalhes sobre a lista de discussão R-br