similar to: R help for a newby

Dear R users, I am trying to produce some plots in a postscript file, but I am experiencing some issues. I open the device with ----------------------------- setPS() postscript (file='gs_mcmc_dust.ps',width=5*3,height=5*3,horizontal = FALSE, paper = "special",family = "ComputerModern",encoding="TeXtext.enc")#, onefile = FALSE -----------------------------

Hola eRReros. Os lo explico con un ejemplo: # Cargamos los paquetes y el ejemplo install.packages("SPEI") library("SPEI") data(wichita) # los primeros 12 datos head(wichita,12) # mi subset de los primeros 12 datos meu<-wichita[c(1:12),] meu # como veis los valores de TMED son iguales en ambos dataframes. # ahora viene el problema # calculamos

I found the need of converting a matrix into a list of its columns (for use with do.call), and was surprised there was no method as.list.matrix, it could easily be a part of as.list default I wrote my.as.list.matrix <- function(mat) { if(!is.matrix(mat))stop("Argument must be a matrix") n <- NCOL(mat) res <- vector(mode="list", length=n)

Script started on Tue Mar 3 07:35:34 2015 doctor.nl2k.ab.ca//usr/source/openssh-SNAP-20150304$ make tests [ -d `pwd`/regress ] || mkdir -p `pwd`/regress [ -d `pwd`/regress/unittests ] || mkdir -p `pwd`/regress/unittests [ -d `pwd`/regress/unittests/test_helper ] || mkdir -p `pwd`/regress/unittests/test_helper [ -d `pwd`/regress/unittests/sshbuf ] || mkdir -p `pwd`/regress/unittests/sshbuf [ -d

Estimada comunidad tengo un problema del que no encuentro datos que me ayuden mucho en la web. Estoy haciendo calculos por grupos con data,table. Tengo un archivo (zp.res) con tres columnas que clasifican los datos (sol, con, dia) y una columna de datos numericos (media), de la siguiente forma: sol con dia media 1: con 0 1 -22.6 2: con 0 1 -36.6 3: con 0 1 -35.6 y

Hi R Users, I am going to run a multiple linear regression with around 57 independent variables. Each time I run the model with just 11 variables, the results are reasonable. With increasing the number of independent variables more than 11, the coefficients will get ?NA? in the output. Is there any limitation for the number of independent variables in multiple linear regressions in R? I attached

Dear Community, I have been using two types of survival programs to analyse a data set. The first one is an R function called aftreg. The second one an STATA function called streg. Both of them include the same analyisis with a weibull distribution. Yet, results are very different. Shouldn't the results be the same? Kind regards, J -- View this message in context:

Dear All, I am newbie to R, and I wanted to plot a barplots with R and in such a way that It will also show me position which I can plot on the bar line. Here is my code that I am using to plot, > chromosome <- c(40.2, 35.6, 36.1, 29.6, 31, 29.6, 31, 29.4, 28.2, 23, 23, 28.2) >barplot (chromosome, col="purple", xlab="Oryza sativa Chromosomes", border = NA, space =

Hello, I am not entirely sure the subject line captures what I am trying to do, but hopefully this description of the problem will help folks to see my challenge and hopefully offer constructive assistance. I have an experimental setup where I measure the decrease in oxygen in small vials as an organism, such as an oyster, consumes the oxygen. Each vial is calibrated before the experiment and

Short: get rid of the loops I use and optimize runtime Dear all, I want to calculate for each row the amount of the month ago. I use a matrix with 2100 rows and 22 colums (which is still a very small matrix. nrows of other matrixes can easily be more then 100000) Table before Year month quarter yearmonth Service ... Amount 2009 9 Q3 092009 A ...

We''ve finally finished up work on our application and are now beginning our beta. The site (SMS+mobile blogging+event planning+LBS+photos) is built using 100% Ruby on Rails. Many thanks to all the people who answered questions over the months on this list. Our experience with Rails was overwhelmingly positive! We''d love to hear any feedback you might have:

Hi, I want to remove the entire row from the dataset if any of the row contains blank or NA my dataset look like State Year Y X2 X3 X4 X5 X6 State1 1960 27.8 397.5 42.2 50.7 78.3 65.8 State2 1960 29.9 413.3 38.1 52 79.2 66.9 State1 1961 29.8 439.2 40.3 54 79.2 67.8 State2 1961 30.8 459.7 39.5 79.2 69.6 State1 1962 31.2 492.9 37.3 54.7 77.4 68.7 State2 1962 528.6 38.1 80.2 73.6 State1 1963

Hi all, sorry to intrude on xen-devel, but I think I need direction from the expertise here. I''ve admin''d Xen servers of various flavors for a couple years, but never seen this before. After a period ranging from several hours to several days, my primary database and development DomU completely locks up. Net disconnects, but CPU(sec) continues to tick in xentop. No errors, and

Hello, I''m getting an URI::InvalidURIError in my controller when I use german characters in the uri. Google accepts them but open-uri not. How to encode them? That''s my controller simplified: def index require ''open-uri'' address = "Bürgerstraße+68+01127+Dresden" api_key = "my api key" json =

Dear All: Urgent help is needed. I have a data set in matrix format of three columns: X, Y and index of four groups (1,2,3,4). What I need to do is the following; 1- How I can subtract the sample mean of each group indexed 1,2,3,4 from the corresponding data values of this group and create new columns say X-sample mean and Y-sample mean? I tried to use the "tapply" but

I asked this before and I am going to try again in more applied terms. I am trying to use R to extract variance components for a two-factor random effects model with both factors crossed. It would also be nice to generate some confidence intervals as well. For example, a data set using SAS Proc Mixed is below followed by the four variance component estimates and the respective confidence

I have noticed a significant performance degradation using merge in 2.9.1 relative to 2.8.1. Here is what I observed: N <- 100000 X <- data.frame(group=rep(12:1, each=N), mon=rep(rev(month.abb), each=N)) X$mon <- as.character(X$mon) Y <- data.frame(mon=month.abb, letter=letters[1:12]) Y$mon <- as.character(Y$mon) Z <- cbind(Y, group=1:12) system.time(Out

I am hoping for some advie regarding the following scenario. I have data from 26 studies that I wanted to analyse using logistic regression. Given the data was from studies and not individuals I was unsure how I would perform this in R. When analysed in SPSS, weighting was used so that each study was included twice. Where "use" occurred, a value of 1 was assigned and was weighted

I have been trying to run a logistic regression using a number of studies. Below is the syntax, error message & data. Any advice regarding what I am doing wrong or solutions are appreciated, regards Bob Green > logreg <- read.csv("c:\\logregtest.csv",header=T) > attach(logreg) > names(logreg) [1] "medyear" "where" "who"

Hello, dear R-community. This is a question about TukeyHSD between factor combinations of a Three-Way ANOVA, which is - since it is a multi measure ANOVA - not a simple ANOVA but a Generalised Linear Mixed Model (GLMM), calculated with "lmer". > growth <- groupedData(length~meas|box_id,outer=~spec*comp*water,data=all.spec) > model <-