reenviando...
myFormat() apenas formata a data para mes/ano... mas agora que estou
acordado, nem precisava do lubridate, ne'?
Abaixo a solucao sem frescura:
theSeq - with(dados, sort(format(seq(min(Data), max(Data),
by='month'), '%m/%Y')))
with(dados, table(factor(format(Data, '%m/%Y'),
Tem certeza que o $k$ não está sendo redefinido em algum lugar e que sempre
length(p)=1?
2012/9/14 Felipe Buchbinder felb...@gmail.com
Fábio,
Muito obrigado pela ajuda. Mas o *if *já não dá conta do caso em que
p.vec não possui um elemento maior do que seed? E como p.vec é construido
dentro
ifelse quando length é igual a 1 não é recomendável:
system.time(replicate(1, if(10) {1} else { 0}))
user system elapsed
0.028 0.000 0.030
system.time(replicate(1, ifelse(10, 1, 0)))
user system elapsed
0.224 0.000 0.224
2012/9/14 Fabio Mathias Corrêa
E o problema não está no which...
resp - function() {
k - 5
p.vec = NULL
for(category in 1:k) # k = 5
{
p = runif(1)
p.vec = c(p.vec,p)
}
##Generates a random number
seed-runif(1,0.2,0.8)
##Selects an answer
temp = which ( p.vecseed )
if( length (temp) 0 ) {
Eduardo e Fábio,
Muito obrigado pela ajuda. Se o problema não está no which, alguma idéia de
onde ele possa estar?
A função *resp* já retorna resultados absurdos, então o problema tem que
estar nela. Só não vejo onde
2012/9/14 Eduardo Leoni e.le...@gmail.com
E o problema não está no
*Felipe, *
*
*
*pelo que entendi você esta considerando modelos de teoria de resposta ao
item certo? E dados os parâmetros e o score simulado irá obter o padrão de
respostas certo?*
*
*
Em 14 de setembro de 2012 23:00, Felipe Buchbinder felb...@gmail.comescreveu:
Eduardo e Fábio,
Muito obrigado
certo
2012/9/14 Leandro Marino leandromar...@leandromarino.com.br
*Felipe, *
*
*
*pelo que entendi você esta considerando modelos de teoria de resposta ao
item certo? E dados os parâmetros e o score simulado irá obter o padrão de
respostas certo?*
*
*
Em 14 de setembro de 2012 23:00,
*Tenho algo que fiz e está bem implementado*
*
*
*
PudT - function(Theta,ai,bi,ci,D){
PudT - ci + (1-ci)/(1 + exp(-D*ai*(Theta-bi)))
PudT
}
parit - data.frame(id= ..., a = ..., b=..., c=)
respo - matrix(ncol=nitem,nrow=nalu)
for(k in 1:nitem){
respo[,k] -
Coloque o código completo, com dados.
2012/9/14 Felipe Buchbinder felb...@gmail.com
Eduardo e Fábio,
Muito obrigado pela ajuda. Se o problema não está no which, alguma idéia
de onde ele possa estar?
A função *resp* já retorna resultados absurdos, então o problema tem que
estar nela. Só não
