similar to: FWER + multiple linear models

Hello, I urgently need a histogram of the vector geo1$ds1 but: > hist(geo1$ds1) Error in hist(geo1$ds1) : unused argument(s) ( ...) I can use density, summary, boxplot and stem with this vector but not hist ? Any suggestions ? Thanks in advance, Detlef Groth > density(geo1$ds1) Call: density(x = geo1$ds1) Data: geo1$ds1 (1412 obs.); Bandwidth 'bw' = 0.2805

The function heatmap() allows to specify row/column labels via the options labRow/labCol. From the code of heatmap(), I understand that when no labels are specified, the row/column labels (or indices) of the input matrix are taken as labels and re-ordered together with the rows and columns of the matrix before plotting, whereas labels supplied via labRow/labCol are plotted in the original order.

Dear all, I wonder, if anybody experiences similar problems and if there are any simple solutions to be suggested. We observe that R causes a lot of network traffic and thus slows down the performance of the whole network. When tracing the network traffic on the machine which serves the R installation via NFS, we see thousands of requests at initialization of R processes and regular calls,

Dear R-list members, I would like to pose a question about the use and results of the glm() function for logistic regression calculations. The question is based on an example provided on p. 229 in P. Dalgaard, Introductory Statistics with R, 2nd. edition, Springer, 2008. By means of this example, I was trying to practice the different ways of entering data in glm(). In his book, Dalgaard

You asked the same question on the Bioconductor mailing list back in August. At that time, you suggested yourself a solution for how the adjusted p-values should be interpreted. I answered your query and told you that your interpretation was correct. So I'm not sure what more can be said, except that you should read the article Wright (1992), which is cited in the help entry for p.adjust(),

Is it any way? This code doesn''t work: %%%%%%%%%%%%%%% \multido{}{10}{ <<fig=T>>= ... @ } %%%%%%%%%%%%%%% Or: %%%%%%%%%%%%%%% <<1, fig=T>>= ... @ <<>>= <<1>> @ %%%%%%%%%%%%%%% Serguei Baranov, MSc. Max-Planck-Institut fur Molekulare Genetik Abteilung Lehrach, room 2.112 Ihnestrasse 73 14195 Berlin-Dahlem Germany Tel. 49-30-8413-1667 Fax.

Dear useRs, last week I asked you about a problem related to isoMDS. It turned out that in my case isoMDS was trapped. Nonetheless, I still have some problems with other data sets. Therefore I would like to know if anyone here has experience with how well isoMDS performs in comparison to other non-metric MDS routines, like Minissa. I have the feeling that for large data sets with a high

Hello, after the line cat("%",as.character(as.name(match.call())), "\n%\n", file=fi, append=T) R complains about Error in as.name(x) : character argument required Is it sufficient to delete as.name()? After doing that, R complains Error: comparison is possible only for vector types Don't know which function call. Can you give me a hint how to debug it? There

Hello all, I used the following script to generate a plot, but when I try to place axes, the axes start at the end of my plot and continue off the plot; the attached image shows what happens when I try to generate axes. library(fields) as.matrix(read.table("Matrix.txt", sep="\t"))->x #the matrix is a square matrix measuring 104x104; it was too large to attach #on left:

Josh, Those speed improvements are great. Quick test on an 365MB wav of Neil Young's eponymous 1st album on a PIII 650E (Coppermine) under FreeBSD 4.2: User Sys MB shorten23 0:46s 4.0 208.3 shorten31 1:24s 5.3 208.6 flac CVS 4:21s 7.5 199.2 It's getting near to be a useful replacement for shorten. Only 3x slower. Any idea why it uses so much more sys time? I would've thought

--- Mark Powell <M.S.Powell@salford.ac.uk> wrote: > Josh, > Those speed improvements are great. Quick test on > an 365MB wav of Neil > Young's eponymous 1st album on a PIII 650E > (Coppermine) under FreeBSD 4.2: > > User Sys MB > > shorten23 0:46s 4.0 208.3 > shorten31 1:24s 5.3 208.6 > flac CVS 4:21s 7.5 199.2 > > It's getting near to be

That article is a bit too dismissive. I agree that one cannot hear the difference between 48KHz/16bit and 192KHz/24bit if you just transfer the data directly to the audio output device. As such, there is no good reason for Opus to support higher than 48KHz (especially since this is lossy compression, anyway). However, in general, that's not all you do with audio data. 192KHz is useful for

Hi, I have just seen Stefan Steinhaus' web page : http://www.uni-franfurt.de/~stst/ncrunch.html I think it would be nice to include "R" as well. I have taken Forrest Young's email on stat-lisp list and changed the stuff for "R" :) Here it is: (someone please check this so we can also send it to Stefan Steinhaus.

Hi all! I do not know openssh patch policy so I am just sending the patch to the mailing list. Sorry for inconvenience. Ssh-agent seems to be too slow if you need to access thousands of servers. This is a simple patch to enable threading in ssh2 authentication. Patch adds "-p numthreads" option and defaults to the number of processors. I've tested it as I could, but

Hi all , I have a dataframe (Hypertension) with following headers :- > Hypertension ID Hypertension(before drug A) Hypertension(On drug A) On drug B? Healthy diet? 1 160 90 True True 2 190

library(glmpath) data(heart.data) # heart.data is a list, $y a vector, $x a matrix data <- data.frame(x=I(heart.data$x), y = heart.data$y) > data[1:2,] x.1 x.2 x.3 x.4 x.5 x.6 x.7 x.8 x.9 y 1 160 12 5.73 23.11 1 49 25.3 97.2 52 1 2 144 0.01 4.41 28.61 0 55 28.87 2.06 63 1 > dimnames(heart.data$x)[[2]] [1] "sbp"

Hi, I am trying a linear regression model where the dependent variable is the size of the heart corrected for the patient's height and weight. This is labelled as LAVI. The independent variables are race (european or non-eurpoean), age, sex (male or female) of the patient and whether they have diabetes and high blood pressure. sample size 2000 patients selected from a community. when I

I am trying to figure out how to calculate p-values for the difference in prevalence of a risk factor between men and women. For example, I find that 277 out of 710 male patients and 125 out of 305 female patients have obesity, what is the p-value for their difference? If there is a package that can calculate this in bulk, I would appreciate to learn about it! Thank you [[alternative HTML

Hi everyone, I've only been using R for a week or so now, and am now required to work with some pretty mixed data at my hospital. Basically, we're looking to see if hypertension status (nominal) affects cardiac function (multiple & continuous). In order to do this, we've collected multiple echocardiograph data reports on each patient. 1. Each patient is coded either as 0 or 1

OpenSSH Developers, I'm researching whether or not OpenSSH is a viable commercial alternative to F-SECURE SSH or SSH.COM's ssh, but I'm not getting the kind of results that I expected from a "Non patent encumbered ssh client". When I attempt to build OpenSSH against an OpenSSL build without rc5, idea, or rsa it bombs since OpenSSL doesn't place the header files in the