Escribir una función que devuelva el valor medio de dos números

Buenas tardes a todos,

estoy teniendo problemas a la hora de exportar un conjunto de datos he probado varias cosas y no sé cómo solucionarlo.

1 - Función write.table

write.table(resultados, "D:/wd.lopc/d1.8b/resultados.txt", sep="\t",
col.names=TRUE, row.names=TRUE, quote=TRUE, na="NA")

Me sale el mensaje "Conjunto de datos activos exportados al fichero D:/wd.lopc/d1.8b/resultados.txt" y sin embargo el archivo resultados.txt no aparece en el directorio.

2 - También he intentado otra solución que encontré en un post de este foro,

sink("resultados.txt")

print(resultados)

En este caso sí que crea el archivo resultados.txt pero al abrirlo lo único que pone es NULL

3. Por último he probado write.csv

write.csv (resultados, file="resultados.csv")

Y en este caso me sale el error: ERROR: objeto 'xj' no encontrado

¿Alguna idea sobre qué puede estar ocurriendo?

Gracias de antemano,

Un saludo

Nerea

Hola a todos, estoy intentando salvar unos gráficos en formato eps y que previamente veo en pantalla perfectamente.

1: In var(as.vector(x), na.rm = na.rm) : NAs introducidos por coerción

Hola, existe algún comando en R para calcular la matriz de información de Fisher? Gracias.

Hola, estoy utilizando la función vglm para calcular el estadísitco de Tarone y me da el siguiente error:

Error en dotFortran(name = "veigen", as.integer(M), as.integer(n), as.double(x), :
NA/NaN/Inf en llamada a una función externa (arg 3)

Alguien tiene idea de a qué se debe este error o como puedo solucionarlo?

El código es el siguiente:

Zvalue=matrix(nrow=100,ncol=M)
pvalue=matrix(nrow=100,ncol=M)
prob=matrix(nrow=100,ncol=M)
rho=matrix(nrow=100,ncol=M)
low=matrix(nrow=100,ncol=M)

tarone=vector(length=M)
sebench=vector(length=M)

for (j in 1:M) {

w=as.numeric(x[j,]<=g)

p=mean(w)

for (l in j0:100) {

k=l+1 #a lo largo de los bloques dependientes

A=vector(length=k) #sj
n11=s%/%k
for (i in 1:k){A[i]=sum(w[((i-1)*n11+1):(i*n11)])}

A[k]=sum(w[((k-1)*n11+1):s])
B=c(rep(n11,k-1),length(w[((k-1)*n11+1):s])) #nj

fit=vglm(cbind(A,B-A)~1,betabinomial)
pmodel=exp(coef(fit)[1])/(1+exp(coef(fit)[1]))

#Tarone statistic:
#p=sum(A)/sum(B)
S=sum((A-p*B)^2)/(p*(1-p))
Zvalue[l,j]=(S-sum(B))/sqrt(2*sum(B*(B-1)))
pvalue[l,j]=1-pnorm(Zvalue[l,j])

prob[l,j]=pmodel
rho[l,j]=exp(coef(fit)[2])/(1+exp(coef(fit)[2]))
low[l,j]=coef(fit)[1]-sqrt(vcov(fit))[1,1]*qnorm(.95)

}

tarone[j]=pvalue[ki-1,j]
sebench[j]=(log(prob[ki-1,j]/(1-prob[ki-1,j]))-low[ki-1,j])/qnorm(.95)

}

Gracias.

Hola, me gustaría saber si alguien conoce la manera de aleatorizar las posiciones de un vector numérico, de tal forma que a partir de un vector que podría ser

x = 1 2 3 4 5 6 7 8 9

Pueda generar n vectores

n1 = 5 4 7 8 9 1 2 3 6
n2 = 4 9 2 1 6 3 5 7 8
...

Gracias!

Rafa

Hola a todos, en mi pc tengo la distribución 2.12.1 de R y he instalado manualmente el paquete tm para la minería de texto, necesito hacer una aplicación para el manejo con la minería de texto usando Symfony 2, cuando se ejecuta con el comando exec() de php algunas instrucciones de R se ejecutan pero cuando ejecuto la plantilla desde php con la libreria tm no me responde, me devuelve una variable en blanco, sin embargo cuando ejecuto el mismo código en la consola pues me da el resultado esperado, quisiera si alguien pudiera me ayudara con este tema

Saludos y Gracias de antemano

Hola,

me gustaria saber si habeis probado la integracion R en CONDOR (con la libreria MPI). Aunque las pruebas hayan salido mal, tambien me vienen bien si teneis conclusiones al respecto.

Gracias
An3

Hola, necesito introducir una matriz por bloques de gran dimensión y no se cómo hacerlo de forma sencilla. Hay alguna función o método que pueda utilizar. También quisiera encontrar una forma de generar una variable normal multidimensional de forma más eficiente que con mvrnorm. Gracias.