similar to: Saving R graphics with version 2.3.1

Hi, I installed the package Hmisc with the command install.packages("Hmisc") without errors. When I try to load the library with command library(Hmisc) I get the error > library(Hmisc) Error in library(Hmisc) : there is no package called 'Hmisc' > version _ platform i386-pc-linux-gnu arch i386 os linux-gnu system i386,

Dear all, I tried to use "step" function to do model selection, but I got an error massage. What I don't understand is that data as data.frame worked well for my other programs, how come I cannot make it run this time. Could you please tell me how I can fix it? ***************************************************************************************************

Hello, I was wondering if I can plot two curves I get from "density(data)" into one plot. I want to compare both. With the following commad, I just get one curve plotted: plot( density(mydata) ) Sorry for this stupid question but I could not find a solution until now... Antje

Dear all, I want to design a histogram and I need to have the frequency at certain points. For example I have the following 2 columns: *X Y* 0.1 25 0.4 22 0.45 11 0.55 21 I want the chart to have 4 columns. First column is from 0.0-0.1 (on X) and frequency is 25. Next colum is wider and form 0.1-0.4 with 22 frequency. Next column is narrow with 11 frequency and the last column is

Hi ! I am trying to get at the covariance of the predictions of a linear model. Suppose the we have: > x<-runif(1000) > y<-2 + 25x*x +rnorm(1000) > lm1 <-lm(y~x, data = data.frame(y = y, x=x)) > x.pred <-runif(10) > y.hat <- predict(lm1, newdata = data.frame(x=x.pred)) I was wondering how to get an estimate of the covariance of y.hat which would be a 10 x 10

I was reading about the project drawer feature in Textmate, which is Mac only. Is there a similar feature in a Windows based text editor that works with R. This feature sounds really useful. Thanks, Graham [[alternative HTML version deleted]]

Hello! I need the Code of the functions runif and rexp. Where can I get them? Can You help me? I thank you for an answer! bye Ilka --

Hi, Suppose I have a multivariate response Y (n x k) obtained at a set of predictors X (n x p). I would like to perform a linear regression taking into consideration the covariance structure of Y within each unit - this would be represented by a specified matrix V (k x k), assumed to be the same across units. How do I use "lm" to do this? One approach that I was thinking of

Hi R users I have a dataset which represents points that are market by patients as the source of pain. Basically the patients indicates by a cross on a chest pictures where he/she thinks is the source of pain. The data was then digitalized by divinding the chest into small squares and each square was give value 1 if it was the center 2 if it was touched by the markings and 3 if it was not

Hi, Sorry for the question, I know it should be basic knowledge but I'm struggling for two hours now. How do I select only the first entry of each list member and ignore the rest? So for > $"121_at" > -113691170 > $"1255_g_at" > 42231151 > $"1316_at" > 35472685 35472588 > $"1320_at" > -88003869

Hi All, I need some help in how one can implement maximumlikelihood estimation for models with discrete hidden variables in EM in R. Regards [[alternative HTML version deleted]]

Hello All, I have 8 variables named a b c d e f g h I need to create four variables from these 8 vraibles in R. the new variables are ab,cd,ef,gh. Can anyone pleas help me thanks, Pratap --------------------------------- [[alternative HTML version deleted]]

hie when i use plot.survfit to plot more than one graph why I only see the last graph how do i see the other graphs.for example n=20 n1=n/2 n2=n/4 a11=4;a12=4 ;a21=4 ;a22=4 t1<-array(1,c(n1)) t2<-array(2,c(n1)) treatgrp=matrix(c(t1,t2))

Hi, I''m trying to configure QoS on my linux in the following manner: I have a main link with 64K, so I divided it in 3 classes of 18K, 14K and 9K with an excess (not used for classified traffic, just to be shared) of 23K. This excess should be distributed proportonally among the 3 classes, that is, the class that has more rate should borrow more bandwidth. What is happening is

Hello, I'm Vitor Afonso Strabello and I need to be added as a member into the Wiki to edit/create some tips and entries on the Wiki. My userid is VitorStrabello. I would like to contribute on network, configuration and installation (to begin, I'll post something in kickstart section on the wiki...) tips/tricks/how-tos, so on at the CentOS wiki. Also I'm a member of the CentOSBR

I would like to buy a basic statistics book (experimental design, sampling, ANOVA, regression, etc.) with examples in R. Or download it in PDF or html format. I went to the CRAN contributed documentation, but there were only R textbooks, that is, textbooks where R is the focus, not the statistics. And I would like to find the opposite. Other text I am trying to find is multivariate data

hi, thanks Dr. Murdoch i'd appreciate if anyone on r-help could help me narrow this down? i believe the segfault occurs because there's a single line with 4GB and also embedded nuls, but i am not sure how to artificially construct that? the lodown package can be removed from my example.. it is just for file download cacheing, so `lodown::cachaca` can be replaced with `download.file`

Hey folks, I'm from Brazil and I have the following doubt. May I use an asterisk box with some cards to act as an PSTN simutator between a little amount of sites? I will start to think about it in a voice lab for my studies. Kind regards, Vitor -- Vitor Afonso Strabello MSN: vstrabello at hotmail.com

hi, i realized that the segfault happens on the text file in a new R session. so, creating the segfault-generating text file requires a contributed package, but prompting the actual segfault does not -- pretty sure that means this is a base R bug? submitted here: https://bugs.r-project.org/bugzilla3/show_bug.cgi?id=17311 hopefully i am not doing something remarkably stupid. the text file

Hi everyone, I'm using the function 'HWE.exact' of 'genetics' package to compute p-values of the HWE test. My data set consists of ~600 subjects (cases and controls) typed at ~ 10K SNP markers; the test is applied separately to cases and controls. The genotypes are stored in a list of 'genotype' objects, all.geno, and p-values are calculated inside the loop over all