Re: [R] FW: logistic regression

Kjetil Halvorsen wrote: Hola! If the original questioner wants a guide as to how variables to measure IN THE FUTURE, when using his model in practice, thwen I think he will be unhappy with any advice which forces him to measure each of the 44 variables when probably a small subset will do!

[R] column-wise z-scores by group

Hi, I have a dataset of historical monthly temperature data that is grouped by weather station. I want to create z-scores of the monthly data using a base period of a subset of years. I subset the dataset first to include only data from the years (V2) that make up the base period so I could calcul

Re: [R] random normally distributed values within range

one could also use the "Truncate"-methods of package distr; cf. http://tolstoy.newcastle.edu.au/R/e5/help/08/09/1870.html If the situation is extreme Peter Daalgard gave some sophisticated code at http://tolstoy.newcastle.edu.au/R/e5/help/08/09/1892.html Best, Matthias Rolf Turner wrote: On 7

Re: [R] fitting a curve to data points

On 7/10/2008, at 2:53 PM, Carlos Guâno Grohmann wrote: Hello all. This is likely to be a silly question, but I have a set of data points and I want to fit a curve to it, like this: http://www.igc.usp.br/pessoais/guano/temp/curve.png. which function should I use? Have a look at smooth.spline(

Re: [R] Load a program at the front end

Em Qui, 2008-10-02 às 14:36 -0400, Gang Chen escreveu: > I want to run a R program, prog.R, interactively. My question is, is > there a way I can start prog.R on the shell terminal when invoking R, > instead of using source() inside R? > > TIA, > Gang Hi Gang I my system just only type: R --n

[R] The Innagural Atlanta useR Group Meeting!

Announcing a new Meetup for Atlanta useR Group (R Programming Language)! > > What: The Innagural Atlanta useR Group Meeting! > > When: November 5, 2008 7:00 AM > > Where: Click the link below to find out! > > Meetup Description: Hello All! > > I am putting this on the books for the first Atlanta us

[R] fitting a curve to data points

Hello all. This is likely to be a silly question, but I have a set of data points and I want to fit a curve to it, like this: http://www.igc.usp.br/pessoais/guano/temp/curve.png. which function should I use? many thanks Carlos -- +---+

[R] libg2c.so with new Debian

Hi R People: I just installed Debian on my desktop and am compiling R from source. However, I'm getting the error message: checking for Fortran 77 libraries of f77... checking how to get verbose linking output from gcc -std=gnu99... -v checking for C libraries of gcc -std=gnu99... -L/usr/local/

Re: [R] Automatic code diagramming for R?

This thread mentioning package mvbutils, function foodweb() might give you a start? http://tolstoy.newcastle.edu.au/R/help/06/04/24587.html Best, Erik Ben Bryant wrote: Greetings - Is anyone aware of an automatic code diagrammer/flow chart creator that works for the R language (either a cont

Re: [R] FW: logistic regression

Hola! If the original questioner wants a guide as to how variables to measure IN THE FUTURE, when using his model in practice, thwen I think he will be unhappy with any advice which forces him to measure each of the 44 variables when probably a small subset will do! What is wrong with first usin

Re: [R] Fisher inverse cumulative probabilities

On Mon, 6 Oct 2008, Raphael Saldanha wrote: Hi! How can I retrieve fisher's inverse cumulative probabilities, like, in Excel, with=INF(0.05,2;6) ? qf(0.95,2,6) [1] 5.143253 which is the same as =FINV(0.05,2,6) in Excel, presuming that is what you are after. Remember that on this list an

Re: [R] Automatic code diagramming for R?

findGlobals in the codetools package can find all the functions that a function calls and you could construct a call graph from that by combining that with RGraphviz in Bioconductor. For example, the proto package can automatically create a graph using RGraphviz showing the parent child relationsh

[R] Fisher inverse cumulative probabilities

Hi! How can I retrieve fisher's inverse cumulative probabilities, like, in Excel, with=INF(0.05,2;6) ? Regards, Raphael Saldanha UFJF - Brazil [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mail

[R] Automatic code diagramming for R?

Greetings - Is anyone aware of an automatic code diagrammer/flow chart creator that works for the R language (either a contributed package, or external software)? I need to explain some code structure of a package I'm working on to non-R users, and would find it extremely helpful to have such a p

Re: [R] random normally distributed values within range

