similar to: INSESGAMIENTO DE ESTIMADORES DE REGRESION

Hola Alexander Quizá esto te pueda ser útil reg<-function(a,b,sigma,x){ data<-a+b*x+rnorm(length(x),0,sigma) lm<-lm(data~x) c(lm$coefficients[1],lm$coefficients[2]) } zz<-replicate(1000,reg(1,2,0.5,100,c(1:100))) hist(zz[1,]) Paco

Hola a todos los R-ayudas estoy tratando de escribir una función en R que me devuelva los estimadores (coef()) de un análisis de regresión no lineal. La estructura de mi base de datos es la siguiente: ##str(nigra) ''data.frame'':   334 obs. of  3 variables:  $ parcela: Factor w/ 55 levels "cu1","cu10","cu11",..: 1 1 1 1 1 1 1 1 12 12 ...  $ edad   :

Hola, Tengo una base de datos con estructura jerárquica en la que quiero clasificar observaciones en distintas categorías. En el caso más simple, tengo una variable con dos categorías (variable Y1) y dentro de cada una de ellas hay otras dos categorías (variable Y2). Además tengo una variable explicativa cuantitativa discreta X. El banco de datos sería de este tipo: X Y1 Y2 5 0 1 9 0 0 2

Hola, buen día. Una consulta, no sé conozco en extenso todas las bondades de R, pero quería saber exactamente si hay algún paquete que permita realizar análisis ecológicos específicamente para comunidades marinas. Existen muchos índices, indicadores, estimadores en general en muchos programas que están a nivel de usuario, pero quería saber si todos estos pueden estar ya disponibles en R.

Saben si el estimador VARHAC de Den Haan y Levin (1994) ha sido implementado en R? Solamente encontre rutinas para RATS y GAUSS en la pagina de Den Haan. Saludos, Jose Jose Luis Iparraguirre [[alternative HTML version deleted]]

Hola:   Tengo un modelo de regresión lineal en el cual las componentes son cosenoidales, y lo construyo del siguiente modo:     modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]" x1[, 1] = cos(2 * pi * t / periods[1]) x2[, 1] = sin(2 * pi * t / periods[1]) for (i in 2:nComp) {   x1[, i] = cos(2 * pi * t / periods[i])   x2[, i] = sin(2 * pi * t / periods[i])   modelo =

Hola buenas tras meses investigando como hacer anovas para factores con efectos fijos y efectos aleatorios, he encontrado una serie de funciones que satisfacen mis pretensiones y creo correctas en cierta medida. Me gustaría compartirlas con vosotros con doble intención, la primera es compartirla para que si otro se encuentra en esta situación que tenga el trabajo hecho y la segunda es que sean

Hi, Does R have any functions implementing such multinomial regression: (S_t^A,S_t^B)~MN(N_t-Y_{t-1},P_t^A,P_t^B) where MN(n,p_1,p_2) is multinomial distribution with parameters n, p_1, p_2. Here P_t^A and P_t^B are nonlinear functions from predictor variables and parameters which need to be estimated. Here A and B are used for notation, they are not parameters. My second question is about

Saludos a todos. Tengo problemas para encontrar la funcion que te permite generar la formula de la regresion lineal entre dos variables y la linea de tendencia en el gráfico de dispersion de puntos, lo que si puedo realizar es el gráfico con los puntos distribuidos en los 2 ejes pero lo demás que mencioné no se como ingresarlo. Saludos, Gabriel Antonio Trujillo Paucar Estudiante de Biologia

Hii all, may I introduce myself. my name is Putro Nugroho i'am astudent of Gadjah Mada University, i'am doing on my final exam about Parametric Regression Model with gamma and log-gamma regression model. I have a problem about the syntax program on R because it doesn't work with gamma distribution with log link function. when I use zelig package it doesn't work when the data is

Hi, I have a dataset which has var1 from 1 sourse and var2 from 2 different methods. I need a new variable such that var2 values from both methods can beused as 1 variable. I believe deming regression can be used to do this. I just don't know how to do it. My data looks like: idvar1var2method1var2method2 11.22.1NA 21.62.4NA 31.52.2NA 41.3NA2.8 51.6NA3.1 61.4NA2.7 Please help. Sincerely,

Hi Hi I am working on credit scoring model using logistic regression. I havd main sample of 42500 clentes and based on their status as regards to defaulted / non - defaulted, I have genereted the probability of default. I have a hold out sample of 5000 clients. I have calculated (1) No of correctly classified goods Gg, (2) No of correcly classified Bads Bg and also (3) number of wrongly

Hi: Could anyone help me to clarify this: are the weights normalized inside lm function (package:stats) before applied to the error term? For example: >lm (cost ~ material, weights=quatity, data=receipt) will lm normalize quatity such that sum(quatity) = 1? I traced to lm.wfit and then the weights get transferred into a precompiled FORTRAN module so I can't figure out. Thanks!

Saludos colegas. Soy nuevo en R y me gustaría saber si es posible embeber código HTML en funciones R para hacer interfaces web. Es decir, lo que quiero hacer es una aplicación web que utilice R para el procesamiento estadístico de los datos y que me los grafique en la web. ¿Es posible hacer esto? . Saludos. --

Dear Mr Peter Dalgaard and Mr Dieter Menne, I sincerely thank you for helping me out with my problem. The thing is taht I already have calculated SENS = Gg / (Gg + Bg) = 89.97% and SPEC = Bb / (Bb + Gb) = 74.38%. Now I have values of SENS and SPEC, which are absolute in nature. My question was how do I interpret these absolue values. How does these values help me to find out wheher my model is

Hola a tod en s, estoy trabajando con modelos de regresión mixta con el paquete "nlme", concretamente con la función "lme". Una vez que hago mi modelo y el summary del mismo obtengo por ejemplo la significación (valor p) de la variable independiente, así como el valor del AIC para comparar mi modelo con otro similar, pero no sé cómo obtener un valor del grado de ajuste de la

Estimados por solicitar su ayuda, estoy intentando cambiar el color de una etiqueta de un gráfico, he realizado múltiples intentos sin un resultado satisfactorio. Debido a que necesito realizar múltiples gráficos he realizado una función que dando parámetros de entrada me grafique las variables o los valores que deseo. Mi problema es en la representación de los colores. Los parámetros de la

While I don't know anything about Box's M test, I googled around and found a Matlab M-file that computes it. Below is my straight-forward translation of the code, without knowing Matlab or the formula (and done in a few minutes). I hope this demonstrates one of Prof. Ripley's point: If you really want to shoot yourself in the foot, you can probably program R to do that for you. [BTW:

Hola a todos...estoy usando MariaDb como base de datos (similares caracteristicas a MySQL), y la conexión entre R y esta con RODBC.¿Existe alguna forma de incrementar la conexión para carga de datos ya que son largos, ejm: 3000 registros en 1 minuto...Si conocen de esto o me pueden sugerir alternativas, estaré agradecido. Patricio. [[alternative HTML version deleted]]

Buenas tardes a todos, Quisiera crear una variable de acuerdo a ciertas condiciones. Me gustaria llegar de "datain" a "dataout". ## entrada datain <- structure(list(REF = c("999", "999", "999", "1099", "731", "731", "731", "731", "1442", "1442", "1442",