Re: [R] How I can read the binary file with different type?

On Fri, 22 Aug 2008, Bingxiang Miao wrote: Hi all, I have a binary file which have 8*100 bytes. The structure of the file is as follows: every eigth bytes represent 2 data:the first four bytes is the little-endian for integer, the next four bytes is the little-endian for floating32. The

Re: [R] Help Regarding 'integrate'

As my post showed, it is a scaling issue. The function has so small a peak that it is effectively 0 -- when scaled sensibly integrate then works out of the box. On Thu, 21 Aug 2008, Thomas Lumley wrote: On Thu, 21 Aug 2008, Moshe Olshansky wrote: The phenomenon is most likely caused by

Re: [R] R memory limitations under Windows vs UNIX

On Fri, 22 Aug 2008, Fiona Giannini wrote: Hi all, Our section at work is looking at buying a high-powered computer (32 Gb RAM, 64-bit) to be able run models on large datasets and have more processing power for the software we commonly use. We mostly use R to fit GLMs and GAMs. Our

[R] Predictive Analytics event Sept 24-25, Chicago

Hi, I wanted to make sure you were all aware of these upcoming events. There is a seminar in Predictive Analytics on Sept 24-25 in Chicago, and two following in D.C. (Oct.) and San Francisco (Nov.). This is intensive training for managers, marketers, and IT people who need to make sense of

[R] subset grouped data with quantile and NA's

I can't quite seem to solve a problem subsetting a data frame. Here's a reproducible example. Given a data frame: dat - data.frame(fac = rep(c(a, b), each = 100), value = c(rnorm(130), rep(NA, 70)), other = rnorm(200)) What I want is a new data frame (with

[R] Parabolic cylinder function

Dear all, I need your advice since I am looking for an implementation of the parabolic cylinder function in R. I found implemantations of the hypergemetric functions (the Whittaker and the confluent hypogeometric functions) in the package fAsianOptions but the parabolic cylinder function was

Re: [R] GAM-binomial logit link

Dear all, I'm using a binomial distribution with a logit link function to fit a GAM model. I have 2 questions about it. First i am not sure if i've chosen the most adequate distribution. I don't have presence/absence data (0/1) but I do have a rate which values vary between 0 and 1. This

[R] Censored Poisson Data

Dear list, I am wondering whether R allows to perform a Poisson regression on counting data which include censored observations, that is, observations of the form 2 events or more or less than 2 events or even 1 or 2 events. Can anyone give me a hint whether this is already possible and how

Re: [R] Simple look up

Please take a look at the manuals before posting (DO read the posting guide!) Bart PS: yourdataframe[yourdataframe$TAZ==103,2] PDXRugger wrote: So i am hoping this solution is simple, which i beleive it is. I would like to look up a value in one column of a data set and display the

Re: [R] read.csv : double quoted numbers

On Thu, Aug 21, 2008 at 12:54 AM, Gabor Grothendieck [EMAIL PROTECTED] wrote: The problem is the commas, not the quotes: First thanks to everyone for taking out time to help me. As Mr. Gabor rightly corrected commas was my problem and not the double quotes (in numerical values). This solves

Re: [R] skipping values when reading binary files

Pieter Beck wrote: Hi, I would like to skip a preset number of values when reading in a binary file (readBin). Is there any way to do this? See ?seek Uwe Ligges kind regards, and thanks in advance for the help. Pieter Beck [[alternative HTML version deleted]]

Re: [R] barplot with anchored bars

Paulo Barata wrote: Dear R list members, How to produce barplots anchored to the x-axis (not floating above the x-axis) with a box around? With both following codes, the lower horizontal line of the box is below the y = 0 line: # first code x - c(1,2,3,4) barplot(x,yaxs='i') box() #

Re: [R] read.csv : double quoted numbers

gsub fails on this file due to a particular value C￴te d'Ivoire which is present in the file. This is when my locale is set to en_US.utf8 but when I set locale to C it (gsub) works. Have a look at ?file, especially the section on Encoding. cu Philipp -- Dr. Philipp Pagel

[R] chisq

Hi, the chisq value for .05 level of signif. and 3 DF is 7.81 (one-sided). Is there a way to calculated it in R (giving level of signif. and DF)? Marco __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] Help on competing risk package cmprsk with time dependent covariate

