similar to: error with summary(vector)??

Hi everyone,   I'm trying to do an Logistic Politomic Regression in R. Because I have my resposes variables and the aswer is 0 and 1 in 3 bacterial genes. Somebody know how to do this in R in a easy way?   Thank so much,   José Bustos Facultad de Ciencias Universidad Catolica Ssma. Concepcion Chile --- El lun, 31/8/09, r-help-request@r-project.org <r-help-request@r-project..org>

I have a data set like this I want to assign "outward" to Y if sc <90 and assign "inward" to Y if sc>=90. then cbind(p1982,Y) to get like these p aa as ms cur sc Y 1 154l_aa ARG 152.04 108.83 -0.1020140 92.10410 inward 2 154l_aa THR 15.86 28.32 0.2563560 103.67100 inward 3 154l_aa ASP 65.13 59.16 0.0312137 7.27311 outward 4 154l_aa CYS 57.20 49.85

>>>>> Martin Maechler <maechler at stat.math.ethz.ch> >>>>> on Thu, 19 May 2016 10:26:44 +0200 writes: >>>>> <frederik at ofb.net> >>>>> on Wed, 18 May 2016 15:03:31 -0700 writes: >> Readline <= 6.2 shouldn't require the SIGWINCH patch, so >> if older versions have trouble finding

Hi, I have 6 career types, represented as a factor in R, coded from 1 to 6. I need to use the effect coding (also known as deviation coding) which is normally done by contr.sum, e.g. contrasts(career) <- contr.sum(6) However, this results in the 6th category being the reference, that is being coded as -1: $contrasts [,1] [,2] [,3] [,4] [,5] 1 1 0 0 0 0 2 0 1 0

I had a regression in config.site so the nightly build didn't. Retrying.... Looks like it will build, but the ctl-R, ctl-C bug is still present on OSX (w/Simon's libs). This _was_ fixed for a while, was it not? (The NEWS entry is also wrong: The issue existed before readline 6.3) -pd On 24 May 2016, at 12:55 , Martin Maechler <maechler at stat.math.ethz.ch> wrote: > > Can

I agree with Martin's summary of the situation, and with the updated NEWS entry. I'm not familiar with Subversion, can you tell me the command to use? (I tried "svn co https://svn.r-project.org/R/" but it seems to be downloading all branches) Frederick On Tue, May 24, 2016 at 04:30:11PM +0200, Martin Maechler wrote: > >>>>> peter dalgaard <pdalgd at

When I run the following function HQ2 <- function(n) { nv <- 6 * sqrt(n) fcn <- function(z) { pchisq(z^2 / 36, n - 1) * dnorm(nv - z) } ## I want the integral from 0 to infinity: f.Inf <- integrate(fcn, 0, Inf) ## Doc: "Don't do this": f.100 <- integrate(fcn, 0, 100) cbind(f.Inf, f.100) } I get, for n = 9 and

I have a data set p1982, its structure is the following Then I take 20 observations from this dataset, and assign to pr. in p1982, p has 1982 levels, in dataset pr, p should have 1 levels. But I do str(pr), it shows that p still has 1982 levels. also for these > pr$aa [1] ARG THR ASP CYS TYR ASN VAL ASN ARG ILE ASP THR THR ALA SER CYS LYS THR ALA LYS Levels: ALA ARG ASN ASP CYS GLN

Hello list, how could this happen? The Server doesn't respond from time to time with a high load avarage. We found a suspicious smbd process: top - 13:43:07 up 1 day, 2:27, 5 users, load average: 32.49, 58.41, 37.95 Tasks: 1196 total, 5 running, 1190 sleeping, 0 stopped, 1 zombie Cpu0 : 14.7% us, 3.8% sy, 0.0% ni, 79.8% id, 1.3% wa, 0.0% hi, 0.3% si Cpu1 : 1.3% us, 84.6%

>>>>> Keith O'Hara <keith.ohara at nyu.edu> >>>>> on Tue, 24 May 2016 12:47:43 -0400 writes: > svn checkout https://svn.r-project.org/R/trunk/ <target-directory> yes, indeed. thank you, Keith. and from then on only cd <target-directory> svn up (which is short for 'svn update'). Another hint: Then do *not*

>>>>> <frederik at ofb.net> >>>>> on Tue, 24 May 2016 15:15:17 -0700 writes: > Thank you, Martin. I linked to your message in a comment here so maybe > other people will know about that useful technique: > http://singmann.org/installing-r-devel-on-linux/#comment-161 > However, when I try it, I get an error: > $ make

Hi All, We have a xen server and using 8 core processor. I can see that there is 99% iowait on only core 0. 02:28:49 AM CPU %user %nice %sys %iowait %irq %soft %steal %idle intr/s 02:28:54 AM all 0.00 0.00 0.00 12.65 0.00 0.02 2.24 85.08 1359.88 02:28:54 AM 0 0.00 0.00 0.00 96.21 0.00 0.20 3.19 0.40 847.11 02:28:54 AM