Código da função

   granulosa<-function(Horário,IndivíduoG1,Grupo1,IndivíduosG2,Grupo2)
             {
Y<-c(Grupo1,Grupo2)
N<-length(Grupo1)
Grupo<-rep(c("Grupo 1","Grupo 2"),each=N)
X<-c(Horário,Horário)
Indivíduos<-c(IndivíduosG1,IndivíduosG2) 
data<-data.frame(Indivíduos,Grupo,X,Y) 

mediaporhora<-with(data, tapply(Y, list(Grupo,X), mean)) 
desviopadraoporhora<-with(data, tapply(Y, list(Grupo,X), sd)) 
medianaporhora<-with(data, tapply(Y, list(Grupo,X), median)) 
minimoporhora<-with(data, tapply(Y, list(Grupo,X), min)) 
maximoporhora<-with(data, tapply(Y, list(Grupo,X), max)) 

lista1<-list(mediaporhora,desviopadraoporhora,medianaporhora,minimoporhora,maximoporhora)
names(lista1)<-c("Média por hora","Desvio padrão por hora","Mediana por hora","Mínimo por hora","Máximo por  hora")

data$Grupo = factor(data$Grupo, levels=c("Grupo 1","Grupo 2")) 
mediadogrupo<-tapply(data$Y,data$Grupo,mean) 
desviooadraoporgrupo<-tapply(data$Y,data$Grupo, sd) 
medianaporgrupo<-tapply(data$Y,data$Grupo, median)
minimoporgrupo<-tapply(data$Y,data$Grupo, min)
maximoporgrupo<-tapply(data$Y,data$Grupo, max)
lista2<-list(mediadogrupo,desviooadraoporgrupo,medianaporgrupo,minimoporgrupo,maximoporgrupo)
names(lista2)<-c("Média por grupo","Desvio padrão por grupo","Mediana por grupo","Mínimo por grupo","Máximo  por grupo")
quantil<-tapply(data$Y,data$Grupo, quantile)
names(quantil)<-c("Quantils Grupo 1", "Quantils Grupo 2") 

aov.out <- aov(X ~ Grupo * Y + Error(Indivíduos/Y), data=data) 
sumario<-summary(aov.out) 

data2<-data.frame(Grupo1,Grupo2,Horário)
g1<-tapply(data2$Grupo1,data2$Horário,mean)
g2<-tapply(data2$Grupo2,data2$Horário,mean)
Y<-c(g1,g2)
X<-Horário 
M<-length(g1)
Grupo<-rep(c("Grupo 1","Grupo 2"),each=M) 
erropadrao1<-(sd(Grupo1))/(sqrt(M)) 
erropadrao2<-(sd(Grupo2))/(sqrt(M)) 
erropadrão <- c(erropadrao1,erropadrao2) 
press<-data.frame(X,Y,erropadrão,Grupo) 

f1 = ggplot(data = press, aes(x = X, y = Y, Grupo = Grupo) ) 
f2 <- f1 + geom_errorbar(aes(ymin = Y - erropadrão, ymax = Y + 
                               erropadrão), width=0.3) + geom_line() + geom_point(aes(shape=Grupo, fill=Grupo), size=5)  
f3 <- f2 +  scale_x_continuous("Tempo") +
  scale_y_continuous("Medida") + 
  scale_shape_manual(values=c(24,21)) + 
  scale_fill_manual(values=c("white","black")) + 
  theme_bw()

lista4<-list(sumario,f3)
names(lista4)<-c("Sumário ANOVA", "Gráfico exploratório") 
listafinal<-list(lista1,lista2,lista3,lista4) 
return(listafinal) 
    }
    

Carla, a sua proposta 1 está um pouco vaga . De que gráficos e medidas descritivas você está falando.? Pelo que pude entender, uma possibilidade seria plotar o tamanho do edema em função do tempo. É isso? Se for esse o caso, uma simples lm() resolveria e sua função não estaria acrescentando nada novo. Já a proposta B está excessivamente simples, pois , pelo que percebi, a sua função estaria simplesmente somando minutos a partir de determinada hora. Por favor, tente desenvolver mais suas propostas para eu poder ajuda-la.

—- Cristiane

Comentários Concordo com os comentários da Cris. DEVE descrever melhor o que a função fará, o formato dos dados de entrada, quais argumentos terá e o tipo e estrutura do objeto e gráficos de saída . Fique com o primeiro, mas tente colocar em um contexto mais geral

Veja algumas dicas de como estruturar a proposta que incluímos no wiki:

http://ecologia.ib.usp.br/bie5782/doku.php?id=bie5782:01_curso_atual:alunos:trabalho_final:start#a_proposta
Aguardamos suas correções! Avise a Cris quando terminar via email acima para ela ou eu darmos uma olhada. — Alexandre Adalardo de Oliveira 2014/04/25 17:02

Carla, deve manter todas as propostas no wiki. Não apague nenhuma delas para termos o histórico documentado. Sua proposta reelaborada é simples e deve ter cuidado para executá-la. Uma função que faça um bom gráfico geral é uma possibilidade, mas deve ser algo que ofereça algum desafio com o risco de fazer uma função muito básica, o que compromete sua avaliação. Na sua descrição não conseguimos reconhecer o desafio que essa função representa. Além disso, uma variavel contínua (y) e uma “temporal” (x) (tb. contínua), não estão no âmbito da ANOVA.

Sugestão: uma função que receba uma variável continua e uma categórica e faça um gráfico com médias e intervalos de confiança para cada nivel da categórica e rode uma anova, retornando uma tabela de ANOVA além do gráfico. Note que para isso ser um desafio vc. teria que fazer tudo usando apenas as funções básicas de gráfico e (points, lines, segments) e fazer a anova tb. usando apenas as funções aritméticas básicas, se usar as funções prontas, a sua função teria apenas duas linhas… o que não é adequado para o trabalho final.

— Alexandre Adalardo de Oliveira 2014/05/08 09:06