[rkward-devel] FALLA en RKward
Juan Manuel Piñeiro
jpineiro at agro.uba.ar
Mon Mar 18 13:18:25 UTC 2013
El motor R se apagará inmediatamente. Esto significa que no podrá usar más
las funciones que dependen de él. Es decir, no podrá hacer prácticamente
nada, ni siquiera guardar el espacio de trabajo (pero, si tiene suerte, R
ya lo hizo). Lo que sí puede hacer es guardar cualquier archivo abierto, la
salida, o copiar datos de los editores abiertos. Después salga de RKWard.
Dado que esto nunca debería pasar, agradeceríamos que nos escribiese a
rkward-devel at lists.sourceforge.net contando lo que intentaba hacer cuando
esto ocurrió. ¡Lo sentimos!
no terminó de correr el siguiente script:
M1<-lme(ivncom~scom+scom2+bm+bm2+area+area2+areaperimetro+areaperimetro2+desnivel+desnivel2+hete+hete2+even+even2+scom.even+scom.even.2+scom.2.even+scom.2.even.2+com,data
= base_com, random = ~1 |fecha/ut, method= "ML", na.action=(na.omit)) # ut
is nested in fecha
d1<-dredge(M1, fixed="scom", evaluate=TRUE, rank=AICc,beta=TRUE)
babscom<-subset(d1,which(abs(d1$AICc)<(min(abs(d1$AICc))+2))) # devuelve el
minimo del criterio y los que estan entre #ese y + 2 unidades del criterio
bestscom<-subset(d1,which(d1$AICc<(min(d1$AICc)+2))) # devuelve el minimo
del criterio y los que estan entre ese y + 2 #unidades del criterio
write.table(d1,"MUMINdredge1.txt",row.names=FALSE)
write.table(babscom,"MUMINbabscom.txt",row.names=FALSE)
write.table(bestscom,"MUMINbestscom.txt",row.names=FALSE)
#escala marco
M2<-lme(ivnum~sm+sm2+bm+bm2+area+area2+areaperimetro+areaperimetro2+desnivel+desnivel2+hete+hete2,data
= base_um, random = ~1 |fecha/ut/com, method= "ML", na.action=(na.omit)) #
ut is nested in fecha
d2<-dredge(M2, fixed="sm", evaluate=TRUE, rank=AICc,beta=TRUE)
babsum<-subset(d2,which(abs(d2$AICc)<(min(abs(d2$AICc))+2))) # devuelve el
minimo EN VALOR ABSOLUTO del criterio #y los que estan entre ese nº + 2
unidades del criterio
bestsum<-subset(d2,which(d2$AICc<(min(d2$AICc)+2))) # devuelve el minimo
del criterio y los que estan entre ese y + 2 #unidades del criterio
write.table(d2,"MUMINdredge2.txt",row.names=FALSE)
write.table(babsum,"MUMINbabsum.txt",row.names=FALSE)
write.table(bestsum,"MUMINbestsum.txt",row.names=FALSE)
## prueba de grafico de un modelo
bo<-mix_lin9(
base_com[,"scom"],base_com[,"com"],base_com[,"even"],base_com[,"ivncom"],base_com[,"fecha"],base_com[,"ut"],
grafica=TRUE)
local({
## Preparar
## Computar
dev.set (2)
dev.print (device=postscript, file="bo.eps", onefile=FALSE,
paper="special", horizontal=FALSE);
## Imprimir el resultado
})
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.kde.org/pipermail/rkward-devel/attachments/20130318/cda871e9/attachment.html>
More information about the Rkward-devel
mailing list