similar to: Convert string to list?

Hello list members, I'm looking for a way to force the results of a crosstabulation to be square - that is, to include 0 values. For example: table(letters[1:4],letters[c(1:3,3)]) yields: a b c a 1 0 0 b 0 1 0 c 0 0 1 d 0 0 1 I would like to return: a b c d a 1 0 0 0 b 0 1 0 0 c 0 0 1 0 d 0 0 1 0 Any suggestions? Thanks! -- Manuel A. Morales

Hello list, I'm trying to fit a model like beta[trt]/(1+alpha*x) where the data include some grouping factor. The problem is that the estimate for alpha is undefined for some of the treatments - any value greater than 20 is equally good and a step function would suffice. Ignoring the grouping structure, I can fit this using nls with the port algorithm by restricting the upper value of alpha

I tried to submit this in R, but not sure if it worked. When running R as sudo, using ctrl-c dumps me to the command line. Hitting exit to exit the terminal window results in R taking 100% of resources. I am using R-2.1.0 on Fedora Core 3. Thanks. Manuel

Hi, I would like to have lineplot.CI and barplot.CI to actually plot confidence intervals , instead of standard error. I understand I have to use the ci.fun option, but I'm not quite sure how. Like this : > qt(0.975,df=n-1)*s/sqrt(n) but how can I apply it to visualize the length of the student's T confidence intervals rather than the stdandard error of the plotted means

Is it possible to add standard error bars to the means on interaction plots? Thanks Murray -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

Dear r-helpers, Spencer Graves and Manual Morales proposed the following methods to simulate p-values in lme4: ************preliminary************ require(lme4) require(MASS) summary(glm(y ~ lbase*trt + lage + V4, family = poisson, data = epil), cor = FALSE) epil2 <- epil[epil$period == 1, ] epil2["period"] <- rep(0, 59); epil2["y"] <- epil2["base"]

Hello: I'm having quite a bit of difficulty adjusting the x-labels using bargraphCI. I've tried using text and srt=45 to rotate the labels or mtext for 2 lines to break up the labels. However, using either method, I cannot line up the labels with the midpoints of the bars (they line up with some sort of tick mark that is off the midpoint of the bars). Any suggestions would be greatly

Hi all, Given a data frame as... > head(veg) genus species trophia type geo zone importance 1 Sphagnum subsecundum M A En 100 2 Sphagnum denticulatum M A En 200 3 Molinia caerulea M A En 300 4 Sphagnum flexuosum M A En 400 5 Juncus squarrosus M A En 500 6

Full_Name: Manuel A. Morales Version: 2.8.1 OS: Linux Submission from: (NULL) (137.165.199.246) When fitting a model in nls using the algorithm port with constraints and the shorthand parameter[factor] in the model, I get the following error message: "Error in nls_port_fit(m, start, lower, upper, control, trace) : (list) object cannot be coerced to type 'double' In addition:

Dear list, A while ago, I posted a question asking how to use data or subset arguments in a user-defined function. Duncan Murdoch suggested the following solution in the context of a data argument: data <- data.frame(a=c(1:10),b=c(1:10)) eg.fn <- function(expr, data) { x <- eval(substitute(expr), envir=data) return(mean(x)) } eg.fn(a,data) I've

Hello list, Is anyone aware of a non-lattice-based alternative to xyplot()? Thanks! Manuel -- http://mutualism.williams.edu

Hi list, I'd like to install an archived version of lmer to compare results with the current version. I guess one way to do this would be to download the source, rename the package and then install it. Is there a better alternative? -- http://mutualism.williams.edu -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type:

Hello list, I'm doing a bootstrap analysis where some models occasionally fail to converge. I'd like to automate the process of restricting AIC to the models that do converge. A contrived example of what I'd like to do is below: resp <- c(1,1,2) pred <- c(1,2,3) m1 <- lm(resp~pred) m2 <- lm(resp~poly(pred,2)) m3 <- lm(resp~poly(pred,3)) # Fails, obviously ## Some

Dear R users, I have now tried out several options of obtaining p-values for (quasi)poisson lmer models, including Markov-chain Monte Carlo sampling and single-term deletions with subsequent chi-square tests (although I am aware that the latter may be problematic). However, I encountered several problems that can be classified as (1) the quasipoisson lmer model does not give p-values when

I'm trying to create a bargraph of means with standard error bars using the function bargraph.CI (in the sciplot package). Like this: bargraph.CI(x.factor, response,data,xlab, ylab, par(family="serif"),font=11) However, an error message comes up when I try to use the par funtion. Does the character font/style need to be changed in a different way for this kind of plot? Or am I

Dear R-lister I have a question about how to create multiple error bar plots. I found that I can use an "errbar" function from Hmisc package to create one error bar plot in which there are multiple data points (x, y) with the error bars. Thus, I know that I can get "one" error bar plot which consists of many data points. However, I did not find a way to put

The package "sciplot" is now available for download from CRAN. This package includes a collection of functions that create graphs with error bars for data collected from one-way or higher factorial designs, as well as a function to plot bifurcation diagrams resulting from analysis with XPPAUTO. The functions in this package replicate some of the functionality of plotmeans() from the

Hola, estoy representando los datos Duncan de la libreria car y cuando hago un gráfico lattice observo que la profesión "architect" aparece en dos subgráficos que no se corresponden con su categoría. Alguien sabe el motivo? Os adjunto un fichero con el código y el gráfico resultante. Gracias, Manuel [[alternative HTML version deleted]]

Hi List, I have a problem plotting data using the lineplot.CI command in the sciplot package. I want to plot the data of 2 experimental cases using different lines (traces). Time is on the X-axis. The tricky thing is that the data collection in the second case started later than for the first case. This is to say: the first n data points for the second case are missing. So far so good. However,

Hi, I've looked through the posts but couldn't find a solution to this. I'd be really grateful if someone could help, I'd like to convert a data file of mutual information that is formatted as a matrix:             TF1    TF2    TF3    TF200... Gene1    0.0    0.2    0.2 Gene2    1.4    0.0    2.8 Gene3    0.3    0.6    1.7 Gene6000.... To a list: Gene1    TF1    0.0 Gene1