thr3ads.net - R help es - [R-es] Consulta [Jun 2011]

If this information is useful, please help other people find it:
Share via:

Ricardo Borda

2011-Jun-07 02:32 UTC

[R-es] Consulta

Mario:

tengo la siguiente duda:

Supongamos que tengo la siguiente base de datos:

factor1

respuesta

a

1

a

23

a

3

a

4

b

2

b

3

b

43

b

3

a

23

c

2

c

3

c

3

c

4

a

3


Como hago para crear en R una base de datos que solo abarque los valores de
del nivel a  con su respectiva respuesta. Así:

 factor1

respuesta

a

1

a

23

a

3

a

4

a

23

a

3

Agradezco su ayuda



-- 
Atentamente:

Ricardo A. Borda Hernández
Docente de Estadística

	[[alternative HTML version deleted]]

Antonio Maurandi López

2011-Jun-07 07:59 UTC

head link

[R-es] Consulta

imagino que están mal tabulados, no?
imagino que tu bbdd es así:
factor1 respuesta
a 1
a 23
a 3
a 4
b 2
b 3
b 43
b 3
a 23

Yo haría esto:

df.datos<-read.table(''bbdatos-inicial.txt'',header=T,sep="
")
nueva_bbdd_solofactor_a<-df.datos[which(df.datos$factor1=="a"),]
datos.a<-nueva_bbdd_solofactor_a

Si está bien tabulado y tus datos están en un asola columna habría que 
ordenarlos primero o hacerlo de otra forma.

Un saludo.

El 07/06/11 04:32, Ricardo Borda escribió:> Mario:
>
> tengo la siguiente duda:
>
> Supongamos que tengo la siguiente base de datos:
>
> factor1
>
> respuesta
>
> a
>
> 1
>
> a
>
> 23
>
> a
>
> 3
>
> a
>
> 4
>
> b
>
> 2
>
> b
>
> 3
>
> b
>
> 43
>
> b
>
> 3
>
> a
>
> 23
>
> c
>
> 2
>
> c
>
> 3
>
> c
>
> 3
>
> c
>
> 4
>
> a
>
> 3
>
>
> Como hago para crear en R una base de datos que solo abarque los valores de
> del nivel a  con su respectiva respuesta. Así:
>
>   factor1
>
> respuesta
>
> a
>
> 1
>
> a
>
> 23
>
> a
>
> 3
>
> a
>
> 4
>
> a
>
> 23
>
> a
>
> 3
>
> Agradezco su ayuda
>
>
>
> -- Atentamente: Ricardo A. Borda Hernández Docente de Estadística 
> [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>    
-- 
Antonio Maurandi López
Serv. Cálculo Científico. Apoyo Estadístico.
Servicio de Apoyo a la Investigación (SAI)
Vicerrectorado de Investigación.
Universidad de Murcia

Edif. CAID. Campus de Espinardo.
30100 Murcia
@. amaurandi@um.es
T. 868 88 7315 F. 868 88 7302
www.um.es/sai www.um.es/ae


	[[alternative HTML version deleted]]

Carlos Ortega

2011-Jun-07 08:05 UTC

head link

[R-es] Consulta

Hola,

Si efectivamente tienes los datos formateados de esta manera, para sólo
conseguir los valores de "a" puedes utilizar este sencillo bucle:

v.gd<-0
c.gd<-0
for (i in 1:length(v)) {
 if(v[i]=="a") {
   c.gd<-c.gd+1
   v.gd[c.gd]<-v[i]
   v.gd[c.gd+1]<-v[i+1]
   c.gd<-c.gd+1
 } else next
}


Saludos,
Carlos Ortega
www.qualityexcellence.es


2011/6/7 Ricardo Borda <ricardobordah@yahoo.es>
> Mario:
>
> tengo la siguiente duda:
>
> Supongamos que tengo la siguiente base de datos:
>
> factor1
>
> respuesta
>
> a
>
> 1
>
> a
>
> 23
>
> a
>
> 3
>
> a
>
> 4
>
> b
>
> 2
>
> b
>
> 3
>
> b
>
> 43
>
> b
>
> 3
>
> a
>
> 23
>
> c
>
> 2
>
> c
>
> 3
>
> c
>
> 3
>
> c
>
> 4
>
> a
>
> 3
>
>
> Como hago para crear en R una base de datos que solo abarque los valores de
> del nivel a  con su respectiva respuesta. Así:
>
>  factor1
>
> respuesta
>
> a
>
> 1
>
> a
>
> 23
>
> a
>
> 3
>
> a
>
> 4
>
> a
>
> 23
>
> a
>
> 3
>
> Agradezco su ayuda
>
>
>
> --
> Atentamente:
>
> Ricardo A. Borda Hernández
> Docente de Estadística
>
>        [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
	[[alternative HTML version deleted]]

Ricardo Borda

2011-Jun-14 19:41 UTC

head link

[R-es] Consulta

Buenas tardes

Quiero preguntarles lo siguiente:

En R, ¿cómo verifico si una matriz es defina positiva y semidefinida
positiva?
-- 
Atentamente:

Ricardo A. Borda Hernández

	[[alternative HTML version deleted]]

Carlos J. Gil Bellosta

2011-Jun-14 19:43 UTC

head link

[R-es] Consulta

?eigen

El día 14 de junio de 2011 21:41, Ricardo Borda
<ricardobordah en yahoo.es> escribió:> Buenas tardes
>
> Quiero preguntarles lo siguiente:
>
> En R, ¿cómo verifico si una matriz es defina positiva y semidefinida
> positiva?
> --
> Atentamente:
>
> Ricardo A. Borda Hernández
>
>        [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

Valeska Yaitul

2011-Jun-17 21:20 UTC

head link

[R-es] consulta

hola
Me pueden ayudar por favor.
Tengo una tabla de contingencia :

Damage-Bien           Damage-Mal            NoDamage-Bien
NoDamage-Mal
    21.153846              17.307692                53.846154
   7.692308

y necesito guardar de forma individual cada valor.

Gracias



Saludos Cordiales

Valeska .

	[[alternative HTML version deleted]]

Carlos Ortega

2011-Jun-17 22:10 UTC

head link

[R-es] consulta

Hola,

Una forma es transformar tu tabla a data.frame.

Mira este ejemplo:
> table(state.division, state.region)                    state.region
state.division       Northeast South North Central West
  New England                6     0             0    0
  Middle Atlantic            3     0             0    0
  South Atlantic             0     8             0    0
  East South Central         0     4             0    0
  West South Central         0     4             0    0
  East North Central         0     0             5    0
  West North Central         0     0             7    0
  Mountain                   0     0             0    8
  Pacific                    0     0             0    5> t.df<-as.data.frame(table(state.division, state.region))
> t.df       state.division  state.region Freq
1         New England     Northeast    6
2     Middle Atlantic     Northeast    3
3      South Atlantic     Northeast    0
4  East South Central     Northeast    0
5  West South Central     Northeast    0
6  East North Central     Northeast    0
7  West North Central     Northeast    0
8            Mountain     Northeast    0
9             Pacific     Northeast    0
10        New England         South    0
11    Middle Atlantic         South    0
12     South Atlantic         South    8
13 East South Central         South    4
14 West South Central         South    4
15 East North Central         South    0
16 West North Central         South    0
17           Mountain         South    0
18            Pacific         South    0
19        New England North Central    0
20    Middle Atlantic North Central    0
21     South Atlantic North Central    0
22 East South Central North Central    0
23 West South Central North Central    0
24 East North Central North Central    5
25 West North Central North Central    7
26           Mountain North Central    0
27            Pacific North Central    0
28        New England          West    0
29    Middle Atlantic          West    0
30     South Atlantic          West    0
31 East South Central          West    0
32 West South Central          West    0
33 East North Central          West    0
34 West North Central          West    0
35           Mountain          West    8
36            Pacific          West    5


Y sobre el data.frame el acceso a cada dato es directo:
> t.df[1,]  state.division state.region Freq
1    New England    Northeast    6



Saludos,
Carlos Ortega
www.qualityexcellence.es


2011/6/17 Valeska Yaitul <valeskayaitul@gmail.com>
> hola
> Me pueden ayudar por favor.
> Tengo una tabla de contingencia :
>
> Damage-Bien           Damage-Mal            NoDamage-Bien
> NoDamage-Mal
>    21.153846              17.307692                53.846154
>   7.692308
>
> y necesito guardar de forma individual cada valor.
>
> Gracias
>
>
>
> Saludos Cordiales
>
> Valeska .
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
	[[alternative HTML version deleted]]

Ricardo Borda

2011-Jun-19 00:17 UTC

head link

[R-es] Consulta

Buenas tardes


Quiero usar la función *corSymm* de la librería *nlme*  para hacer las
estimaciones de los parámetros  de  un *Linear mixed-effects model. *Esta
metodología es propuesta por  por Pinheiro, J.C., and Bates, D.M. (2000)
"Mixed-Effects Models in S and S-PLUS", Springer.


Quiero obtener una estimación similar a la que genera el siguiente  script:

*library(nlme)*

*spatDat <- data.frame(res=c(2,3,4,2,4), x = c(0,0.25,0.5,0.75,1), y
c(0,0.25,0.5,0.50,0.75), pob=c(12,23,34,45,32), otro3=c(1,1,1,1,1))*

*cs1Exp <- corExp( 1, form = ~ x + y )*

*cs1Exp <- Initialize( cs1Exp, spatDat )*

*pruebaq <- glmmPQL(res~ x + y + offset(log(pob)), data = spatDat, family
poisson, random = ~1 |otro3,  correlation = cs1Exp)*

*pruebaq*



*Mi interés recae sobre el siguiente script:
*

*
*

mat<-corMatrix( cs1Exp ) # Estoy partiendo de una matriz de covarianzas
cualquiera. Por comodidad parto de la matriz definida en el anterior
ejemplo.

superior<-t(mat)[lower.tri(t(mat))]

super <- corSymm( value = superior, form = ~1|otro3)

super<- Initialize(super, data = spatDat)

pruebaq2 <- glmmPQL(res~ x + y + offset(log(pob)), data = spatDat, family
poisson, random = ~1 |otro3,  correlation = super)


*El error que me genera es el siguiente: *

*iteration 1*

*Error en lme.formula(fixed = zz ~ x + y + offset(log(pob)), random = ~1 |  :
*

*  nlminb problem, convergence error code = 1*

*  message = false convergence (8)*

*
*

*Quiero saber cómo corregir este error o algún guía sobre este. *

*
*

*Agradezco su atención *

	[[alternative HTML version deleted]]

Ricardo Borda

2011-Jun-30 16:47 UTC

head link

[R-es] Consulta

Buenas tardes

Quisiera saber cómo puedo sacar los valores p utilizando la opción $. En el
archivo adjunto explico con detalle mi consulta.


Agradezco su ayuda



Atentamente:


Ricardo Borda
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20110630/bd79a0cf/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: Consulta.docx
Type: application/vnd.openxmlformats-officedocument.wordprocessingml.document
Size: 54306 bytes
Desc: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20110630/bd79a0cf/attachment-0001.bin>

Carlos Ortega

2011-Jun-30 23:26 UTC

head link

[R-es] Consulta

Hola

Supongo que lo quieres extraer tras hacer un ajuste lineal (lm).
Aquí tienes la solución:

http://tolstoy.newcastle.edu.au/R/e6/help/09/04/10389.html

O mira este ejemplo:


N <- 1000
u <- rnorm(N)
x1 <- rnorm(N)
x2 <- 1 + x1 + rnorm(N)
y <- 1 + x1 + x2 + udf <- data.frame(y,x1,x2)


fit <- lm(y ~ x1 + x2, data = df)summary(fit)

> val<-summary(fit)
> names(val) [1] "call"          "terms"         "residuals"  
"coefficients"
"aliased"
 [6] "sigma"         "df"            "r.squared"  
"adj.r.squared"
"fstatistic"
[11] "cov.unscaled"> cf.v<-val$fstatistic
> cf.v   value    numdf    dendf
2630.194    2.000  997.000> pf(cf.v[1], cf.v[2], cf.v[3], lower.tail=F)value
    0


En este caso, el p-value del modelo que proporciona summary(fit) es de
2.2e-16. Que puede considerarse equivalente al valor de pf... de "0".

Saludos,
Carlos Ortega
www.qualityexcellence.es


2011/6/30 Ricardo Borda <ricardobordah@yahoo.es>
> Buenas tardes
>
> Quisiera saber cómo puedo sacar los valores p utilizando la opción $. En el
> archivo adjunto explico con detalle mi consulta.
>
>
> Agradezco su ayuda
>
>
>
> Atentamente:
>
>
> Ricardo Borda
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
	[[alternative HTML version deleted]]

Jorge Ivan Velez

2011-Jul-01 00:49 UTC

head link

[R-es] Consulta

?str

Saludos,
Jorge


2011/6/30 Ricardo Borda <>
> Buenas tardes
>
> Quisiera saber cómo puedo sacar los valores p utilizando la opción $. En el
> archivo adjunto explico con detalle mi consulta.
>
>
> Agradezco su ayuda
>
>
>
> Atentamente:
>
>
> Ricardo Borda
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
	[[alternative HTML version deleted]]

R help es - Jun 2011 - Consulta

[R-es] Consulta

[R-es] Consulta

[R-es] Consulta

[R-es] Consulta

[R-es] Consulta

[R-es] consulta

[R-es] consulta

[R-es] Consulta

[R-es] Consulta

[R-es] Consulta

[R-es] Consulta