Dear R users, I d like to assess the effect of treatment covariate on a disease relapse risk with the package cmprsk. However, the effect of this covariate on survival is time-dependent (assessed with cox.zph): no significant effect during the first year of follow-up, then after 1 year a

Re: [R] GAM-binomial logit link

Dear all, I'm using a binomial distribution with a logit link function to fit a GAM model. I have 2 questions about it. First i am not sure if i've chosen the most adequate distribution. I don't have presence/absence data (0/1) but I do have a rate which values vary between 0 and 1. This

[R] Retrieving sample points

Hi everyone. My problem is about multivariate analysis. Given a multivariate data X, of dimension (n by p), n is the number of observations while p is the number of variables. I wanted a sample of size p+1 that will have its product of the eigen values of the covariance matrix to be the minimum

Re: [R] psychometric functions

Paul, thanks a lot for your advice! I got that kuss et al. paper and started to read, this is more than I had expected! Very cool paper, and an algorithm implemented in R that seems to be even superior to psignifit in matlab. So thanks again! mario

[R] swap

Hello everybody,   I wonder if there is any swap function in R that does the following: x - seq(1,10,12) x1 - swap(x) x1 1 8 3 4 5 6 7 2 10 Thank you very much!   Amor __ Schutz gegen Massenmails. [[alternative HTML version deleted]]

Re: [R] Very confused with class

Hi Rich, Richard M. Heiberger wrote: Dan, The real problem is the use of csv files. csv files don't handle missing values (#VALUE is most likely from Excel), dates, or other complications very well. Read your Excel file directly into R with one of the packages designed

Re: [R] Help on competing risk package cmprsk with time dependent covariate

Hello, Something i don't understand in your question. Is treatment a time-dependent covariate? That is, do patients receive the treatment at the beginning of the study or later? cmprsk cannot handle time-dependent covariates. But if treatment is a baseline covariate, but has a time-varying

Re: [R] swap

I wonder if there is any swap function in R that does the following: x - seq(1,10,12) This just makes x equal to 1. I'm guessing you meant something like x - 1:10 x1 - swap(x) x1 1 8 3 4 5 6 7 2 10 It's not terribly clear what you want swap to do. You can reorder the elements of x

[R] Coordinate systems for geostatistics in R

Hi, I read somewhere that when carrying out geostatistical analysis in R you should not use latitude and longitude...can anyone expand on this a little for me, and what would be the best coordinate system to use? I have my data in a geographic coordinate system, WGS84, decimal degreesis

Re: [R] Large data sets with R (binding to hadoop available?)

Hi Apart from database interfaces such as sqldf which Gabor has mentioned, there are also packages specifically for handling large data: see the ff package, for instance. I am currently playing with parallelizing R computations via Hadoop. I haven't looked at PIG yet though. Rory

Re: [R] call perl

Maybe ?system Might be what you need here? -Rory -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jay an Sent: 22 August 2008 03:08 To: r-help@r-project.org Subject: [R] call perl Hi, It may be the old question. can anyone tell me how to call perl in

Re: [R] help needed for HWE.exact in library genetics

You could follow the advice given when loading the library (see the code you posted for details) and use the enhanced genetics packages to cross-validate the results (ideally you should get the same answer). The results not weird though. Your working with SNPs and having a homozygote with a

[R] Bias reduction method

Dear all, I'm doing a simulation study for an MM1 queue and I would like to know if there are bias correction methods in R for time series (estimations of warm-up points or something else). Hope you can help me.

Re: [R] barplot with anchored bars

On Thu, 2008-08-21 at 15:47 -0300, Paulo Barata wrote: Dear R list members, How to produce barplots anchored to the x-axis (not floating above the x-axis) with a box around? With both following codes, the lower horizontal line of the box is below the y = 0 line: # first code x -

[R] simple generation of artificial data with defined features

Dear R-colleagues, I am quite a newbie to R fighting my stupidity to solve a probably quite simple problem of generating artificial data with defined features. I am conducting a study of inter-observer-agreement in child-bronchoscopy. One of the most important measures is Kappa according to

[R] random effects in lmer

