Buenas tardes.  Quisiera su ayuda con el siguiente problema.

Debo ajustar un conjunto de datos a una ley de potencias cuya forma es f(x)=C(k) x^{-k}, donde C(k) se obtiene a partir de la condición de completez, es decir, que la integral de f(x) es igual a 1.  El problema es que C(k) depende del parámetro K y ese parámetro se obtiene a partir del método de MLE.

¿Cómo puedo definir en R una forma específica para una función de distribución (nueva) y obtener su MLE para k?

NOTA: es visto que, al parecer, se puede aproximar mi ley de potencias a una Distribución de Pareto pero no veo cómo.

Les agradezco mucho su ayuda

Hola tengo un problema con unos datos ya me pude conectar a mongodb con Rstudio  usando la libreria mongolite. de ahi tomo los datos  la columna date  viene asi 01/14/2013 07:00:00 AM"  quiero saber que día de la semana fue (Lunes, martes, Miercoles, jueves, viernes,  sabado ,domingo) con la libreria openair  con la funcion CutData se supone que me dice que dia fue  solo que al tratar de ejecutarla me marca ese error. 

pues tengo unas funciones que se llaman:
sheet.x1,sheet.x2,sheet.x3 y sheet.x4

Lo que quiero hacer es un bucle que me llame a esas funciones ejemplo

for(i in 1:4){
sheet.x||i
}

Pero no se la manera de concatenar eso, en phyton seria

for(i in 1:4){
sheet.x${i}
}

me podeis echar una mano? gracias.

Pues veamos tengo un excel de donde cargo unos datos de 3 años y genero una matrid con esos datos y luego le agrego una prediccion de arima y exponencial.

Me gustaría añadir en ese mismo excel despues de las predicciones 3 columnas apartes que contengan los datos del primer año , del segundo y del 3 año.(de los datos originales.)

Esta es la parte del codigo donde unifico la prediccion y la secuencia original. (la seguencia original se la entrego en una matriz de 1columnax42filas.

no.pred   <- rep(NA, times = pred.times)

data.pred <- data.frame(no.pred, pred.arima.prec, pred.arima.auto, 

                        pred.HW.add.beta.gamma, pred.HW.add.gamma,

                        pred.HW.mul.beta.gamma, pred.HW.mul.gamma,

                        pred.HW.no.seasonal, pred.HW.exponential.smooth)

colnames(data.pred) <- c(series.name, name.arima.prec, name.arima.auto, 

                         name.HW.add.beta.gamma, name.HW.add.gamma,

                         name.HW.mul.beta.gamma, name.HW.mul.gamma,

                         name.HW.no.seasonal, name.HW.exponential.smooth)

data.original <- data.frame(rep(values, times = length(colnames(data.pred))))

colnames(data.original) <- colnames(data.pred)

data.complete <- rbind(data.original, data.pred)

first.date    <- as.Date(paste(year.begin, month.begin, '1', sep = '/'))

dates         <- as.character(seq(from = first.date, by = 'month', length.out = nrow(data.complete)))

data.complete <- data.frame(dates, data.complete)

Hola soy nueva en R, tengo un csv que contiene información de películas dos de estas columnas son director_name  y Genero(Estos están divididos por este caracter | algunos tienen un solo genero y otros varían) quiero obtener el director que mas se repita y el genero que mas se repita 

Alguien que me ayude con diseños geneticos en R, o como calcular los componenes de varianza.

gracias

Hola, soy nuevo con R y me gustaria saber como puedo crear una gráfica a partir de unos datos en un archivo CSV. Primero subo los datos y a continuación pongo la siguiente instruccion: Protein<-gvisBubbleChart(proteins, idvar="prote", xvar="datandb",yvar="name", colorvar="measure", sizevar="datanutrients_id", options=list(hAxis='(minValue:75,maxValue:150)')

+Plot(protein)

Pero me da error.

saludos,

Soy  nuevo en el mundo de R y quería saber si me pueden ayudar,

quiero ejecutar un script para eliminar outliers el escript lo descargue de la siguiente web http://datascienceplus.com/identify-describe-plot-and-removing-the-outliers-from-the-dataset/ ahí especifica que debo reemplazar "dat" con el nombre de mi base de datos, y "var" con el nombre de mi variable pero no se a que se refiere con eso xq no tengo ninguna variable si me pudieran ayudar de antemano muchas gracias.

Hola,

Quise  instalar el paquete summary.gamlss pero sale el siguiente mensaje: is not available (for R version 3.3.1)

¿Cómo puedo saber cuál es la versión adecuada de R para instalar este paquete?

Gracias por la ayuda

Hola muy buenas, estoy haciendo un trabajo en R y me piden lo siguiente...

En una mesa de dados de un casino se echa un par de dados no trucados 180 veces por hora.

En 24 horas se tiran los dados 4320.

Los parámetros de la distribución que me dan son (4320, 1/6)

Calcular con ayuda de R la probabilidad de que haya entre 700 y 750  sumas de 7 en 24 horas.

¿Cómo puede calcular esta probabilidad en R?

Un saludo y muchas gracias.