Foro de discusión y soporte para usuarios de R

Un foro abierto a todos los temas relacionados con R que quiera.
Nota: Una url estable a este foro es https://knuth.uca.es/R-foro.

test de Friedman

de Angel Pascual Asura -

Se trata de ver s el nivel sociocultural (A, B y C) influye en el interés por el estudio.
La tabla de datos es la siguiente:

Conjuntos

A

B

C

1

17

15

12

2

15

16

14

3

18

15

17

4

17

18

12

5

14

17

13

6

13

18

19

7

16

16

19

8

16

17

16

9

19

17

15

10

17

18

14

La tabla se ha importado desde Excel con formato csv

> Friedman <- read.table("I:/FicherosEstadistica/Friedman.csv", header=TRUE, sep=";", na.strings="NA", dec=".",

+   strip.white=TRUE)

> attach(Friedman)

 

Me da este error al intertar realizar el test de Friedman y no encuentro la razón

> friedman.test(A,B,C)

Error in friedman.test.default(A, B, C) :

  not an unreplicated complete block design


Paralelización R

de Carmen Gonzalez Donquiles -

Buenas, 

Estoy usando un script que contiene una parte que me tarda muchisimo tiempo en ejecutar.

Querría saber si alguien me puede ayudar a paralelizarlo. Copio la parte del Script.

> example.randPathDS <- sapply(1:10000,     

                                    function(run){

                                      podaDS (genotypes, status, 

                                             pathGeneList, geneSnpList, 

                                             NsamplePerms=1000, resamplePaths=TRUE)

                                    }

)


Gracias

problema con summary

de Angel Pascual Asura -

He realizado el siguiente análisis en RStudio: AnovaVariableX <- summary(aov(AnteriorX ~ Grupo, data=Anacova)), que da como resultado

> AnovaVariableX
            Df Sum Sq Mean Sq F value  Pr(>F)   
Grupo        2  36.75  18.375   6.785 0.00533 **
Residuals   21  56.88   2.708             

Quiero extraer la suma de cuadrados residual y no encuentro la forma de hacerlo; cualquier instrucción me da como resultado NULL; por ejemplo AnovaVariableX$`Sum Sq`.

Sin embargo en otras listas sí que lo consigo

HOLA, SOY PRINCIPIANTE EN R

de Zoe Navarro Apodaca -

Mi duda es ... que significamente estla linea de instruccion

plot(feb[,2])
Exactamente el 2 en corchetes a que se refiere?

He cambiado el los y hecho de nuevo la grafica, pero aun no entiendo que es lo que modifica el hecho de cambiar el dos por otros numeros


Muchisimas gracias

Solución sin bucles

de Roberto Caride -
Hola a todos.

Necesito vuestra ayuda.

Tengo un data.frame con más de 100.000 filas y 31 variables.


Precios
x1 x2 x3 x...31
y1 1.3465 1.3464 1.3471 1.3466
y...+100000 1.3443 1.3446 1.3451 1.3442
Necesito hacer que Precios[y1,x1]= Precios [y1,x2 - x1]. Por ejemplo Preciostemp[1,1]= 1.3464-1.3465 = -1.

Pero el problema es que lo tengo que hacer si bucles, ahora mismo lo tengo con bucles y es inviable.

Nota:sólo necesito las 30 primera variables de cada fila, alimento la fila con 31 para poder calcular el número 30.


Un saludo y muchas gracias.



Acceso a lista

de David Tarin Bernad -


Tengo una duda, quiero coger numero por numero en una lista de 100 datos, os paso el script a ver si me podeis ayudar.Quiero que Rho, vaya cogiendo datos uno a uno para que me calcule W, ya que me da mal.

rm(list=ls())

r=matrix(c(0.095, 0.070, 0.090, 0.075), nrow=4)

x=matrix(c(0.038, 0.0085, 0.0089, 0.0066, 0.0085, 0.0331, 0.0156, 0.0039, 0.0089, 0.0156, 0.0334, 0.007, 0.0066, 0.0039, 0.007, 0.0240),nrow=4)