Dear all I am running several generalized mixed model using lmer. The models typical look like this: model2xx-lmer(numbers~Treatment+b+(1|Genotype)+(1|Field)+(1|Genotype:Treatment), family=quasipoisson) All factors are categorical. And the output looks like this: Generalized linear mixed

Re: [R] subset grouped data with quantile and NA's

This will also remove the NAs from the output; you will have to change it to also keep the NAs. Wasn't sure what you wanted to do with them. dat - data.frame(fac = rep(c(a, b), each = 100), value = c(rnorm(130), rep(NA, 70)), other = rnorm(200)) # split the

Re: [R] Combining multiple datasets

Here is the way that I usually do it using lapply: myData - lapply(list.files(pattern=.RWL$), function(.file) read.rel(.file, header=TRUE)) # combine into a single dataframe myDF - do.call(rbind, myData) On Thu, Aug 21, 2008 at 10:42 PM, Alison Macalady [EMAIL PROTECTED] wrote: Hi, I've

[R] nesting environments

Hello, I got the hint to use nesting environments (for a GUI problem), but I'm not sure whether I understand what he means (and I did not really succeed with reading documentation - probably I looked for the wrong keywords?). That's what he wrote: Another good way to control the scoping is

[R] swap

Hello Richie, I would like to do three (or k) swap steps in each step just 2 ID recursive swaping x - 1:10 swap - function(x){   a - sample(x,2)   x[x==a[1]] - swap[2]   x[x==a[2]] - swap[1]   return(x)   }   swap(swap(swap(x))) - mix   Is this possible? Thanks you in advance! Amor

[R] Re : Help on competing risk package cmprsk with time dependent covariate

Hello again, I m trying to use timereg package as suggested (R2.7.1 on XP Pro). here is my script based on the example from timereg for a fine gray model in which relt = time to event, rels = status 0/1/2 2=competing, 1=event of interest, 0=censored random = covariate I want to test  

Re: [R] swap

Is this what you want: swap - function(z){ + a - sample(length(z), 2) + z[a] - z[rev(a)] + z + } swap(1:10) [1] 2 1 3 4 5 6 7 8 9 10 swap(1:10) [1] 5 2 3 4 1 6 7 8 9 10 swap(1:10) [1] 8 2 3 4 5 6 7 1 9 10 swap(1:10) [1] 2 1 3 4 5 6 7 8 9

Re: [R] nesting environments

Look at the example section of ?new.env Also the proto package on CRAN defines proto objects (which are just environments with slightly different semantics) and in the Applications section of the home page has links to examples of using proto objects to creates GUIs with gWidgets:

Re: [R] nesting environments

Thanks a lot! I'll give it a closer look :-) Ciao, Antje Gabor Grothendieck schrieb: Look at the example section of ?new.env Also the proto package on CRAN defines proto objects (which are just environments with slightly different semantics) and in the Applications section of the home page

Re: [R] Coordinate systems for geostatistics in R

imicola schreef: Hi, I read somewhere that when carrying out geostatistical analysis in R you should not use latitude and longitude...can anyone expand on this a little for me, and what would be the best coordinate system to use? I have my data in a geographic coordinate system, WGS84, decimal

Re: [R] swap

Hello Richie, I would like to do three (or k) swap steps in each step just 2 ID recursive swaping x - 1:10 swap - function(x){ a - sample(x,2) x[x==a[1]] - swap[2] x[x==a[2]] - swap[1] return(x) } swap(swap(swap(x))) - mix I tried my best with a response before, but if you

[R] How to read malformed csv files with read.table?

Hi, how do I read files that have two header fields less than they have columns? The easiest solution would be to insert one or two additional header fields, but I have a lot of files and that would be quite a lot of awful work. Any ideas on how to solve that problem? ### R stuff:

[R] Re : Help on competing risk package cmprsk with time dependent covariate

Hello again, I m trying to use timereg package as you suggested (R2.7.1 on XP Pro). here is my script based on the example from timereg for a fine gray model in which relt = time to event, rels = status 0/1/2 2=competing, 1=event of interest, 0=censored random = covariate I want to test

[R] draw a function over a histogram

Sorry, I have some troubles with the graph device. How can I draw a function over a histogram? Thank's so much. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] mapping a vector into predefined bins

