similar to: Función ARIMA

Buenas tardes, ¿Cómo interpretarías el intercepto que da R en un modelo de ceros inflados? Por un lado en la parte de conteo tengo un intercepto de -4.2 y en la parte de ceros de 102, ambos salen significativos (***). ¿Qué me dirían? Gracias

Hola a todas, Quiero realizar un subset selection usando el paquete leaps, entre mis variables explicativas tengo rezagos de las mismas, por tanto tienen datos NA. sin embargo, al tratar de realizar cross validation me pide que los datos esten en formato data.frame con lo que me arroja un error al ejecutar el algoritmo. Mi pregunta es cómo puedo forzar en un data.frame en que tengo variables

Muy buenas, quiero calcular correlaciones de Pearson entre dos variables (a,b) teniendo en cuenta una tercera (c). Para ello estoy usando una función llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce una función similar en alguna librería de R? Por otro lado, para ver si me cuadraban los resultados,

Buenos días, *Muchas gracias, todas las aportaciones han sido bien útiles.* Las he tenido en cuenta y he pasado los datos con el R, siguiendo el siguiente comando: *modeloRTUN2<-lm(AVE.~ Tariff + d1 + d2 + d3 + d4 + d5 + d6 + d7 + d8 + d9+ d10 + d11+ d12+ d13+ d14+ d15+ d16+ d17+ d18+ d19+ d20 +d21 + Tariff*d1 + Tariff*d2 + Tariff*d10)* *summary(modeloRTUN2)* Siendo: AVE. = Variable

Hola a todos, quería preguntaros un medio para obtener los valores p usando lmer. He tratado con pvals.fnc, que es lo que me habían recomendado, pero por algún motivo no está ya disponible etc. Ésta es la función que tengo, pero da las "t", sin los valores p. Aunque Baayen indica que valores por encima de 2 son significativos necesito saber las p. resultado = lmer(rt_ln ~ (fre_ln *

I have problem with my shorewall. We are now doing some stress test with a http application behind the shorewall. Firstly we send 10.000 requests to a http based application with no firewall. It can received 100% requests. But when we put shorewall in front of it then it stats to loose requests. Is there any packet limitation from shorewall all it''s about conntrack? Thanks for the reply.

Hi, I would like to save many trees created in a loop as different ones. I can plot them, but I can not get the whole tree properties saved. I paste the script below Also, How can I perform multivariate trees?, that is predicting a vector o values, not a simple scalar observation. Thanks in advance Angel for (i in 1:7){#loop para hacer arb arb=arb[i]#contador

Hello Dovecoters, I am using dovecot-2.2.10-8.el7.x86_64 on RHEL 7.4 When I backup mdboxes (using SIS) using this command doveadm -Dv -o plugin/quota= sync -1f -u user1 at user.com mdbox:/mailbackup/user1 mails get backedup but all m files are not copied. Last file copied is m.4 If I run the command again I get below error I tried with doveadm backup command as well (many times).

Hi everybody. I have an HP StorageWorks array connected to one of my servers, connected via an LSI controller (lspci reports: SCSI storage controller: LSI Logic / Symbios Logic SAS1068E PCI-Express Fusion-MPT SAS (rev 08)) on an HP DL385 running CentOS 5.5 x86_64. I have set up multipath and all is working but, when I create a snaphost and export it as a new LUN from the HP web interface, it

Hola Gente, Como están? necesito ayuda con modelos aditivos generalizados, GAMs. Estou usando el paquete mgcv. tengo una tabla de datos con lagunas y pozos de erosión, cada ambiente con 6 muestreos no tengo problema en disponibilidad los datos para que me ayuden. a variable y= Distmod que es el resultado de distancia de una matriz de especie x= Tpico: y es el tiempo que ocurrió la

Dear all, We're trying to use Nessus to monitor our system and we're having a hard time mapping the package versions from RedHat to CentOS in order to be able to report against CVE (vulnerability reports). After some research, we think that the mapping is as follows (using HTTPD as the example): CentOS httd-2.2.3-43.el5.centos.3 is equivalent to RedHat httpd-2.2.3-43.el5_5.3 So, it

Hola Manuel lo he tratado de hacer pero me sale Error: unexpected string constante in: "anova(a,as,test=Chisq") no tengo ni idea de por qué... Me resulta alucinante no poder contar ya con pvals.fnc. ¿Será imposible hacerse con ello? Saludos, Miguel -------------------------------------------- El vie, 13/6/14, Manuel Azcárate <mazcarategarcia en gmail.com> escribió:

Thank you all!!! The problem was the decimal symbol! My data was saved in a txt file, so I?ve introduced the dec="," in ?read.table? and it worked. What I?ve done was coeficientes<-read.table("coeficientes.txt",sep="\t",header=T,dec=",") Then, subset worked fine coeficientesWanted<-subset(coeficientes,b1>0) Thanks again, Cec?lia Carmo

Hi R-helpers! I have the following object: > head(coeficientes) caedois b1 b2 b3 1 1 0,033120395 -20,29478338 -0,274638864 2 2 -0,040629634 74,54239889 -0,069958424 3 5 -0,001116816 35,2398622 0,214327185 4 10 0,171875 5 14 0,007288399 40,06560548 -0,081828338 6 15 0,027530346 0,969969409 0,102775555

--- the forwarded message follows --- -------------- next part -------------- An embedded message was scrubbed... From: "Cecilia Carmo" <cecilia.carmo at ua.pt> Subject: Re: [R] subset dataframe/list Date: Mon, 01 Jun 2009 21:33:15 +0100 Size: 3657 URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090601/921f7638/attachment-0002.mht>

Estimad en s errer en s He hecho este modelo glm m1.pile<-glm(ger~tem+pot+time+I(tem^2)+I(tem^2):pot ,family="binomial" ,data=long.PILE ) Que nos da la probabilidad de germinación de una semilla en función de tem (Temperatura), pot (Humedad del suelo) y time (Tiempo que la semilla pasa en esas condiciones). Ahora quiero, para diferentes tem, pot

Un ejemplo con un modelo más simple: He especificado este modelo: >formula(m2.pile) ger ~ tem + pot + time Si hago predict me da: >predict(m2.pile,newdata=data.frame(tem=25,pot=0,time=3),type="response") 0.08243262 Extraigo los coeficientes: > coef(m2.pile) (Intercept) tem pot time -1.89521331 -0.02303313 4.74499714 0.02043222 Ahora calculo la

¿Quieres decir que para un nivel de una variable categorica todas las observaciones de la variable respuesta sean ceros? Gracias El Lun, 9 de Abril de 2018, 19:59, Carlos J. Gil Bellosta escribió: > ¿Podría ser que para algún nivel de alguna variable independiente > categórica solo hubiese ceros? En ese caso, casi seguro, aparecería ese > tipo de warning. > > El lun., 9 abr. 2018 a

Buenas tardes, Estoy estimando un modelo binomial negativo de ceros inflados (ZINB) utilizando el comando zeroinfl() del paquete pscl. Al ejecutarlo me da el siguiente aviso: Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred ¿Sabéis que significa y si puedo usar el modelo aún con ese aviso? ¿Los coeficientes son fiables? Muchas gracias, Miriam

Muchas gracias por la respuesta. He mirado y los coeficientes no son altos pero sí tengo una gran cantidad de ceros en la variable dependiente (más del 90%). Sin embargo, al incluir otro tipo de variables independientes no me da ese aviso, dejando la misma variable dependiente. ¿Cómo podría utilizar stan/rstan de forma sencilla para diagnosticar el modelo? Muchas gracias El Lun, 9 de Abril de