On 7/10/2008, at 11:54 AM, Achaz von Hardenberg wrote: Hi all, I need to create 100 normally distributed random values (X) which can not exceed a specific range (i.e. 0 Presumably you want a truncated normal distribution. Duncan Murdoch posted some neat code on this list, back in the end of J

Re: [R] random normally distributed values within range

Hi Achaz, Maybe you are interested in the generalized beta distribution? To the best of my knowledge, there is no way to restrict the values of normal deviates, since one may end up with a different distribution. Regards, Pedro -Original Message- From: [EMAIL PROTECTED] [mailto:[

Re: [R] querry in simple "for" loop in Box plot

worked fine for me (after removing 'data=xx' which may have caused caused problems since it is supposed to be a dataframe or a list, not a matrix) > par(mfrow=c(3,2)) > for(i in 2:7) + { + boxplot(xx[,i]~xx[,1],xlab="lev") + print(i) + } [1] 2 [1] 3 [1] 4 [1] 5 [1] 6 [1] 7 On Mon, Oct 6, 2008 at

[R] random normally distributed values within range

Hi all, I need to create 100 normally distributed random values (X) which can not exceed a specific range (i.e. 0With rnorm I cannot specify Max and min values among which values have to stay, like in runif so does some other simple way exist to do this with normally distributed random valu

Re: [R] Graphics window BUG

I think I managed to figure out the problem (for the record I use Windows XP 2002): One is required to set options in Edit-GUI Preferences, then change "Pager style" to single window. In the graphics window (x11()) choose History-Clear History, then History-Record each time you wish to generate

[R] R on the Freakonomics Blog

http://freakonomics.blogs.nytimes.com/2008/10/06/free-super-crunching-so ftware/ >>> This e-mail and any attachments are confidential, may contain legal, professional or other privileged information, and are intended solely for the addressee. If you are not the intended recipient, do not use the

[R] Graphics window BUG

I thought this problem would be resolved when I switched to R version 2.7.0 (for Windows), but no - anytime I plot something that produces more than one page of graphics, the graphics window starts by showing the first page, until such time as I hit enter to show me the next page, at which time it

Re: [R] Sweave and echoing code chunks conditionally

On Mon, Oct 6, 2008 at 3:30 PM, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > Try > > \SweaveOpts{echo=FALSE} > or > \SweaveOpts{echo=TRUE} Perfect (and obvious). Thanks. Sean > to set the default for echo one way or the other. > > On Mon, Oct 6, 2008 at 3:27 PM, Sean Davis <[EMAIL PROTECTED

Re: [R] Sweave and echoing code chunks conditionally

Try \SweaveOpts{echo=FALSE} or \SweaveOpts{echo=TRUE} to set the default for echo one way or the other. On Mon, Oct 6, 2008 at 3:27 PM, Sean Davis <[EMAIL PROTECTED]> wrote: > I have started to use Sweave to prepare reports for "end-users" of > statistical analyses. For these users, there is no

[R] Sweave and echoing code chunks conditionally

I have started to use Sweave to prepare reports for "end-users" of statistical analyses. For these users, there is not generally interest in the code used to accomplish the analysis, but only the results. However, sometimes I would like to prepare the same document, but with code chunks echoed, f

Re: [R] 3 curves / 1 plot

You can use the points() and lines() functions to add points and lines to an existing plot. Julian Michel PETITJEAN wrote: I am a new user of R. Please does somebody knows how to plot 3 datasets (x1,a1),...,(xn,an), (x1,b1),...,(xn,bn), and (x1,c1),...,(xn,cn) on a single x,y plot, each of the

Re: [R] Lattice and Axis height

On Mon, Oct 6, 2008 at 10:37 AM, Saptarshi Guha <[EMAIL PROTECTED]> wrote: > Hello, > My plots consist of 2 rows, 1 column, many pages > I've managed to selectively turn of strips for the bottom panel and > roll my own strip for the top > > par.strip.text = list(lines = c(0,3),cex=0.6) > > I have

Re: [R] Color of title in legend()

not reproducible On Mon, Oct 6, 2008 at 2:01 PM, Michael Kubovy <[EMAIL PROTECTED]> wrote: > Dear R-helpers, > > This produces a red title. How do I make it black, without making the > labels black? > > legend(title = 'Channel Integration', 0.7, 0.3, > c('loudness','luminance'), col = c('red', 'bl

Re: [R] package question

library() viginette's usually describe how to use the packages with examples. On Mon, Oct 6, 2008 at 2:24 PM, ANJAN PURKAYASTHA <[EMAIL PROTECTED]> wrote: > hi, > am new to R. what command do i use to check of certain packages have been > loaded? > Also, it seems that most packages come with "vign

[R] querry in simple "for" loop in Box plot

Dear R users, I am new users of this software. I want to make box plot. Here, i have simulated data set with following commands: x<-matrix(rnorm(90),nrow=10,ncol=9) x a<-matrix(c(1,1,1,1,1,2,2,2,2,2),nrow=10,ncol=1) xx<-cbind(a,x) colnames(xx)<-c("a","b","c","d","e","f","g","h","i","j") rownames

[R] package question

hi, am new to R. what command do i use to check of certain packages have been loaded? Also, it seems that most packages come with "vignette"- is this a document that describes the usage of the package? thanks for your help, anjan -- = anjan purkayastha, phd bioinformat

Re: [R] question on lmList

eugen pircalabelu yahoo.com> writes: > Using the lmList function from “nlme package” I get the following error >> message: > > “Error in !unlist(lapply(sum.lst, is.null)) : invalid argument type” > try str(zz) to check if your data are strange. And plot the data, maybe there is something fun

[R] Color of title in legend()

Dear R-helpers, This produces a red title. How do I make it black, without making the labels black? legend(title = 'Channel Integration', 0.7, 0.3, c('loudness','luminance'), col = c('red', 'blue'), text.col = c('red', 'blue'), lty = 1, xjust = 0.5, bty = 'n') _

Re: [R] Stacked area chart and legends

On Mon, Oct 6, 2008 at 10:54 AM, mbr <[EMAIL PROTECTED]> wrote: > > I have found lots of good advice on this forum about stacked area charts but > I've run into problems with the 2 recommended options: stackploy in plotrix > or qplot in ggplot2. I have a many page report that will be in a 2x2 page

[R] Lattice and Axis height

Hello, My plots consist of 2 rows, 1 column, many pages I've managed to selectively turn of strips for the bottom panel and roll my own strip for the top par.strip.text = list(lines = c(0,3),cex=0.6) I have my own custom prepanel function too. What i'd like to do is turn off the x-axis for th

[R] Need to calculate within- and between- run CV

Dear R-helpers, I have a dataset named "qu", organized as follows: Sample Run Replicate Value 1 1 1 25 1 1 2 40 1 1 3 33 1 1 4 29 1 2 1 37 1 2

Re: [R] Convert Binary to Decimal

"Jason Thibodeau" <[EMAIL PROTECTED]> wrote in message news:[EMAIL PROTECTED] > Hello, > > A few weeks ago I asked about converting a decimal number to binary. This > time, I need to do the opposite, take a binary number, and convert ti to > its > decimal equivalent. I have googled, but didn't fi

[R] Computing on the language redux

Warning: This is only for those interested in R language minutiae A recent post on this list asked if there was a simple way to change the R language object: ex1 <- expression(x < a) ## just the part to the right of the <- assignment to the object expression( x < a & y < b) ## or something l

[R] Convert Binary to Decimal

Hello, A few weeks ago I asked about converting a decimal number to binary. This time, I need to do the opposite, take a binary number, and convert ti to its decimal equivalent. I have googled, but didn't find much in the way of elegant solutions. Since you all were such a help last time, I figure

Re: [R] easier way to do this without a loop? (successive euclidean distances between points)

Or just: diag(rdist(z, z)[,-1]) On Mon, Oct 6, 2008 at 1:58 PM, Henrique Dallazuanna <[EMAIL PROTECTED]> wrote: > Try this: > > diag(rdist(rbind(z, 0), rbind(0, z))[,-(1:2)]) > > On Mon, Oct 6, 2008 at 1:49 PM, stephen sefick <[EMAIL PROTECTED]> wrote: >> a <- c(1:10) >> b <- c(.5, .6, .9, 10, .4

Re: [R] easier way to do this without a loop? (successive euclidean distances between points)

Try this: diag(rdist(rbind(z, 0), rbind(0, z))[,-(1:2)]) On Mon, Oct 6, 2008 at 1:49 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > a <- c(1:10) > b <- c(.5, .6, .9, 10, .4, 3, 4, 9, 0, 11) > d <- c(21:30) > > z <- data.frame(a,b,d) > library(fields) > results <- c() > for(i in 1:(length(rowname

[R] easier way to do this without a loop? (successive euclidean distances between points)

a <- c(1:10) b <- c(.5, .6, .9, 10, .4, 3, 4, 9, 0, 11) d <- c(21:30) z <- data.frame(a,b,d) library(fields) results <- c() for(i in 1:(length(rownames(z))-1)){ results[i] <- rdist(z[i,], z[(i+1),]) } results.1 <- data.frame(results) f <- rownames(z) r <- f[-1] rownames(results.1

[R] Computationally singular [provides coefficients but not covariance matrix]

Hi, I am estimating a regression but the summary command is unable to provide me results, while the coefficients are available from the coefficients value. I suppose that it cannot estimate the covariance matrix. Is there any command that I can relax the tolerance so it can estimate the covarianc

Re: [R] stepplr

stephen sefick wrote: you are going to have to install it from source- look at the type argument in ?install.packages. There is no windows binary. The reason is that stepplr fails the checks under Windows, hence we do not provide a binary version. Please ask the package maintainer to fix

Re: [R] stepplr

you are going to have to install it from source- look at the type argument in ?install.packages. There is no windows binary. On Mon, Oct 6, 2008 at 10:35 AM, Samor Gandhi <[EMAIL PROTECTED]> wrote: > Hello everybody, > > I am trying to install the library stepplr under windows > (http://www.mat

[R] Stacked area chart and legends

I have found lots of good advice on this forum about stacked area charts but I’ve run into problems with the 2 recommended options: stackploy in plotrix or qplot in ggplot2. I have a many page report that will be in a 2x2 page format “par(mfrow=c(2,2))”” and need one of the page components to be a

[R] Extracting text from html code using the RCurl package.

Dear R-help, I want to download the text from a web page, however what i end up with is the html code. Is there some option that i am missing in the RCurl package? Or is there another way to achieve this? This is the code i am using: > library(RCurl) > my.url <- 'https://stat.ethz.ch/mailman/list

[R] Error in X11

Hello everyone, I'm trying to plot a graphic in Linux, when I type X11() then I have an error which is the next: Error in X11(d$display, d$width, d$height, d$pointsize, d$gamma, d$colortype,  :  unable to start device X11cairo   Why?, what I must do to fix it?. Thanks so much   Carlos _

Re: [R] horizontal boxplot + xlim

> True, I made a mistake here. Still, I have problems to visualize my data (not > the example code I used). I just see a flat line instead a proper plot... > > Another example code with creating a strange plot: > > x <- rnorm(100) + 100 > maxval <- max(x) > boxplot(x, notch=TRUE, xlim = c(0,ma

[R] stepplr

Hello everybody, I am trying to install the library stepplr under windows (http://www.maths.bris.ac.uk/R/web/packages/stepPlr/index.html), in order to use the function plr, but I still have problem to find the right link for this purpose! I am very thankful for your help! Samor __

Re: [R] lmer: random factor nested in a fixed factor

Dear Agnes, I think your model specification should look like this: YourModel1 <- lmerlmer(y ~ poptype*matingtype + (1|poptype:pop) + (1|poptype:fam), data = ...) The "1" in front of "|" refers to models that are random intercepts models as opposed to general random coefficients models in whi

Re: [R] horizontal boxplot + xlim

[EMAIL PROTECTED] schrieb: I get a strange behaviour of a boxplot with the following code. There seems to be a problem with the xlim-parameter. Did I do anything wrong? What else can I do to force the boxplot to have a defined x-range? x <- rnorm(100) boxplot(x, notch=TRUE, xlab=parameter, x

[R] question on lmList

Hi list, Using the lmList function from “nlme package” I get the following error message: “Error in !unlist(lapply(sum.lst, is.null)) : invalid argument type” # this is the syntax used # fm_cream <-lmList(MULTDV~TIME|cod, data=zz, na.action=na.omit) # summary(fm_cream) When I generate a si

[R] Rmpi memory error

Hi List, I am testing to see if Rmpi work on our cluster. I load Rmpi and try a simple test: > library(Rmpi) > mpi.spawn.Rslaves(nslaves = 4) *** caught segfault *** address 0x248cba8, cause 'memory not mapped' Traceback: 1: .Call("mpi_comm_spawn", as.character(slave), as.character(slavear

[R] Query about computational demand

Hi all; I've programmed a couple of C libraries which are loaded dynamically into R (Linux). With one of these, I'm conducting Monte Carlo analysis, but every individual execution of my model is about 15'. So, I'm running 1000 executions in about 11 days. This is not enough for my needings, as I n

[R] lme and lmer df's and F-statistics again

Dear R-users, I did do a thorough search and read many articles and forum threads on the lme and lmer methods and their pitfalls and problems. I, being not a good statistician but a mere "user", came to the conclusion that the most correct form of reporting statistics for a mixed linear model wou

Re: [R] horizontal boxplot + xlim

> I get a strange behaviour of a boxplot with the following code. > There seems to > be a problem with the xlim-parameter. Did I do anything wrong? What > else can I > do to force the boxplot to have a defined x-range? > > x <- rnorm(100) > boxplot(x, notch=TRUE, xlab=parameter, xlim <- c(-4,4

Re: [R] Read.table fail to recognize (+/-) sign

I am able to replicate your problem, using R 2.7.2: R version 2.7.2 (2008-08-25) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 ... > source <- read.table("clipboard") > source V1 V2V3 V4 V5 V6V7 1 HIT1 hg18 chr8 + 175

[R] lmer: random factor nested in a fixed factor

Hi all, I try to build a model using lmer, with 2 crossed fixed factors (poptype and matingtype) and 2 random factors (pop and family) which I want to nest within poptype. Which of these formulae should I use ? 1. lmer(y~poptype*matingtype + (1|poptype/pop/fam)) 2. lmer(y~poptype*matingtype +

[R] horizontal boxplot + xlim

Hi there, I get a strange behaviour of a boxplot with the following code. There seems to be a problem with the xlim-parameter. Did I do anything wrong? What else can I do to force the boxplot to have a defined x-range? x <- rnorm(100) boxplot(x, notch=TRUE, xlab=parameter, xlim <- c(-4,4), ho

[R] redirected to r-help (was RE: metafiles under windows)

My experience is that this could happen when I switch between an external display and the laptop LCD that have different resolution. Can't seem to fix it short of reboot. If anyone has better solution, I'm all ears. Best, Andy From: [EMAIL PROTECTED] > > Probably not a bug, not tested on a re

Re: [R] Read.table fail to recognize (+/-) sign

Please do show us the courtesy of following the posting guide, and make sure you use a current version of R. R 2.7.2 and 2.8.0 beta do not do this, and you were asked to update *before* posting. [My guess is that you are using the long-obsolete 2.7.0.] On Mon, 6 Oct 2008, Gundala Viswanath wr

Re: [R] aggregate empty row for pretty appearance also subtotal if possible

Thanks Jim. Dhruv -Original Message- From: jim holtman [mailto:[EMAIL PROTECTED] Sent: Friday, October 03, 2008 10:41 PM To: Sharma, Dhruv Cc: r-help@r-project.org Subject: Re: [R] aggregate empty row for pretty appearance also subtotal if possible You can try reshape: > n <- 2000 >

[R] Odp: Read.table fail to recognize (+/-) sign

Hi works for me, which version of R? > read.table("clipboard") V1 V2V3 V4 V5 V6V7 1 HIT1 hg18 chr8 + 1759549 1894206 Human 2 HIT5 hg18 chr11 - 8593662 8636959 Human 3 HIT7 hg18 chr5 + 31675274 32146794 Human 4 HIT00011 hg18 c

[R] Read.table fail to recognize (+/-) sign

Dear all, I have the following dataset: # HIT or acc, UCSC genome, chromosome, strand, start, end, species name HIT1hg18chr8+ 1759549 1894206 Human HIT5hg18chr11 - 8593662 8636959 Human HIT7hg18chr5+ 3167527432146794Human HIT

Re: [R] t.test() on a list

Hi maybe not an answer you like but apply(do.call(rbind, lapply(myList, c)), 2, t.test) shall give you desired results Regards Petr [EMAIL PROTECTED] napsal dne 02.10.2008 23:25:23: > I appreciate your suggestion. The example I provided was fabricated > because I was only focusing on the pro

Re: [R] ignore error in for-loop

Erich STRIESSNIG wrote: > Hi all, > > Can anybody tell me what to do in case an error occurs inside a > for-loop and I don't want the program to exit the loop? Rather instead > I would like it to just go to the next number in the loop and try > again with the new number. Is there any function like

[R] Odp: an unsophisticated question about recoding in a data frame with control structure if {}

Hi yourm question has two aspects. one is simple data.frame$thevector[dataframe$factor=='3'] <- an arithmetic.mean > zdrz doba otac sklon 1 189.0 0.6 110 2 256.0 0.680 3 286.0 0.650 4 105.0 1.250 zdrz$otac==0.6 zdrz$sklon[zdrz$otac==0.6]<-11 > zdrz doba otac sklon

Re: [R] ignore error in for-loop

Dear Erich, Have a look at the try()-function. HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie, methodologie en kwaliteitszorg / Sect

[R] ignore error in for-loop

Hi all, Can anybody tell me what to do in case an error occurs inside a for-loop and I don't want the program to exit the loop? Rather instead I would like it to just go to the next number in the loop and try again with the new number. Is there any function like "on error go back to" or "skip

Re: [R] Problem with Grep Under Loop

On Mon, Oct 6, 2008 at 2:16 PM, Gundala Viswanath <[EMAIL PROTECTED]> wrote: > Dear all, > > This is the sample of the source data (tab delimited) > > http://dpaste.com/82660/plain/ > > - Gundala Viswanath > Jakarta - Indonesia > Hi, using this data, your original code works flawlessly - nothing

Re: [R] Problem with Grep Under Loop

Dear all, This is the sample of the source data (tab delimited) http://dpaste.com/82660/plain/ - Gundala Viswanath Jakarta - Indonesia On Mon, Oct 6, 2008 at 9:05 PM, Gustaf Rydevik <[EMAIL PROTECTED]> wrote: > On Mon, Oct 6, 2008 at 1:37 PM, Gundala Viswanath <[EMAIL PROTECTED]> wrote: >> De

Re: [R] Problem with Grep Under Loop

On Mon, Oct 6, 2008 at 1:37 PM, Gundala Viswanath <[EMAIL PROTECTED]> wrote: > Dear all, > > I have no problem with this individual grep command: > >> datk <- grep("XM_528056", source$V1) >> dat2 <- source[datk,] >> print(dat2) > V1 V2 V3 V4 V5 V6 V7 > 35995 XM_

Re: [R] NA's in segmented

Dear Tyler, Yes the problem is with NA.. There are two solutions: 1) You can use lm() + segmented (you fit a gaussian model, so why do you use glm()?) 2)If you want to use glm()+ segmented(), use na.omit() to pass your dataframe to the data argument of glm, glm(.., data=na.omit()) Also, if

Re: [R] Problem with Grep Under Loop

Seems to work fine for me. You need to read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. You did not provide 'source', so I created one and it seems to work. So I have to conclude you have something wrong in your d

Re: [R] color code from csv

Paul Hiemstra skrev: kerfuffle wrote: hi folks, this is driving me up the wall. Apologies for posting twice in the same week, I'm writing up a thesis. I wish to color-code some dots in an xy plot. I've got a csv file with various elements, one of which is the color-key (with the header 'co

Re: [R] Margin on trellis plot

> I have a problem when I make a trellis plot - I would like to reduce > the white margin surrounding the plot. > > I used S some years age, as I remember it was possible to use the > "par" and "oma" to control margin size even on a trelllis plot (see > the plot below). > > However, it does not

[R] Problem loading package created with package.skeleton

Hi, i'm trying to build an r-package (Windows Vista, R 2.7.2), so i created one with the package.skeleton() command. After that i zipped it and tried to load it in R. Attaching the package via the Windows R-console menu 'Packages/Install package(s) from local zip-files' failed: /> utils:::menuIn

Re: [R] color code from csv

kerfuffle wrote: hi folks, this is driving me up the wall. Apologies for posting twice in the same week, I'm writing up a thesis. I wish to color-code some dots in an xy plot. I've got a csv file with various elements, one of which is the color-key (with the header 'color'). If the color-key

[R] Problem with Grep Under Loop

Dear all, I have no problem with this individual grep command: > datk <- grep("XM_528056", source$V1) > dat2 <- source[datk,] > print(dat2) V1 V2 V3 V4 V5 V6 V7 35995 XM_528056 panTro2 chr8 + 1775569 1896107 Chimpanzee BUT, when I run them under the loop i

Re: [R] propensity score adjustment using R

the Matching() package by Jasjeet Sekhon does propensity score matching in a very user friendly way. (as you said you don't want to reinvent the wheel...) just feed it with the fitted values from a glm model (fitted$myglmmodel). afaik, you may additionally match on some covariates directly. HTH m

[R] Margin on trellis plot

Dear all I have a problem when I make a trellis plot - I would like to reduce the white margin surrounding the plot. I used S some years age, as I remember it was possible to use the "par" and "oma" to control margin size even on a trelllis plot (see the plot below). However, it does not seem

Re: [R] Bioconductor installation on 64 bit R problem

Bioconductor has its own list, but note that many of us (including the Bioconductor test farm) are happliy testing all the Bioconductor packages on 64-bit platforms (including Linux and Solaris -- some at least have been used on 64-bit Mac OS 10.5). So please read the posting guide and follow

Re: [R] building packages: "R Help for package foo" vs. "HTML Help"?

Uwe Ligges wrote: Michael Friendly wrote: Uwe Ligges wrote: Michael Friendly wrote: In building a package, what are the settings in the package files or the build commands that determine whether the compiled HTML help windows have the window title "R Help for package foo" vs. "

[R] R 2.7.2 upgrade in Ubuntu, tcltk does not work any more

Hello, I had no problems running R with tcltk until the recent upgrade to R-2.7.2-2gustsy1 (I use Ubuntu 7.10 Gutsy Gibbon). Since the upgrade from R-2.7.2-1, I got this error (see hereunder). Obviously, I have now a problem with tcltk-related .so files. What should I do? Thanks. Philippe [

Re: [R] Error in Q-Q plot

Megh Dal yahoo.com> writes: > > > library(sn) > > library(car) > > dat1 = rst(1000, 0, 1, 0, 2) > > qq.plot(dat1, "st", 0, 1, 0, 9) > Error in plot.window(...) : > invalid value specified for graphical parameter "las" You have set parameter las (the 6th) to 9, and only values 0 and 1 are all

[R] Bioconductor installation on 64 bit R problem

Hello, I am attempting to install Bioconductor on 64 bit R. The installation of some bioconductor packages fails with the errors below. I would be very grateful for any help on how to resolve the issue. Best regards, Octavio Espinosa - Most of the bioconductor packages

[R] RBGL installation

Hi I'm trying to install the RBGL package. However I am getting the following error; can anybody shed some light on what I'm missing? /usr/lib/gcc/i486-linux-gnu/4.1.2/../../../../include/c++/4.1.2/cstdlib:135: error: ‘::system’ has not been declared Thanks. Iain > sessionInfo() R vers

[R] package flexmix: models for ordinal data?

Dear List, does anyone have experience with fitting finite mixture regression models to ordinal/ordered dependent data using Friedrich Leisch's flexmix package? Best, Klaus __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

[R] Error in Q-Q plot

Hi, I am trying to draw a Q-Q plot, however got following error. > library(sn) > library(car) > dat1 = rst(1000, 0, 1, 0, 2) > qq.plot(dat1, "st", 0, 1, 0, 9) Error in plot.window(...) : invalid value specified for graphical parameter "las" Can anyone please tell me why this error is coming?

Re: [R] Tinn-R command line argument to increase memory limits.

You can add the arguments here: Options --> Application --> R --> Path Regards, Yihui -- Yihui Xie <[EMAIL PROTECTED]> Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.name School of Statistics, Room 1037, Mingde Main Building, Renmin University

[R] Tinn-R command line argument to increase memory limits.

I searched the FAQ and found that if I set --max-mem-size=512M (or above) on the command line I can increase the. What I was unable to find in the archives was how to pass this command line argument when using Tinn-R. I know I have seen this just recently. Any suggestions? Thank you. Kevin __

Re: [R] building packages: "R Help for package foo" vs. "HTML Help"?

Michael Friendly wrote: Uwe Ligges wrote: Michael Friendly wrote: In building a package, what are the settings in the package files or the build commands that determine whether the compiled HTML help windows have the window title "R Help for package foo" vs. "HTML Help"? Michael, can yo

Re: [R] How to model this? (R-help Digest, Vol 67, Issue 25)

Hello Senne: I am currently teaching a simulation and modeling course in which we use R for most purposes and Sim Py (a python module for discrete event simulation). I am looking for a real-life case such as yours to approach in class- if you wish, i would be happy to look into the case you ar