Hi, I have a mapping M consisting of a vector of bin centers (say 0.01, 0.02, 0.03 ..., 1) with associated values (say red, orange, yellow, ...), stored as a 2-column matrix with column 1 the bin center and column 2 the value. I also have a vector D of data (say 0.9056, 0.4118,

[R] filtering out data

Greetings, Apologies for such a simple question, but I have been trying to figure this out for a few days now (I'm quite new to R), have read manuals, list posts, etc., and not finding precisely what I need. I am accustomed to working in Stata, but I am currently working through the

[R] filtering out data

Greetings, Apologies for such a simple question, but I have been trying to figure this out for a few days now (I'm quite new to R), have read manuals, list posts, etc., and not finding precisely what I need. I am accustomed to working in Stata, but I am currently working through the

Re: [R] draw a function over a histogram

Daniela Garavaglia wrote: Sorry, I have some troubles with the graph device. How can I draw a function over a histogram? Thank's so much. Daniela, What function? Here is one example using density() and using dnorm() x - rnorm(1000,2,2) hist(x,prob=TRUE) lines(density(x,na.rm=TRUE),col=red)

Re: [R] How to read malformed csv files with read.table?

Try this. It will read the file and see if there is a difference and add in the extra headers: x -time/ms C550.KMSCyt_b559.KMSCyt_b563.KMS Cyt_f_.KMS P515FR.KMS Scatt.KMS Zea2.KMSPC P700 0 Point1 -599.500 0.000 0.000

Re: [R] How to read malformed csv files with read.table?

on 08/22/2008 10:19 AM Martin Ballaschk said the following: how do I read files that have two header fields less than they have columns? The easiest solution would be to insert one or two additional header fields, but I have a lot of files and that would be quite a lot of awful work. Any

Re: [R] filtering out data

Exactly what do you want this statement to do: if(merged.dataset$model1 == 1) model1.coeff - lm(dv ~ iv1 + iv2 + iv3, data = merged.dataset) The result in the if is a vector of logical values (assuming that merged.dataset is longer than 1); that is the source of the error message. Do you want

[R] testing for differences in sums - how to?

If I want to test for a salesman effect on average sale, I model it as Sale ~ Salesman. But what if I want to test for a salesman effect on total sales? How do I model that? Given a significant salesman effect, how do I follow up with the equivalent of Tukey's HSD to determine who the

Re: [R] How to read malformed csv files with read.table?

On Fri, 22 Aug 2008, Daniel Folkinshteyn wrote: on 08/22/2008 10:19 AM Martin Ballaschk said the following: how do I read files that have two header fields less than they have columns? The easiest solution would be to insert one or two additional header fields, but I have a lot of files and

[R] Test of Homogeneity of Variances

I am testing the homogeneity of variances via bartlett.test and fligner.test. Using the following example, how should I interpret the p-value in order to accept or reject the null hypothesis ? set.seed(5) x - rnorm(20) bartlett.test(x, rep(1:5, each=4)) Bartlett test of homogeneity

Re: [R] HELP: how to add weight to a [x,y] coordinate

There are several options. You can use the sunflowerplot function. You can plot solid points with a small alpha value (not all graphics devices support this), so that just a few points show up pale, but a lot show up as more opaque. You can use the symbols function (or bubbleplot, or

Re: [R] Test of Homogeneity of Variances

What are your hypotheses? Once you state what they are, interpretation should be straightforward. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Daren Tan Sent: Friday, August 22, 2008 11:18 AM To: [EMAIL PROTECTED] Subject: [R] Test of Homogeneity

[R] lme questions re: repeated measures covariance structure

Hello, We are attempting to use nlme to fit a linear mixed model to explain bird abundance as a function of habitat: lme(abundance~habitat-1,data=data,method=ML,random=~1|sampleunit) The data consist of repeated counts of birds in sample units across multiple years, and we have two questions:

Re: [R] Null and Alternate hypothesis for Significance test

The problem with trying to prove that 2 distributions are equal is that there is exactly one way in which they can be equal and an infinite number of ways that they can be different (and its rather a large infinity). The traditional use of equality as the Null hypothesis works, because we can

Re: [R] swap

I'm guessing that the following (untested) does what is wanted: function(x) { pos - sample(length(x), 2, replace=FALSE) x[pos] - x[ pos[2:1] ] x } Patrick Burns [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User)

Re: [R] lme questions re: repeated measures covariance structure

Dear Mark, I would include the repeated measure as the smallest stratum in the random effects specification: random=~1|sampleunit/year Setting up user-defined variance structures should be possible using for example: weights=varPower(form=~habitat) or also try out the available corStruct()

Re: [R] mapping a vector into predefined bins

Have a look at ? cut Best wishes Wolfgang -- Wolfgang Huber EBI/EMBL Cambridge UK http://www.ebi.ac.uk/huber 22/08/2008 14:14 Dr. G. Lawyer scripsit Hi, I have a mapping M consisting of a vector of bin centers

Re: [R] filtering out data

If I understand correctly the result you are trying to achieve, I think what you may be looking for is this: model1.coeff - lm(dv ~ iv1 + iv2 + iv3, data = merged.dataset[merged.dataset$model1 == 1, ]) This will give you a regression using only the data for which model1 == 1. on 08/22/2008

Re: [R] Test of Homogeneity of Variances

I am testing whether the sample variances are equal. When p-value 0.05 (alpha), should accept null hypothesis (sample variances are equal) or reject it ? The two new examples with each having same sample variances also puzzle me. Why are the p-values different ?

Re: [R] filtering out data

Thank you, thank you everyone! This last worked the best for what I need to do. The other options were definitely helpful in that they help me understand how R is structured (something I haven't been able to get my head around yet) and gave me good ideas about how to manipulate data, however

[R] WinBUGS with R

Dear Users, I am new to both of things, so do not blame me too much... I am busy with semiparametric regression and use WinBUGS to sample posteriors. The code to call Winbugs is as follows: data- list(y,X,n,m) #My variables inits.beta - rep(0,K) inits.beta0 - 0 inits

Re: [R] Combining multiple datasets

Here's a function that does what you asked, you may need to adjust the column names of your data frames before using it. If all your data frames are similar (same number of rows, same years) then try do.call('cbind', yourList). #This function takes a list of data frames and merges them into one

Re: [R] Quantile regression with complex survey data

Just as a follow-up on T. Lumley's post, 2 citations that may be useful in reference to application of quantile regression with survey samples are: Bassett and Saleh. 1994. L_1 estimation of the median of a survey population. Nonparametric Statistics 3: 277-283. (L_1 estimation of median

[R] grep, gsub and metacharacters

Hello, I have an expression that I want to substitute for something else. myvector-c(ajkhfkiuwe,Variable(kg/min)) if I use the following code to extract variable(kg/min) and substitute it for va gsub(myvector[grep(var, myvector, ignore=T)], va, myvector) grep identifies the element of

Re: [R] Survey Design / Rake questions

I *think* I'm making progress, but I'm still failing at the same step. My rake call fails with: Error in postStratify.survey.design(design, strata[[i]], population.margins[[i]], : Stratifying variables don't match To my naïve eyes, it seems that my factors are in the wrong order. If so,

Re: [R] WinBUGS with R

artimon wrote: Dear Users, I am new to both of things, so do not blame me too much... I am busy with semiparametric regression and use WinBUGS to sample posteriors. The code to call Winbugs is as follows: data- list(y,X,n,m) #My variables inits.beta - rep(0,K) inits.beta0

Re: [R] How to read malformed csv files with read.table?

Hi folks, thank you for your friendly and immediate help! Am 22.08.2008 um 17:14 schrieb Prof Brian Ripley: Or, better, use header=FALSE, skip=1 and the col.names arg of read.table(). My solution is reading the files without the headers (skip = 1) and seperately reading the headers with

Re: [R] grep, gsub and metacharacters

Try this: myvector-c(ajkhfkiuwe,Variable(kg/min)) gsub( \\(kg/min\\), va, myvector ) Does that come close to what you want? If not, maybe an example of what you want the result to look like, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED]

Re: [R] How to read malformed csv files with read.table?

On Fri, 22 Aug 2008, Martin Ballaschk wrote: Hi folks, thank you for your friendly and immediate help! Am 22.08.2008 um 17:14 schrieb Prof Brian Ripley: Or, better, use header=FALSE, skip=1 and the col.names arg of read.table(). My solution is reading the files without the headers (skip =

Re: [R] simple generation of artificial data with defined features

On the general question on how to create a dataset that matches the frequencies in a table, function as.data.frame can be useful. It takes as argument an object of a class 'table' and returns a data frame of frequencies. Consider for example table 6.1 of Fleiss et al (3rd Ed): birth.weight -

Re: [R] How to read malformed csv files with read.table?

Hi Brian, Am 22.08.2008 um 19:22 schrieb Prof Brian Ripley: On Fri, 22 Aug 2008, Martin Ballaschk wrote: My solution is reading the files without the headers (skip = 1) and seperately reading the headers with scan (scan(myfile.CSV, what = character, sep = \t, nlines = 1). After throwing

Re: [R] Test of Homogeneity of Variances

Have you read the documentation to either of the functions you are using? ?bartlett.test Performs Bartlett's test of the null that the variances in each of the groups (samples) are the same. This explicitly tells you what is being tested, i.e. the null tested is that var1 = var2. ?rnorm

Re: [R] Test of Homogeneity of Variances

You don't need to test that the _sample_ variances are different: They already are. Statistical tests of hypotheses are not about sample statistics, but distributional charateristics. It seems to me that reinforcement of some basic stat concept may do you quite a bit of good. If you don't have

[R] inner margins for lattice

I would like to control the inner margins of a lattice graph. The graph is a single superposed panel, but there is too much white space around the data for my liking. I've played around with some of the layout options, but so far have found only how to control outer margins. R version 2.7.1

[R] R and GNUPLOT

Hello, I am trying to send commands to GNUPLOT to load a .plt file that was generated with a C++ software module called with R, and to replot. I am able to open the program with shell.exec (below) but I am at a loss at what do next. Any suggestions?? shell.exec(C:\\ gnuplot \\ bin \\

Re: [R] R and GNUPLOT

There are some interface commands in the TeachingDemos package (see ?gp.open) for communicating between R and Gnuplot. These commands are pretty basic, but looking at the code may give you some ideas of what to do next in your project (or maybe what not to do). Hope this helps, -- Gregory

[R] Sending ... to a C external

I'm trying to figure this out with Writing R Extensions but there's not a lot of detail on this issue. I want to write a (very simple really) C external that will be able to take ... as an argument. (It's for optimizing a function that may have several parameters besides the ones being optimized.)

Re: [R] lme questions re: repeated measures covariance structure

Great, thanks for your reply! I've just tracked down a copy of Pinheiro Bates (2000) so I'll look at that, too. Thanks again, Mark On Fri, Aug 22, 2008 at 9:56 AM, Christoph Scherber [EMAIL PROTECTED] wrote: Dear Mark, I would include the repeated measure as the smallest stratum in the

[R] sd(NA)

Hi All: This was discussed in the R-developers list (see the thread starting at http://tolstoy.newcastle.edu.au/R/e3/devel/ 07/12/0560.html). It has to do with the behavior of sd() when the entire vector is NA. The behavior has changed between 2.6 and 2.7.1 as follows: Run in

Re: [R] sd(NA)

Here's one approach: try_default - function(expr, default = NA, quiet = FALSE) { result - default if (quiet) { tryCatch(result - expr, error = function(e) {}) } else { try(result - expr) } result } failwith - function(default = NULL, f, quiet = FALSE) { function(...)

Re: [R] Writing Rcmdr Plugins

Bernhard and Irina, While tailoring the Rcmdr menu directly does work, I strongly do not recommend it. It will lead to mass confusion, most likely for the author, several months from now when it is installed on a different machine or when R or Rcmdr is updated. It is much cleaner to have a

Re: [R] History pruning

We have different starting points. Please be sure that your modularity allows a cleaned region as well as a history log to be the input to your next step. The history log is incomplete; lines sent to the *R* buffer by C-c C-n are explicitly excluded from history. Lines picked up from a saved

[R] Building colorspace on RHEL5

Dear all, I'm having problems installing the colorspace package on Red Hat Enterprise Linux 5: * Installing *source* package 'colorspace' ... ** libs gcc -std=gnu99 -I/usr/lib/R/include -I/usr/lib/R/include -I/usr/local/include-fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2 -fexceptions

Re: [R] Building colorspace on RHEL5

Hi Hadley, on 08/22/2008 03:52 PM hadley wickham wrote: Dear all, I'm having problems installing the colorspace package on Red Hat Enterprise Linux 5: * Installing *source* package 'colorspace' ... ** libs gcc -std=gnu99 -I/usr/lib/R/include -I/usr/lib/R/include -I/usr/local/include

[R] Newbie programming help

All - Not sure if this is a real programming question, but here goes: I have data that looks like LakeLength Weight 1 158 45 1 179 70 1 200 125 1 202 150 1 206 145 1 209 165 1 210 140 1 215 175 1 216

Re: [R] help needed for HWE.exact in library genetics

Sun, Ying asked: I have a genotype data for both case and controls and would like to calculate the HW p-value. However, since the number of one genotype is 0, I got wierd result. Would someone help me to figure it out? Or confirm it's right? Thanks a lot. ... HWE.exact(g1) Exact Test

Re: [R] Newbie programming help

On 8/22/2008 5:34 PM, Ranney, Steven wrote: All - Not sure if this is a real programming question, but here goes: I have data that looks like Lake Length Weight 1 158 45 1 179 70 1 200 125 1 202 150 1 206 145 1 209 165 1 210

[R] Using interactive plots to get information about data points

I have been experimenting with interactive packages such iplots and playwith. Consider the following sample dataset: A B C D 1 5 5 9 3 2 8 4 1 7 3 0 7 2 2 6 Let's say I make a plot of variable A. I would like to be able to click on a data point (e.g. 3) and have a pop-up window

[R] How to handle ~ character after csv importation

Dear R users, I have to import some csv files in which column headers contain the character ~. Following the following import call, the character seems to be replaced by dots in the column names of my data frame. Plus, I cannot query names(mydata) to find the column index which header should

[R] dimension values of an array

Question: I created an array with 164 columns and 70 rows. I would like every column to have unit value 1 and the rows 0.1, 0.2, 0.3 up to 7. How can I define that? Thanks for your help! Good weekend! Akko _

[R] nlme gls front end crash

Hi, when I try to run GLS, R crashes. I have seen on the mailing list that somebody has encountered the problem before ( http://tolstoy.newcastle.edu.au/R/help/05/01/10915.html ), but there seemed to be no solution / the thread appears incomplete. The program runs gls1 (without correlation

Re: [R] Visualization of two-mode matrix/ networks/ graphs

Jakob, Could you please tell me how you got your 2-mode data INTO R/statnet/sna, etc... ? I can do many things with this data in ucinet, and even export it to pajek, but I cannot get it to import even using the read.paj routine. I have about 620 actors involved in 230 events. When I convert

Re: [R] Sending ... to a C external

Le vendredi 22 août 2008 à 15:16 -0400, John Tillinghast a écrit : I'm trying to figure this out with Writing R Extensions but there's not a lot of detail on this issue. I want to write a (very simple really) C external that will be able to take ... as an argument. (It's for optimizing a

[R] Deleting NA in list()

Dear useRs, I would like to know the way of deleting NA in list(). Following is a example. lt - list(a = 1:3, b = NA, c = letters[1:3], d = NA) for(i in length(lt):1) { if(is.na(lt[[i]])) lt[[i]] - NULL } How to simplify for() loop by using (l)apply family? Thank you in advance.

[R] R

Hi my name is leidy. I am a student of sistem ingenier, i do not speak english. Now I am working in a proyect using R but with a interface of java, I have a problem because a need acces to p s-box.cox.powers(resp,2) s Box-Cox Transformation to Normality Est.Power Std.Err.

[R] problem with rbind

I am trying to use rbind to have the two data on the top of each other but I am getting an extra X on the column header and the rows are numberd , How to get rid of this problem? I appreciate your help x1- read.table(file=data1.txt, header=T, sep=\t) x2-read.table(file=data2.txt, header=T,

Re: [R] Simple look up

WTF mate. you think im about wasting time. I been seeking a solution for this for two weeks. excuse my lack of programming savvy, but this is why i posted. bartjoosen wrote: Please take a look at the manuals before posting (DO read the posting guide!) Bart PS:

  1   2   >