varcov=solve(x)

t<-matrix(c(1,1,1,1),nrow=4)

rho <- seq(0.05, 0.15, 0.001)

lrho=length(rho)

a=t(t)%*%solve(x)%*%t

ba=t(t)%*%solve(x)%*%r

c=t(r)%*%solve(x)%*%r

d=a%*%c-ba^2

g=varcov%*%((t%*%c-r%*%ba))%*%solve(d)

h=varcov%*%((r%*%a-t%*%ba))%*%solve(d)

i=seq(1,lrho)

w=g+h*rho[i]


empezando con R - me da este error

de Ernesto Szini -


Buenos Días Estimados Coforeros 

Estoy empezando a utilizar R, y ya me estoy dando contra las paredes ...

Estoy siguiendo unas prácticas de un libro que utiliza R, y en este contexto, bajé

código ya preparado.


Estoy usando la versión 3.3.2, sobre windows ..  e intento instalar el package "e1071"

inicialmente puse

options(repos=structure(c(CRAN="http://cran.stat.auckland.ac.nz/")))

y luego

install.packages("e1071", dependencies = "TRUE", type = "source")

y me da este error :


> install.packages("e1071", dependencies = "TRUE", type = "source")

Installing package into ‘C:/Users/DATAIQ/Documents/R/win-library/3.3’

(as ‘lib’ is unspecified)

Error in available[p1, dependencies, drop = FALSE] : 

  subíndice fuera de  los límites

No se por donde buscar que es lo que está mal !!

Alguien me podrá orientar ??


desde ya muchas gracias,


Ernesto


eszini at gmail.com








Pregunta

de Miguel Angel Hernandez Martinez -

Hola, estoy comenzando a utilizar R, tengo una base de datos de encuestas aplicadas, las respuestas en su mayoría son Si, No, algunas tienen A,B....D, mi duda es como elaborar la base de datos que se analizara en R,

Es decir, capturo como SI, NO, A,B,....D, o pongo algún código como: 1 para Si, 2 para No, 1 para A, 2 para B, 3 para C y 4 para D, en sus respectivas variables???

Agradeceré mucho alguna orientación, gracias y excelente día!!!



Duda

de David Tarin Bernad -

Tengo una duda, quiero coger numero por numero en una lista de 100 datos, os paso el script a ver si me podeis ayudar.Quiero que Rho, vaya cogiendo datos uno a uno para que me calcule W, ya que me da mal. 

rm(list=ls())

r=matrix(c(0.095, 0.070, 0.090, 0.075), nrow=4)

x=matrix(c(0.038, 0.0085, 0.0089, 0.0066, 0.0085, 0.0331, 0.0156, 0.0039, 0.0089, 0.0156, 0.0334, 0.007, 0.0066, 0.0039, 0.007, 0.0240),nrow=4)

varcov=solve(x)

t<-matrix(c(1,1,1,1),nrow=4)

rho <- seq(0.05, 0.15, 0.001)

lrho=length(rho)

a=t(t)%*%solve(x)%*%t

ba=t(t)%*%solve(x)%*%r

c=t(r)%*%solve(x)%*%r

d=a%*%c-ba^2

g=varcov%*%((t%*%c-r%*%ba))%*%solve(d)

h=varcov%*%((r%*%a-t%*%ba))%*%solve(d)

i=seq(1,lrho)

w=g+h*rho[i]


Examen R-Studio

de Alejandro Solis -

Hola a todos,

Mañana tengo un examen de R y me gustaría saber si hay alguien que vaya a estar disponible durante las horas del examen para poder ayudarme con este.

El examen es de 10.00 a 13.00 (España) y tenemos permiso de usar internet, apuntes y cualquier recurso que creamos necesario y por eso publico esto a ver si alguien pudiese ayudarme ya que sus exámenes suelen ser bastante difíciles.

Si prefieren que se haga por correo envienme un correo a benjanieto@gmail.com

Muchas gracias y un saludo.