dilandou
27 de Mayo del 2005

hola , a todos tengo el siguiente problema el cual consiste en que tengo que realizar un programa que me detecte cuando los tomates ya no sirve, se ingresan los tamtes por fecha y peso de alli el porgrama desecha lo que ya pasaron el limite de dias al hacerlo tengo varios problemas con el codigo el cual lo pongo en este mensaje para ver si alguien me podria ayudar se lo agradeceria mucho,
tomates <-read.table ( file="tomates.csv " , header=TRUE)
tomates [ 1 : 3 , ]
s t r ( tomates )
tomates [ , "Fecha " ] <-as . Date ( as . character ( tomates [ , "Fecha " ] ) ,
format="%Y_%m_%d")
tomates [ 1 : 3 , ]
s t r ( tomates )
#
# Hagamos ahora un breve an_alisis descriptivo de los
# datos en la data frame tomates
#
attach( tomates )
summary( tomates ) # Media, desv. t__pica, etc.
length (unique ( Planta ) ) # >Cu_antas plantas hay?
length (unique ( Fecha ) ) # >Cu_antas fechas de recogida?
xtabs ( Peso ~ Planta ) # >Cu_anto ha dado cada planta?
xtabs ( Peso ~ Fecha ) # >Cu_anto se ha recogido por fecha?

boxplot ( Peso ~ as . Date ( Fecha ) , # Evoluci_on de los pesos el el tiempo.
main="Evoluc ion de l peso de l o s f r u t o s " ,
xlab="Fecha " , ylab="Peso en gramos ")

# Necesitaremos ahora etiquetas para
# TODAS las fechas en la muestra (no
# una para cada fecha distinta).


Grs <-aggregate (Peso , # Peso total recolectado en cada
by=l i s t ( as . Date ( Fecha ) ) ,
sum)
# fecha. El argumento "by="ha de
# ser una lista.

dimnames(Grs ) [ [ 2 ] ] <-c ( "Fecha " , "Total Gramos ")
print (Grs [ 1 : 5 , ] )

GrsMedio <-aggregate (Peso ,
by=l i s t ( Fecha ) , # Peso medio por fruto en cada
# fecha. El argumento "by="ha de
mean) # ser una lista.

dimnames(GrsMedio ) [ [ 2 ] ] <-c ( "Fecha " , "Media Gramos ")
print (GrsMedio [ 1 : 5 , ] )
# Especi_camos ahora cu_ales de las
# plantas recibieron un determinado
# tratamiento.

Tratadas <-c ( 1 , 2 , 5 , 6 , 9 , 1 0 , 1 3 , 1 4 , 1 7 , 1 8 )
Cont rol e s <-(1:20) [ _ Tratadas ]

# "Status"es un factor dando el
# status de cada observaci_on: "Tratada"

# si corresponde a una planta tratada,
# control si corresponde a un control.

Status <-i f e l s e ( Planta %in % Tratadas , "Tratada " , "Cont rol ")
Status <-as . factor ( Status )

t .t e s t ( Peso [ Status=="Tratada " ] , Peso [ Status=="Cont rol " ] )

boxplot ( Peso ~ Status ,
main="Comparacion de pe sos de l o s f r u t o s " ,
xlab="Clas e de pl ant a s " , ylab="Peso en gramos ")

FechaFactor <-as . factor ( as . Date ( Fecha ) )
coplot ( Peso ~ FechaFactor j Status ,
xlab="Fecha ")
# Evoluci_on del peso con el
# paso del tiempo para los dos grupos.

Fac torPlanta <-as . factor ( Planta ) # Creamos un factor para obtener
# desgloses por planta.

coplot ( Peso ~ FechaFactor j FactorPlanta ,
xlab="Fecha ")
# Evoluci_on del peso de los frutos
# con el paso del tiempo para cada
# planta.


plot ( FactorPlanta , Peso ,main="Pesos de l o s f r u t o s de cada planta " ,
ylab="Gramos " ,
xlab="Plantas ") # Un plot cuyo eje X es un factor se
# transforma autom_aticamente en en
# boxplot. Este permite examinar la
# distribuci_on de los pesos de los
# frutos de cada planta.

l ibrary ( l a t t i c e ) # Ofrece gr_a_cos alternativos.

d e n s i t y p l o t (~ Peso j FactorPlanta ,
layout=c ( 5 , 4 ) , xlab="Peso " , ylab="Densidad est imada " ,
main="Di s t r ibuc i on de l o s pe sos por planta ")
d e n s i t y p l o t (~ Peso j Status ,
layout=c ( 2 , 1 ) , xlab="Peso " , ylab="Densidad est imada " ,
main="Di s t r ibuc i on de l o s pe sos por t ratami ento ")

bar char t ( as . character ( as . Date ( Fecha ) ) ~ Peso j Status ,
layout=c ( 2 , 1 ) , xlab="Peso " , ylab="Fecha " ,
main="Di s t r ibuc i on de l o s pe sos en e l tiempo por t ratami ento ")

# q()