similar to: Reading ASCII files

Dear all, I have a surprising problem with the representation of frequencies in a histogram. Consider, for example, the R code: b<-rnorm(2000,3.5,0.3) hist(b,freq=F) When I plotted the histogram, I expected that values in the y-axis (the probability) varied between 0 and 1. Instead, they varied within the range 0-1.3. Have you got any suggestion for obtaining a correct graph with

Dear all, is there a automatic method to obtain Odds Ratio estimates and their confidence intervals from a GLM Logistic model? Thanks Cristian ============================================= Cristian Pattaro ============================================= Unit of Epidemiology & Medical Statistics Department of Medicines and Public Health University of Verona cristian@biometria.univr.it

Dear all, I'm using the package "Genetics", and I'm interested in the computation of D' statistics for Linkage Disequilibrium, for which the LD() command has been realised. Unfortunately I don't find any reference on "how" the D' is computed by the LD() function. In the package documentation it is generally referred as "MLE" estimation, but

Hi! How can I fill with colors a portion of a graph (e.g.: I want fill in red the area within two confidence intervals)? Thank you very much Cristian ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Cristian Pattaro ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Unit of Epidemiology & Medical Statistics University of Verona Tel +39 045 8027668 fax +39 045 505357 cristian at biometria.univr.it

Hi! After installing R 2.0.0, I have problem with the "help.search" as below: > help.search("table") Error in help.search("table") : could not find package 'file23525' > help.search("mean") Error in help.search("mean") : could not find package 'file23525' I have had the same problem with other request... What is

Hi people, I have a text file like this one posted: snp_id gene chromosome distance_from_gene_center position pop1 pop2 pop3 pop4 pop5 pop6 pop7 rs2129081 RAPT2 3 -129993 "upstream" 0.439009 1.169210 NA 0.233020 0.093042 NA -0.902596 rs1202698 RAPT2 3 -128695 "upstream" NA

Dear Sir/Madam, I found your email address and your correspondence with R-users. I hope you could help me with this question about the function "ur.ers" in the package of "urca". It is an improved unit root test (Elliott et al. 1996 Econometrica). Do you know how to extract the value of the test statistic from the output? The only thing I can get is the print-out of all

I'm new to R (previously used SAS primarily) and I have a genetics data frame consisting of genotypes for each of 300+ subjects (ID1, ID2, ID3, ...) at 3000+ genetic locations (SNP1, SNP2, SNP3...). A small subset of the data is shown below: SNP_ID SNP1 SNP2 SNP3 SNP4 Maj_Allele C G C A Min_Allele T A T G ID1 CC GG CT AA ID2 CC GG CC AA ID3 CC GG nc AA

I have a big data frame and I have extracted a bit by doing: > y<-d[1:10,6] > y [1] Headings 0 -49 -98 -49 -41 -120 -155 -204 -169 92329 Levels: -0 -1 -10 -100 -1000 -10000 -10001 -10002 -10003 -10004 -10005 -10006 -10007 -10008 -10009 -1001 -10010 -10011 -10012 -10013 -10014 -10015 -10016 -10017 -10018 -10019 -1002 -10020 -10021 -10022 -10023

Hi Cristian, I don't know about a R routine for exact CIs, but I found a function called "ageadjust" using a gamma distribution approximation some time ago -take a look at http://medepi.org/epitools/rfunctions/index.html. That Webpage seems a bit outdated now, but the calculations given in the ageadjust functions should still work, I think. Maybe this helps; regards Heinrich.

Hello,     Looking up a users group membership I'm showing different results on each DC. UID and GID mapping appears consistent but not all group membership is displayed. I've verified idmap.ldb is backup up and copied over to the other DC's. I do notice when taking a hot backup of idmap.ldb, the file size is dramatically smaller than the original. Using Microsoft RSAT to view

Dear R helpers, I was doing a genetic project with two datasets X and Y. There are some IDs in both data sets, and others in either data set. I used "merge(x,y,by="ID",all=TRUE)". The data set Y contains a variable (a genotype) which is also in data X. When I merge X with Y, these two variables were automatically re-named by appending .x and .y to the original variable names.

After some more screwing around with leaving and rejoining the ADS domain I was finally able to access a share with "valid users =" set to a domain group I was a member of. The _only_ change I made after this was to add yet another group to the valid users on the share and restart samba...after that I could no longer access the share. I removed the additional group, restarted samba and

Hi all! Do you know if there is any R function/package that can be used to estimate "logit" models with panel data and forecasting? Thanks, Moysés. -- Moysés Nascimento Bacharel em Estatística/UFES Mestre em Estatística Aplicada e Biometria/UFV Doutorando em Estatística e Experimentação Agropecuária/UFLA moysesnascim@gmail.com [[alternative HTML version deleted]]

Dear all, I'm working with R 1.2.2 on win98 (256Mb RAM) with ESS 5.1.18. It happens that the start-up is very slow. I don't know exactly how to measure, but with a simple watch I recorded 15 seconds, while the startup of S+ 2000 (Sqpe+4 under ESS) takes around 3 seconds. Tha fact is that sometimes it takes almost a minute (58.30 sec recorded) for the prompt to appear! And ESS gives me an

Dear R users, I would like to employ count data as covariates while fitting a logistic regression model. My question is: do I violate any assumption of the logistic (and, more in general, of the generalized linear) models by employing count, non-negative integer variables as independent variables? I found a lot of references in the literature regarding hot to use count data as outcome,

Hi, I want to ask you if asterisk, when I use the command park(), gives me for example a variable that contains the slot position where it parks the call or if it only tells me (audio) in the channel this position number? In other words, is there a way to obtain and use the value of the slot position when the call is parked? Thanks.

Hi all, how can I catch the event generated when a parked call is hung up? In my dialplan, when arrives a call to a specific number, Asterisk parks the call and announces the parking slot to a number. But if the user hangs up the parked call, I don't know how to catch the event, from dialplan, that removes the call from the parking slot. I want to know if there is a method for do this.

Hi all, I have a problem to park the call with #700 when the SIP phone and the asterisk PBAX run on the same machine. Call parking works well if the SIP phone that I use is on a remote pc in comparison with the asterisk on the which it is registered. If the sipphone is in the same pc where asterisk is installed the command #700 doesn't work, and in the asterisk's CLI I can see no

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20071208/5409f1a7/attachment.pl