R: diferència entre les revisions

895 bytes afegits ,  2 juny 2005
Línia 1.287: Línia 1.287:
  summary(glm(individual$numero~individual$EDAT+individual$SEXE+individual$CONTR,family=poisson))
  summary(glm(individual$numero~individual$EDAT+individual$SEXE+individual$CONTR,family=poisson))


'''5-Agregamos para tener los datos por ambulatorio'''
'''5-Creación de los datos agregados, con un código manifiestamente mejorable'''
 
baixestotals<-aggregate(cbind(individual$UP,individual$numero,individual$dies),by=list(individual$UP),FUN=sum)
subset(individual,UP=="BR017")
subset(baixestotals,Group.1=="BR017")
baixestotals<-baixestotals[,c(1,3,4)]
names(baixestotals)<-c("UP","Baixes","Dies")
edatm<-aggregate(cbind(individual$UP,individual$EDAT),by=list(individual$UP),FUN=mean)
subset(individual,UP=="BR017")
subset(edatm,Group.1=="BR017")
edatm<-edatm[,c(1,3)]
names(edatm)<-c("UP","EdatM")
treballadors<-aggregate(individual,by=list(individual$UP),FUN=length)
subset(individual,UP=="BR017")
subset(treballadors,Group.1=="BR017")
treballadors<-treballadors[,c(1,2)]
names(treballadors)<-c("UP","Treballadors")
agregadapre<-merge(baixestotals,edatm)
agregada<-merge(agregadapre,treballadors)
names(agregada)
agregada
 
'''6-Algún análisis sobre los datos agregados'''


== Sistema de ayuda ==
== Sistema de ayuda ==