similar to: Multiple missing values

Hi, Is the following function built in somewhere? concat = function(v) { res = "" for (i in 1:length(v)) res = paste(res,v[i],sep="") res } e.g. > concat(c("12","3","45")) [1] "12345" Cheers, Andy -- Andy Fugard, Postgraduate Research Student Psychology (Room F3), The University of Edinburgh, 7 George

Hi there, Am I correct to believe that tetrachoric correlation is a special case of polychoric correlation when there are only two levels to the ordered factor? Thus it should be okay to use hetcor from the polycor package to build a matrix of correlations for binary variables? If this is true, how can one estimate 95% confidence intervals for the correlations? My guess would be mat =

Dear all, I'm struggling to find the best (set of?) function(s) to do repeated measures logistic regression on some data from a psychology experiment. An artificial version of the data I've got is as follows. Firstly, each participant filled in a questionnaire, the result of which is a score. > questionnaire ID Score 1 1 6 2 2 5 3 3 6 4 4 2 ...

I would like to be able to randomise presence-absence (i.e. binary) matrices whilst keeping both the row and column totals constant. Is there a function in R that would allow me to do this? I'm working with vegetation presence-absence matrices based on field observations. The matrices are formatted to have sites as rows and species as columns. The presence of a species on a site is

Dear all, The background is that I'm trying to fix this bug in the geometry package: https://r-forge.r-project.org/tracker/index.php?func=detail&aid=1993&group_id=1149&atid=4552 Boiled down, the problem is that there exists at least one matrix X for which det(X) != 0 and for which solve(X) fails giving the error "system is computationally singular: reciprocal condition

Here's a patch (against 2.2.2) which adds battery.temperature to the ietf-mib snmp driver. While looking at ietfmib.h I noticed that the value for IETF_OID_IN_LINES is actually that for upsOutputNumLines.0. Does anyone know whether that's deliberate or a typo? -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: application/x-patch

Perhaps someone else could confirm this, I have an APC UPS with a built-in network management card with environmental sensor for temp/humidity. The OID index for the sensors appears to be .1 (the last component for each of the OIDs listed below). The original source code has them at .0 drivers/apccmib.h #define APCC_OID_IEM_TEMP ".1.3.6.1.4.1.318.1.1.10.2.3.2.1.4.1" #define

Hello list. I'm contemplating buying new batteries and putting an older APC Smart-UPS 3000XL back into service. Would that work well with NUT? Any drawbacks to consider? Thank you, Tobias -- Tobias Balle-Petersen Kontrapunkt Group / System Administrator

Saludos, tengo una base de datos con tres variables: especie, tratamiento y abundancia. Quiero hacer comparaciones múltiples con Tukey. Ya había usado el comando TukeyHSD(aov(x~y)). La cuestión esque ahora tengo varios niveles del factor especie y varios niveles del factor tratamiento. Si hago TukeyHSD(aov(abundancia~especie:tratamiento)) me salen todas las comparaciones posibles, pero yo sólo

We recently had to replace an old APC SmartUPS 3000 with a newer 3000VA (the 2U one). It has both serial and USB in the back (with an RJ-45 on the end of the serial cable). I can't get the apcsmart driver to connect. I figured it should work the exact same as the old UPS, with /dev/ttyS0 as the port (I'm using the same port on the server that the last one was plugged into). I even

Hello all, There's a new mailing list for researchers in psychology who are learning and using R. New users of R are especially welcome. To join, venture to http://www.jiscmail.ac.uk/archives/psych-r.html and click "Join or leave the list (or change settings)" You can also peruse the archives from that page. Best wishes, Andy -- Andy Fugard, Postgraduate Research

Hi all, Is there an R function to generate random integers? Thanks in advance. Sincerely Anup --------------------------------- Now that's room service! Choose from over 150,000 hotels [[alternative HTML version deleted]]

I maintain the geometry package, which integrates the Qhull C library (http://qhull.org) into R. The Qhull function I hook into requires an open FILE handle as one of its arguments. I had set this file handle to stdout, but now R check NOTEs the presence of stdout, and the CRAN maintainers asked me to get rid of these NOTEs. Including the following defines means the checks are passed on CRAN:

Hi folks, I am new to lme in R, and I have a question regarding to the effect of scale function on the lme. When I use the function to scale and centre the levels of the fixed effects (e.g., X and Y; both have two levels) and write them to new columns: ex: dat$cX<-scale(as.numeric(dat$X),center = TRUE, scale = FALSE) dat$cY<-scale(as.numeric(dat$Y),center = TRUE, scale = FALSE) and compare

Hi, I'm trying to set up a prediction software, now i testing the performance of my method, so i need to calculate a ROC curve, specially auc, cut-off, sens and spec, i just looking at ROCH package, but it's a mass for me, i'm not a math guy and I'm getting lost Could any of you recommend me an easy-to-use package to do this task? i just have a list of positive/negative samples

Hello, I implemented the Beta binomial distribution following the patterns of the binomial distribution code and inspired by JAGS' code [1]. I have studied the code carefully but it's my first run in the R internals. Can somebody review the code and if everything it's ok commit to the repository? [1]

Dear all I have Kubuntu linux and have updated to the latest version (6.06 dapper). I do not know why but now I can not make no plots. For instance, when I type >hist(...) this is the message I get: >can't find X11 font > Error in X11 (display, width, heigth......) > unable to start devide X11 > In addition: Warning messages: > 1:locale not supported by Xlib: some X ops

Dear all, I'm working on some data from an experiment on the breeding behavior of birds. In short, I have been measuring how the time spent on performing a certain task (variable 'mean_on_active') differs over time (variable 'day', 2 levels) across three experimental categories (variable 'treat'; levels 'C', 'R', 'E'). The model shows a

Hi, I need Beta binomial and Beta negative binomial functions but in R there is only SuppDists package which provide this distributions using a limited parameter space of the generalized hypergeometric distribution (dghyper & Co.) which provide a limited parameter space for Beta binomial and Beta negative binomial functions (e.g. alpha + beta <1 in the Beta negative binomial). I've

Hi, My name is Laura. I'm a PhD student at Queen's University Belfast and have just started learning R. I was wondering if somebody could help me to see where I am going wrong in my code for estimating the parameters [mu1, mu2, lambda1] of a 2-phase Coxian Distribution. cox2.lik<-function(theta, y){ mu1<-theta[1] mu2<-theta[2] lambda1<-theta[3]