Hola colegas: Siguiendo la idea que me hizo el favor de darme Carlos Ortega busqué en "tseries" algo que me ayudara a determinar la probabilidad de que un evento "sep" se diera en "tpo" y envié una serie como la que adjunto. La columna "A" tiene el orden en que se llevaron a cabo los eventos. Las columnas "B" y "C" tienen los datos de la sucesión original de eventos "sep" y "tpo", donde "sep" varía entre 1 y 20, mientras que "tpo" es tipo factor que codificado es 1, 2, 3 y 7. Las columnas "E" y "G" siguen teniendo los datos originales. La "F" acumula mientras no cambie "tpo". Por último en las columnas "I" y "J" se han puesto los datos que estaban en rojo con objeto de llegar a obtener mínimo, primer cuartil, mediana, media, desviación típica, tercer cuartil y máximo. Hay alguna(s) funciones en R que puedan hacer lo anterior y alguna(s) con las que se pueda obtener la probabilidad de "tpo". Gracias anticipadas. MANOLO MÁRQUEZ P. ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130116/cd8ca1fb/attachment.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: SERIE 5.xls Type: application/vnd.ms-excel Size: 16384 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130116/cd8ca1fb/attachment.xls>
Manuel No se si entiendo correctamente su pregunta, pero podría hacer (usando irónicamente pero para ser claro y no como falta hacia usted). resultado.final <- data.frame(seno, coseno, suma, original, edad_abuela, probabilidad_de_c, distancia_del_sol) Donde todo lo que coloca en el data frame son resultados calculados o definidos por separado. Disculpe la redacción algo irrespetuosa, pero es para colocar en térmicos entendibles lo que posiblemente no entendí de su pregunta. Javier El mié, 16-01-2013 a las 20:50 -0600, Manuel Máquez escribió:> Hola colegas: > Siguiendo la idea que me hizo el favor de darme Carlos Ortega busqué > en "tseries" algo que me ayudara a determinar la probabilidad de que > un evento "sep" se diera en "tpo" y envié una serie como la que > adjunto. > La columna "A" tiene el orden en que se llevaron a cabo los eventos. > Las columnas "B" y "C" tienen los datos de la sucesión original de > eventos "sep" y "tpo", donde "sep" varía entre 1 y 20, mientras que > "tpo" es tipo factor que codificado es 1, 2, 3 y 7. > Las columnas "E" y "G" siguen teniendo los datos originales. La "F" > acumula mientras no cambie "tpo". > Por último en las columnas "I" y "J" se han puesto los datos que > estaban en rojo con objeto de llegar a obtener mínimo, primer cuartil, > mediana, media, desviación típica, tercer cuartil y máximo. > Hay alguna(s) funciones en R que puedan hacer lo anterior y alguna(s) > con las que se pueda obtener la probabilidad de "tpo". > Gracias anticipadas. > MANOLO MÁRQUEZ P. > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es
Hola Manuel,
Antes de proyectar nada, he mirado simplemente como es tu conjunto para
cada nivel de "tpo" (mira la imagen adjunta) y como ves el
comportamiento
de cada grupo (tpo) es irregular, es muy difícil ver un patrón.
Entre las dudas a las que entiendo que se refiere Javier está una y es que
¿quieres proyectar "sep" teniendo en cuenta el valor de
"tpo"? Si es así
entonces realmente estás queriendo proyectar cuatro series diferentes (una
por cada "tpo") y si es así, estas series son muy, muy irregulares. Si
quieres hacer esto es más complejo y los resultados por la poca cantidad de
datos en algunas series (la 3 y la 7 por ejemplo) ofrecerán unos resultados
de proyección con gran error.
[image: Imágenes integradas 1]
¿O lo que quieres es proyectar "sep" cómo si fuese una única serie sin
tener en cuenta el valor de "tpo"?. En este caso parece que pueda ser
más
fácil... (es la imagen en la que hay un sólo gráfico, también adjunta).
[image: Imágenes integradas 2]
El código para generar lo anterior es este:
#########################
Lines <- "id sep tpo
1 2 1
2 5 1
3 1 3
4 8 2
5 11 2
6 13 1
7 10 1
8 2 3
9 4 7
10 1 1
11 3 1
12 3 1
13 20 2
14 1 1
15 1 1
16 4 2
17 8 1
18 1 1
19 5 3
20 3 2
21 11 1
22 7 2
23 5 7
24 17 1
25 5 1
26 3 3
27 2 1
28 16 3
29 4 2
30 13 1
31 1 1
32 4 1
33 4 1
34 4 1
35 15 7
36 20 1
37 11 2
38 8 2
39 4 1
40 2 2
41 3 1
42 5 1
43 1 1
44 1 2
45 1 1
"
t.con<-textConnection(Lines)
dat.in<-read.table(t.con, header=T, as.is=T)
close(t.con)
library(lattice)
barchart(
sep ~ id | as.factor(tpo)
,horizontal=F
,data=dat.in
,xlab="Tipo"
,ylab="Sep"
,layout=c(1,4)
,as.table=T
)
barchart(
sep ~ id
,horizontal=F
,data=dat.in
,xlab="Tipo"
,ylab="Sep"
)
#########################
Saludos,
Carlos Ortega
www.qualityexcellence.es
El 17 de enero de 2013 07:45, Marcuzzi, Javier <
javier.ruben.marcuzzi en gmail.com> escribió:
> Manuel
>
> No se si entiendo correctamente su pregunta, pero podría hacer (usando
> irónicamente pero para ser claro y no como falta hacia usted).
>
> resultado.final <- data.frame(seno, coseno, suma, original, edad_abuela,
> probabilidad_de_c, distancia_del_sol)
>
> Donde todo lo que coloca en el data frame son resultados calculados o
> definidos por separado.
>
> Disculpe la redacción algo irrespetuosa, pero es para colocar en
> térmicos entendibles lo que posiblemente no entendí de su pregunta.
>
> Javier
>
> El mié, 16-01-2013 a las 20:50 -0600, Manuel Máquez escribió:
> > Hola colegas:
> > Siguiendo la idea que me hizo el favor de darme Carlos Ortega busqué
> > en "tseries" algo que me ayudara a determinar la
probabilidad de que
> > un evento "sep" se diera en "tpo" y envié una
serie como la que
> > adjunto.
> > La columna "A" tiene el orden en que se llevaron a cabo los
eventos.
> > Las columnas "B" y "C" tienen los datos de la
sucesión original de
> > eventos "sep" y "tpo", donde "sep" varía
entre 1 y 20, mientras que
> > "tpo" es tipo factor que codificado es 1, 2, 3 y 7.
> > Las columnas "E" y "G" siguen teniendo los datos
originales. La "F"
> > acumula mientras no cambie "tpo".
> > Por último en las columnas "I" y "J" se han puesto
los datos que
> > estaban en rojo con objeto de llegar a obtener mínimo, primer cuartil,
> > mediana, media, desviación típica, tercer cuartil y máximo.
> > Hay alguna(s) funciones en R que puedan hacer lo anterior y alguna(s)
> > con las que se pueda obtener la probabilidad de "tpo".
> > Gracias anticipadas.
> > MANOLO MÁRQUEZ P.
> > _______________________________________________
> > R-help-es mailing list
> > R-help-es en r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20130117/f4678428/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: no disponible
Type: image/png
Size: 14372 bytes
Desc: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20130117/f4678428/attachment-0002.png>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: no disponible
Type: image/png
Size: 24206 bytes
Desc: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20130117/f4678428/attachment-0003.png>
Estimado Manuel Máquez Mirando la respuesta de Carlos, hay un punto que es la cantidad de registros o ¿274 pruebas son el resumen de muchos datos?... R informará algo pero ... Porque si son solo 274 datos cuándo usted dice pruebas, estamos en problemas. Por otro lado observando su aclaración, posiblemente hay algunas alternativas, pero no se exactamente cuál puede ser. Le propongo escribir en R lo siguiente ?by y ver dentro de la ayuda. Javier Marcuzzi El jue, 17-01-2013 a las 12:31 +0100, Carlos Ortega escribió:> ######################### > > > > Lines <- "id sep tpo 1 2 1 2 5 1 3 1 3 4 8 2 5 11 2 6 13 1 7 10 1 8 > 2 3 9 4 7 10 1 1 11 3 1 12 3 1 13 20 2 14 1 1 15 1 1 16 4 2 17 8 1 18 > 1 1 19 5 3 20 3 2 21 11 1 22 7 2 23 5 7 24 17 1 25 5 1 26 3 3 27 2 1 > 28 16 3 29 4 2 30 13 1 31 1 1 32 4 1 33 4 1 34 4 1 35 15 7 36 20 1 37 > 11 2 38 8 2 39 4 1 40 2 2 41 3 1 42 5 1 43 1 1 44 1 2 45 1 1 " > > > > t.con<-textConnection(Lines) dat.in<-read.table(t.con, header=T, > as.is=T) close(t.con) > > library(lattice) > > barchart( sep ~ id | as.factor(tpo) ,horizontal=F > ,data=dat.in ,xlab="Tipo" ,ylab="Sep" > ,layout=c(1,4) ,as.table=T ) > > > > barchart( sep ~ id ,horizontal=F ,data=dat.in ,xlab="Tipo" > ,ylab="Sep" ) > > ####
Hola, Realmente son sólo 47 medidas. Otros 47 son de codificación (los grupos) y los otros 47 son el número de orden de las medidas... Saludos, Carlos Ortega www.qualityexcellence.es El 17 de enero de 2013 15:42, Marcuzzi, Javier < javier.ruben.marcuzzi@gmail.com> escribió:> Estimado Manuel Máquez > > Mirando la respuesta de Carlos, hay un punto que es la cantidad de > registros o ¿274 pruebas son el resumen de muchos datos?... R informará > algo pero ... > > > Porque si son solo 274 datos cuándo usted dice pruebas, estamos en > problemas. > > Por otro lado observando su aclaración, posiblemente hay algunas > alternativas, pero no se exactamente cuál puede ser. > > Le propongo escribir en R lo siguiente ?by y ver dentro de la ayuda. > > Javier Marcuzzi > > El jue, 17-01-2013 a las 12:31 +0100, Carlos Ortega escribió: > > ######################### > > > > > > > > Lines <- "id sep tpo 1 2 1 2 5 1 3 1 3 4 8 2 5 11 2 6 13 1 7 10 1 8 > > 2 3 9 4 7 10 1 1 11 3 1 12 3 1 13 20 2 14 1 1 15 1 1 16 4 2 17 8 1 18 > > 1 1 19 5 3 20 3 2 21 11 1 22 7 2 23 5 7 24 17 1 25 5 1 26 3 3 27 2 1 > > 28 16 3 29 4 2 30 13 1 31 1 1 32 4 1 33 4 1 34 4 1 35 15 7 36 20 1 37 > > 11 2 38 8 2 39 4 1 40 2 2 41 3 1 42 5 1 43 1 1 44 1 2 45 1 1 " > > > > > > > > t.con<-textConnection(Lines) dat.in<-read.table(t.con, header=T, > > as.is=T) close(t.con) > > > > library(lattice) > > > > barchart( sep ~ id | as.factor(tpo) ,horizontal=F > > ,data=dat.in ,xlab="Tipo" ,ylab="Sep" > > ,layout=c(1,4) ,as.table=T ) > > > > > > > > barchart( sep ~ id ,horizontal=F ,data=dat.in ,xlab="Tipo" > > ,ylab="Sep" ) > > > > #### > >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]