<style type="text/css">
p, li { white-space: pre-wrap; }
</style>
<p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px">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.</p>
<p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px">Dado que esto nunca debería pasar, agradeceríamos que nos escribiese a <a href="mailto:rkward-devel@lists.sourceforge.net">rkward-devel@lists.sourceforge.net</a> contando lo que intentaba hacer cuando esto ocurrió. ¡Lo sentimos!</p>
<p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px"><br></p><p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px">no terminó de correr el siguiente script: <br>
</p><p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px"><br></p><p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px">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</p>
<p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px">d1<-dredge(M1, fixed="scom", evaluate=TRUE, rank=AICc,beta=TRUE)<br><br>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<br>
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<br>write.table(d1,"MUMINdredge1.txt",row.names=FALSE)<br>write.table(babscom,"MUMINbabscom.txt",row.names=FALSE)<br>
write.table(bestscom,"MUMINbestscom.txt",row.names=FALSE)<br><br>#escala marco<br>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<br>
d2<-dredge(M2, fixed="sm", evaluate=TRUE, rank=AICc,beta=TRUE)<br>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<br>
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<br>write.table(d2,"MUMINdredge2.txt",row.names=FALSE)<br>write.table(babsum,"MUMINbabsum.txt",row.names=FALSE)<br>
write.table(bestsum,"MUMINbestsum.txt",row.names=FALSE)<br></p><p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px"><br></p><p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px">
## prueba de grafico de un modelo<br>bo<-mix_lin9( base_com[,"scom"],base_com[,"com"],base_com[,"even"],base_com[,"ivncom"],base_com[,"fecha"],base_com[,"ut"], grafica=TRUE)<br>
local({<br>## Preparar<br>## Computar<br>dev.set (2)<br>dev.print (device=postscript, file="bo.eps", onefile=FALSE, paper="special", horizontal=FALSE);<br>## Imprimir el resultado<br>})<br><br></p>