similar to: ENMCA in EnQuireR problems!

Hi, I am trying to use WinBUGS using R, but i am having some problems using bugs function: out <- bugs(data = win.data, inits = inits, parameters.to.save = params, model.file = "model.txt", n.thin = nt, n.chains = nc, n.burnin = nb, n.iter = ni, debug = TRUE, DIC = TRUE, working.directory = getwd()) ls() WinBUGS file is in working directory=D:/Line/documents but R looks for

I would like to set the y axis limit of an xyplot using the object 'ylimit', but receive this error: [1] 990 Error in extend.limits(limitlist[[i]], axs = axs) : improper length of lim I get the same error if I use ylim. library(lattice) trellis.device(col = FALSE, theme = lattice.getOption("col.whitebg")) name <- "Variable name" symbols <-

Hi all, I've got a density graph made with the following commands: win.graph(width=13,height=6) par ( fin=c(13,3) ,mai=c(1,1,0.5,0.5) ,mfrow=c(1,2) ,cex.axis=1.5 ,cex.lab=1.5) dens<-density(DATA1.y[2,]-mean(DATA1.y[2,]),kernel="gaussian") xlimit<-range(dens$x) ylimit<-range(dens$y) hist( DATA1.y[2,]-mean(DATA1.y[2,]) ,xlim=1.1*xlimit ,xlab=expression(q[e])

Hi all, Please let me know if this isn't the correct place to report bugs in contributed code. Otherwise, the following code contains one additional line and a short comment above it. The extra line of code catches the unusual situation where a data column contains all 0 values. The corrected code is below the signature. Thanks for your help. Tom Thomas S. Dye, Ph.D.

Dear R users, I m quite a novice in using R for factor analysis and I would need some help to choose the right function. I have a contingency table and I would like to perform a Correspondence analysis on this table, followed by a hirarchical clustering of my variables projected in on the first principal components. Here are my question : - what is the more appropriate function to do so ...

Hi, I am Using R.2.9.2 and I am trying to make nicer some graphics, which are constructed with plotmer.fnc. In particular, I have two questions: 1. Is is a way to decide in which order we want the labels on the x-axis. For now, the reference value label is the first one, and the other labels appear in the alphabetical order. I would like them to be sorted in the ascending way. Here is my

Thanks for the insight. > -----Original Message----- > From: Prof Brian Ripley [mailto:ripley at stats.ox.ac.uk] <snip> > dots <- list(...) > haveYlim <- "ylim" %in% names(dots) > > is the sort of thing we still understand 5 years later. > I didn't say "understand", I said "easily follow". Obviously how "easily" is

Hi, I have the following file, and I need to work out the linear regression for each sample. I tried the model(*) and receive the error message (**): > data=split(mydata,rep(1:(nrow(mydata)/6),each=6)) > arrang.linear=lapply(data,lm,formula=KA~PA) Erro em storage.mode(y) <- "double" : invalid to change the storage mode of a factor Al?m disso: Warning message: In

Se não pode visualizar este e-mail corretamente, clique aqui <http://www.almasoma.pt/newsletter--terapia-para-bebes-aprenda-a-ajuda-los> . <http://www.almasoma.pt/> Formação em Terapia para Bebés Lisboa, Outono 2013 - Primavera 2016 Workshop: "Baby Clinic", 6 e 7 de Abril de 2013 Conferência de Apresentação, 5 de Abril de 2013 Karlton Terry Pioneiro da assistência aos

Dear R helpers, I am quite new working with R, and that might be a very easy question for some of you. So I hope you can help me!! I used a hirarchical cluster analysis to create dendrograms, like that: hc <- hclust(dist(USArrests), "ave") plot(hc) plot(hc, hang = -1) Here I want to compare the number of bifurcations for each species, but therefor I need this numbers first. There

Good day all, I'm trying to plot a figure and ensure that the minimum and maximum values are always displayed. However, the code below does not display the minimum value, no matter what I try. Could someone please help? Thanking you in advance, George. Code below for reproduction (apologies if the paste is too long) >blah [1] 0.41955516 0.37330949 0.37934005 0.38013805 0.40092939

Consider the following: library(lattice) library(latticeExtra) temp <- expand.grid( subject = factor(paste('Subject', 1:3)), var = factor(paste('Variable', 1:3)), time = 1:10 ) temp$resp <- rnorm(nrow(temp), 10 * as.numeric(temp$var), 1) ylimits <- by(temp$resp, temp$var, function(x) range(pretty(x))) useOuterStrips(xyplot( resp ~ time

I wrote a little function called bgfun that adds gridlines and a background, but it's not working with I plot using the formula. I have some theories on what's happening, but even if my theory is right, I don't know how to fix it. Someone have a straightforward silver bullet? Thank you, Gene bgfun = function(color='honeydew2',linecolor='grey45',

On Tue, 2016-05-17 at 20:12 -0400, Jonathan Billings wrote: > On May 17, 2016, at 7:56 PM, Always Learning <centos at u68.u22.net> wrote: > > (1) I would change the port from 22 to something more difficult to > > guess, perhaps 49026 (for example) and then block port 22 in the > > firewall. > > If you?re going to change the port, change it to something <1024.

Hello!   I am trying to merge two xy-plot with different ylimits.   Can someone please give me possible way of achieving it..   Thank you Regards Dimple [[alternative HTML version deleted]]

Hi all, I'm trying to log chart but with natural looking tick marks. My specifications are very specific -- it must indicate the lowest number's tick as well as the maximum. I've attached sample code and data for a particular case (and there are a few more like this) where the bottom tickmarks on the chart are not set to where I want them to be and yet they fit in the ylim range.

Hi everyone! I'm digging into ggplot for some while now, I must say it's great! But - as some others have posted before, and hadley knows very well himself - the documentation is lacking some bits... So I have to pose that question directly: I'd like to produce a series of maps with different data on, but exactly the same extent in each plot. Is there a way of switching the

Hi R users, Does anybody knows for the following erro after running geweke.diag(MCMC.sampled, frac1=0.1, frac2=0.5) Erro em glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : la??o interno 1; n??o ?? poss??vel corrigir o tamanho do passo Al??m disso: Warning messages: 1: algoritmo n??o convergiu in: glm.fit(x = X, y = Y, weights = weights, start = start,

Hi, I make this model using lme m.lme <- lme(Glycogen~Treatment,random=~1|rTrt/Liver) How to make this using lmer? I try > m.lmer <- lmer(Glycogen~Treatment+(1|rTrt/Liver)) Erro em lmer(Glycogen ~ Treatment + (1 | rTrt/Liver)) : entry 0 in matrix[0,0] has row 2147483647 and column 2147483647 Al??m disso: Mensagem de aviso: / not meaningful for factors in: Ops.factor(rTrt, Liver)

Dear useR's, I wanted to load the package GeoXP in my computer but I couldn't because this warning message appeared: > require(GeoXp) Carregando pacotes exigidos: GeoXp Carregando pacotes exigidos: tcltk Loading Tcl/Tk interface ... Erro em fun(...) : couldn't connect to display ":0" Al?m disso: Warning message: pacote 'GeoXp' foi compilado na vers?o do R 2.4.1