similar to: Is there a faster way to do it?

I'm using the parameters "mfrow" and "mfg" to display some graphics(plots) at the same time. Although, because of the parameter "mfg", the parameter "bg" don´t change the background color. What can I do to solve this? Here is the code I'm using... ... i <- 1 for(j in 1:4){ for(k in 1:2){ limiares <-

Estimado Matias, Reenvio a la lista. Es tarde aqui y hay que descansar. Quizas alguien en otra parte del mundo pueda darte una mano. En unas horas te envio mi opinion. Saludos, Jorge.- 2013/6/11 Matias Ledesma <> > Hola Jorge, > > Muchas gracias por las referencias, ya consegui el libro asi que me voy > aponer a leerlo aver si encuentro una solucion. > > La data

hello, i'm pete ,how can i order rows of matrix by max to min value? I have a matrix of membership degrees, with 82 (i) rows and K coloumns, K are clusters. I need first and second largest elements of the i-th row. for example 1 0.66 0.04 0.01 0.30 2 0.02 0.89 0.09 0.00 3 0.06 0.92 0.01 0.01 4 0.07 0.71 0.21 0.01 5 0.10 0.85 0.04 0.01 6 0.91 0.04 0.02 0.02 7 0.00 0.01 0.98 0.00 8 0.02

After ordering the table of membership degrees , i must get the difference between the first and second coloumns , between the first and second largest membership degree of object i. This for K=2,K=3,....to K.max=6. This difference is multiplyed by the Crisp silhouette index vector (si). Too it dependending on K=2,...,K.max=6; the result divided by the sum of these differences I need a final

Estimado Javier, Gracias por tus comentarios. Si, puede ser que me haya expresado en forma incorrecta. Mi idea fue tratar de ser lo mas conciso posible y limitarme al problema estadistico que tengo, por lo cual obvie mucha informacion. La data proviene de un monitoreo que se hace en el Mar Baltico anulamente hace mas de 20 años, debido a que la poblacion a decendido drasticamente en los

Estimado Matías Ledesma Leí en la lista lo que usted desea realizar, en la misma expresa: La idea del estudio es mostrar un valor medio con un intervalo de confianza por estación de la frecuencia de malformaciones por individuo. Creo que se expreso en forma incorrecta, yo no se cuál es su estudio, yo no soy estadístico, soy veterinario, pero lo que usted expresa en su idea, aunque en

Hola Jose, si CONCATENAR significa APILAR, es decir, concantenar verticalmente, por decirlo de algun modo, podrias hacerlo con rbind(): nuevovector <- rbind(vector1,vector2) Si ademas quieres que cada valor de los vectores originales sea identificado en el nuevovector, puedes usar: nuevovector <- stack(vector1,vector2) en este ultimo caso se agrega una columna adicional tipo factor, con

From: matutetote@hotmail.com To: javier.ruben.marcuzzi@gmail.com CC: r-help-es@r-project.org Subject: RE: [R-es] calculo de poder estadistico en glm Date: Wed, 12 Jun 2013 07:24:13 +0000 Estimado Javier, Gracias por tus comentarios. Si, puede ser que me haya expresado en forma incorrecta. Mi idea fue tratar de ser lo mas conciso posible y limitarme al problema estadistico que tengo, por lo

Hola, Si, lo eh pensado pero no eh tenido tiempo de ondar y tengo que profundisar en la biologia y en el procediminto del monitoreo. Dado que todos los años se muestera la misma estacion talvez habria que utilizar un glmm en vez de glm... En realidad lo que estoy aplicando es el modelo de Hurdle, pero haciendolo en dos partes separadas ya que primero analizo la possiblidad que hay de encontrar

Estimado Matías Epidemiología, no es un paquete, me refería a la parte de medicina, que por suerte en R hay herramientas que nos facilitan algunas de estas cuestiones. En http://cran.r-project.org/doc/contrib/Epicalc_Book.pdf hay un libro gratuito de R y epidemiología, aunque hay otros paquetes que tienen algunas cosas que pueden ser mejores, aunque muchas veces eso es a gusto personal o de

por favor como se calculan las probabilidades para conformar el segundo vector SI.Distr, en este ejemplo puse la existente en SARS2003 #Análisis paso a paso rm(list = ls()) #existe información de la incidencia de casos de una enfermedad incidencia <-c (1,0,0,0,0,0,0,0,2,1,1,0,1,2,0,2,0,2,2,2,3, 1, 2, 0,2, 0, 1, 1, 2, 1, 1, 2, 3, 0, 1, 2, 0, 2, 1, 3, 1, 6, 3, 3,7, 4, 4, 3, 6,

Alguien me podría sugerir un paquete en R para generar pronóticos con modelos robusto de series de tiempo.   Saludos Enrique RAMOS [[alternative HTML version deleted]]

Ahh yes, sorry about that. Here's the corrected snippet: # Create an Example Data Frame Containing Car x Color data carnames <- c("bmw","renault","mercedes","seat") carcolors <- c("red","white","silver","green") datavals <- round(rnorm(16, mean=10, sd=4),1) data <- data.frame(Car=rep(carnames,4),

You can use the graphics parameter "srt" to rotate displayed text by a specified number of degrees, e.g. srt=45 to put it on an angle, srt=90 to put it vertical. If you do this, may need to modify the call to text to increase ylim and change the plot location to give you more room. I'm working to update the 'balloonplot' function in the gregmisc package now to handle this

I've just uploaded gregmisc_0.8.2.tar.gz to CRAN. It should show up in the package repository in a day or two. This version of the gregmisc package provides an enhanced 'balloonplot' function with 'rowsrt', 'colsrt' arguments to control rotation of the labels, and 'rowmar', 'colmar' to control the amount of space reserved for the labels. Here's

I have setup a Zentyal 3.2 Samba4 AD. I am using RSAT to administer the AD server. Windows7 and Linux (CentOS) clients are able to join the AD; uid/gid of directories/files show up as per the definitions in Samba4. With Mac OS X clients, there is a problem. The OS X nodes are able to join the domain (quite a few blogs on this subject). The AD connection shows green button. The AD user is

Hi I have some data that i want to plot but i don't find how to do it. I have car types (bmw,renault,mercedes,seat ...), colors and a number for each car type-color relation.I want to come up with a matrix representation of cars vs colors where in each intersection i could set a dot proportional in size to my third variable. Can anybody give me a clue of hoe to come up with such

Normally `codes' is not the appropriate function to use with an unordered factor. Use `unclass' or `as.numeric' to extract the codes used in the internal representation of the factor, as these do not assume that the codes are sorted. and this is one of the `normally' cases. Your code will only work correctly if the levels are in alphabetical order (in the locale

Hi all, This occurred in R-2.11.0 (WinXP). The R-help page of .[ says that: "Character indices can in some circumstances be partially matched (see pmatch) to the names or dimnames of the object being subsetted (but never for subassignment). Unlike S (Becker et al p. 358)), R has never used partial matching when extracting by [, and as from R 2.7.0 partial matching is not by default used by

I want only to ignore the arp request for alias ip ,at the same time I don't want disable the arp function of the interface ? How do ? thanks! --------------------------------- DO YOU YAHOO!? ÑÅ»¢Ãâ·ÑGÓÊÏ䣭ÖйúµÚÒ»¾øÎÞÀ¬»øÓʼþɧÈų¬´óÓÊÏä From vova at fbsd.ru Fri Jul 1 13:12:10 2005 From: vova at fbsd.ru (Vladimir Grebenschikov) Date: Fri Jul 1 13:12:12 2005 Subject: how to ignore