[R] odfWeave cross-reference

2007-12-21 Thread Chris.H. Snow
How can I insert cross-references to odfWeave generated figures in my source odf before the graphic has been created with odfWeave? Many thanks, Chris __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] R appearance under linux

2007-12-21 Thread Tyler Smith
On 2007-12-21, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > > Dear R users, > > I have just moved to R2.6.1 under Opensuse linux 10.3. I used to > work with R under XPpro. Is it "normal" to have a visual aspect of R > under linux different ? Yes, that's normal. Under windows, you get a GUI int

[R] regarding lack of quadratic term in solve.qp

2007-12-21 Thread Mark Leeds
I was thinking about my solve.qp problem on my way home tonight and I think I can fix it by making Dmat the identity matrix. I'll check this weekend to make sure but my thinking is that doing this will make all the "variances" the same so that they shouldn't come into play during the evaluation of

Re: [R] Apparent discordant solution using NLS() vs Optim or Excel for IWLS problem

2007-12-21 Thread Gabor Grothendieck
Check out the drc package. Also, please read the last line to every message to r-help regarding reproducible code and note that function names are case sensitive in R so NLS is not the same as nls. On Dec 21, 2007 7:36 PM, <[EMAIL PROTECTED]> wrote: > I am writing a program for automated (i.e. n

[R] mle

2007-12-21 Thread Antonio Gasparrini
Dear all, I'm trying to estimate the parameters of a special case of a poisson model, where the specified equation has an integral and several fixed parameters. I think that the MLE command in STATS4 package could be a good choice, but it's a little complicated. I've got some problems with the off

[R] mle

2007-12-21 Thread Antonio Gasparrini
Dear all, I'm trying to estimate the parameters of a special case of a poisson model, where the specified equation has an integral and several fixed parameters. I think that the MLE command in STATS4 package could be a good choice, but it's a little complicated. I've got some problems with the off

Re: [R] Efficient way to find consecutive integers in vector?

2007-12-21 Thread Henrik Bengtsson
In the R.utils package there is seqToIntervals(), e.g. print(seqToIntervals(1:10)) ## from to ## [1,]1 10 print(seqToIntervals(c(1:10, 15:18, 20))) ## from to ## [1,]1 10 ## [2,] 15 18 ## [3,] 20 20 There is also seqToIntervals(), which uses the above, e.g. print(seqToHuman

[R] using solve.qp without a quadratic term

2007-12-21 Thread markleeds
I was playing around with a simple example using solve.qp ( function is in the quadprog package ) and the code is below. ( I'm not even sure there if there is a reasonable solution because I made the problem up ). But, when I try to use solve.QP to solve it, I get the error that D in the quad

Re: [R] Efficient way to find consecutive integers in vector?

2007-12-21 Thread Tony Plate
Martin Maechler wrote: >> "MS" == Marc Schwartz <[EMAIL PROTECTED]> >> on Thu, 20 Dec 2007 16:33:54 -0600 writes: > > MS> On Thu, 2007-12-20 at 22:43 +0100, Johannes Graumann wrote: > >> Hi all, > >> > >> Does anybody have a magic trick handy to isolate directly consec

Re: [R] assigning and saving datasets in a loop, with names changing with "i"

2007-12-21 Thread Tony Plate
Marie Pierre Sylvestre wrote: > Dear R users, > > I am analysing a very large data set and I need to perform several data > manipulations. The dataset is so big that the only way I can play with it > without having memory problems (E.g. "cannot allocate vectors of size...") > is to write a batch s

Re: [R] Numerical precision problem

2007-12-21 Thread dlakelan
[EMAIL PROTECTED] wrote: > Dear All, > >> (1-p)==1 > > [1] TRUE > > The main problem I have is that I need to evaluate the inverse of the normal > Cumulative Distribution Function near 1 (but not exactly 1) and the PC cannot > recognize it. Evaluate it at 1, evaluate the density at 1, an

[R] Apparent discordant solution using NLS() vs Optim or Excel for IWLS problem

2007-12-21 Thread rsposto
I am writing a program for automated (i.e. no user intervention - for biologists) iteratively reweighted least-square fit to a function of the form "reading ~ exp(lm2)/(1 + (dose/exp(lm3))^exp(lm4)" using case weight proportional to the mean, e.g., E(reading).Because for some datasets the s

Re: [R] New Cran Task View: ExperimentalDesign

2007-12-21 Thread Ulrike Grömping
Thanks, Hadley! Regards, Ulrike -- Original Message --- From: "hadley wickham" <[EMAIL PROTECTED]> To: "Ulrike Grömping" <[EMAIL PROTECTED]> Cc: r-help@r-project.org Sent: Fri, 21 Dec 2007 22:49:38 + Subject: Re: [R] New Cran Task View: ExperimentalDesign > And the url is:

Re: [R] Question about which kind of plot to use

2007-12-21 Thread Max
hadley wickham presented the following explanation : >> Perhaps as long as you're learning a new plotting system, you might also >> check out whether ggplot2 might be an option. >> >> I did a quick and dirty version (which I'm sure Hadley can improve and >> also remind me how to get rid of the leg

Re: [R] Finding overlaps in vector

2007-12-21 Thread Gabor Grothendieck
If we don't need any plotting we don't really need rect.hclust at all. Split the output of cutree, instead. Continuing from the prior code: > for(el in split(unname(vv), names(vv))) print(el) [1] 0.00 0.45 [1] 1 [1] 2 [1] 3.00 3.25 3.33 3.75 4.10 [1] 5 [1] 6.00 6.45 [1] 7.0 7.1 [1] 8 On Dec 21,

Re: [R] Diagonal matrix with off diagonal elements

2007-12-21 Thread Chris Stubben
Also try the odiag function in the demogR package odiag( 1:5, -1) [,1] [,2] [,3] [,4] [,5] [,6] [1,]000000 [2,]100000 [3,]020000 [4,]003000 [5,]000400 [6,]000

Re: [R] New Cran Task View: ExperimentalDesign

2007-12-21 Thread hadley wickham
And the url is: http://cran.r-project.org/src/contrib/Views/ExperimentalDesign.html Hadley On 12/21/07, Ulrike Grömping <[EMAIL PROTECTED]> wrote: > Dear UseRs, > > the new Task View "ExperimentalDesign" (Title: Design of Experiments (DoE) & > Analysis of Experimental Data) has just been uploaded

Re: [R] Diagonal matrix with off diagonal elements

2007-12-21 Thread Gabor Grothendieck
On Dec 21, 2007 3:58 PM, Jonas Malmros <[EMAIL PROTECTED]> wrote: > I wonder if there is a function in R with which I can create a square > matrix with elements off main diagonal (for example one diagonal below > the main diagonal). Try this and adjust the formula for other patterns: > mm <- matr

[R] RMySQL installation

2007-12-21 Thread Janet Young
Hi, I am having trouble getting RMySQL running on a solaris machine. [43] bedrock:/home/jayoung/source_codes/R/other_packages> uname -a SunOS bedrock 5.10 Generic_118833-36 sun4v sparc SUNW,Sun-Fire-T200 I thought I had finally managed to get it installed, albeit with some warnings that I didn'

Re: [R] Numerical precision problem

2007-12-21 Thread Thomas Lumley
On Fri, 21 Dec 2007, [EMAIL PROTECTED] wrote: > Dear All, > > > > I have a numerical problem: > > R, as other statistical software, can tell the difference between very > small numbers and 0, while apparently cannot distinguish a number close > to 1 and 1. In the example below, is it possible to

[R] New Cran Task View: ExperimentalDesign

2007-12-21 Thread Ulrike Grömping
Dear UseRs, the new Task View "ExperimentalDesign" (Title: Design of Experiments (DoE) & Analysis of Experimental Data) has just been uploaded to CRAN. It covers packages that have functionalities for designing experiments or specific tools for analyzing data from designed experiments. I hope that

Re: [R] Diagonal matrix with off diagonal elements

2007-12-21 Thread Domenico Vistocco
Jonas Malmros wrote: > Hi, everyone > > I wonder if there is a function in R with which I can create a square > matrix with elements off main diagonal (for example one diagonal below > the main diagonal). > > Thanks in advance! > > You could combine rbind, diag and cbind: rbind(rep(0,3),cbind(d

Re: [R] Forestplot

2007-12-21 Thread francogrex
Sorry here's the good code: forestplot=function(mean,std,conf,threshold){ z=-qnorm((1-conf)/2) CI.H <- mean+z*std # Calculate upper CI CI.L <- mean-z*std # Calculate lower CI plot(mean,1:length(mean),xlim=c(min(CI.L),max(CI.H)),pch=15,cex=1, ylab="",xlab="mean & CIs",main="Forest plot") arrows(

Re: [R] post hoc in repeated measures of anova

2007-12-21 Thread Dieter Menne
Kazumi Maniwa uni-konstanz.de> writes: > > Hallo, I have this dataset with repeated measures. There are two > within-subject factors, "formant" (2 levels: 1 and 2) and "f2 Ref" (25 > levels: 670, 729, 788, 846, 905, 1080, 1100, 1120, 1140, 1170, 1480, > 1470, 1450, 1440, 1430, 1890, 1840,

[R] Diagonal matrix with off diagonal elements

2007-12-21 Thread Jonas Malmros
Hi, everyone I wonder if there is a function in R with which I can create a square matrix with elements off main diagonal (for example one diagonal below the main diagonal). Thanks in advance! -- Jonas Malmros Stockholm University Stockholm, Sweden _

Re: [R] Forestplot

2007-12-21 Thread francogrex
Answering my own query here's a solution (if any reader was curious): forestplot=function(mean,std,conf,threshold){ z=-qnorm((1-conf)/2) CI.H <- mean+z*sqrt(std) # Calculate upper CI CI.L <- mean-z*sqrt(std) # Calculate lower CI plot(mean,1:length(mean),xlim=c(min(CI.L),max(CI.H)),pch=15,cex=1,

Re: [R] Finding overlaps in vector

2007-12-21 Thread Johannes Graumann
Hm, hm, rect.hclust doesn't accept "plot=FALSE" and cutree doesn't retain the indexes of membership ... anyway short of ripping out the guts of rect.hclust to achieve the same result without an active graphics device? Joh >> # cluster and plot >> hc <- hclust(dist(v), method = "single") >> plot(h

Re: [R] Finding overlaps in vector

2007-12-21 Thread Johannes Graumann
Jim, Although I can't find the post this code stems from, I had come across it on my prowling the NG. It's not the one you had shared with me to eliminate overlaps (and which I referenced below: http://tolstoy.newcastle.edu.au/R/e2/help/07/07/21286.html). That particular solution you had come up w

Re: [R] Finding overlaps in vector

2007-12-21 Thread Johannes Graumann
Thank you very much for this elegant solution to the problem. The reason I still hope for an extension of Jim's code (not the one re responded with in this thread, but the one I actually reference) is that windows of overlap can be asymetric with that: one can check e.g. whether values overlap give

Re: [R] Numerical precision problem

2007-12-21 Thread Duncan Murdoch
On 21/12/2007 2:52 PM, [EMAIL PROTECTED] wrote: > Dear All, > > > > I have a numerical problem: > > R, as other statistical software, can tell the difference between very small > numbers and 0, while apparently cannot distinguish a number close to 1 and 1. > In the example below, is it possi

[R] Numerical precision problem

2007-12-21 Thread C.Rosa
Dear All, I have a numerical problem: R, as other statistical software, can tell the difference between very small numbers and 0, while apparently cannot distinguish a number close to 1 and 1. In the example below, is it possible to instruct R to recognize that q is different from 1? > p=.

Re: [R] adding lines to a barchart

2007-12-21 Thread Deepayan Sarkar
On 12/21/07, Spilak,Jacqueline [Edm] <[EMAIL PROTECTED]> wrote: > Thanks Deepayan, this helped a lot and gave me exactly what I wanted > however I want a few changes and am not sure how to do them. In each > panel the bars for the years are the same color with the bottom axis > (the x-axis) labell

[R] matrix to gal object

2007-12-21 Thread bernardo lagos alvarez
useR's I need transform the matrix wdat [,1] [,2] [,3] [,4] [,5] [,6] [1,]010110 [2,]100110 [3,]000011 [4,]110010 [5,]111100 [6,]001000 to gal ob

Re: [R] adding lines to a barchart

2007-12-21 Thread Spilak,Jacqueline [Edm]
Thanks Deepayan, this helped a lot and gave me exactly what I wanted however I want a few changes and am not sure how to do them. In each panel the bars for the years are the same color with the bottom axis (the x-axis) labelled with the years. I would like each year to have a specific color with

Re: [R] Finding overlaps in vector

2007-12-21 Thread Charles C. Berry
On Fri, 21 Dec 2007, Johannes Graumann wrote: > > > Dear all, > > I'm trying to solve the problem, of how to find clusters of values in a > vector that are closer than a given value. Illustrated this might look as > follows: > > vector <- c(0,0.45,1,2,3,3.25,3.33,3.75,4.1,5,6,6.45,7,7.1,8) > > Wh

Re: [R] R script to start session (without automatically finishing)

2007-12-21 Thread Prof Brian Ripley
So you have to write a script that will not terminate until you 'press a button or something like this'. That's easy to do, but as you haven't told us your OS. E.g. on Windows, call winDialogString at the end. On Fri, 21 Dec 2007, Philipp Fechteler wrote: > Hello R friends > > I am quite impr

[R] R appearance under linux

2007-12-21 Thread phguardiol
Dear R users, I have just moved to R2.6.1 under Opensuse linux 10.3. I used to work with R under XPpro. Is it "normal" to have a visual aspect of R under linux different ? I mean without a window with menus etc... What I do is just open a console window and type R then I just obtain a:?? >??

[R] R script to start session (without automatically finishing)

2007-12-21 Thread Philipp Fechteler
Hello R friends I am quite impressed by the power of R, I am using it only since some weeks now. But its visualizing capabilities are outstanding! But one thing I couldn't solve: I have programs producing lots of data, most times 3D. In R I am using the library rgl to visualize nicely the 3D d

Re: [R] install R on Unix without Admin Privilege

2007-12-21 Thread Gabor Csardi
On Fri, Dec 21, 2007 at 10:46:39AM -0500, Wensui Liu wrote: > HI, Gabor, > what program do i need to compile R on unix? At least a C compiler, maybe a Fortran compiler too. Some standard Unix tools like make as well. See also: http://cran.r-project.org/doc/manuals/R-admin.html > local dir? do you

Re: [R] Finding overlaps in vector

2007-12-21 Thread Gabor Grothendieck
This may not be as direct as Jim's in terms of specifying granularity but will uses conventional hierarchical clustering to create the clusters and also draws a nice dendrogram for you. I have split the dendrogram at a height of 0.5 to define the clusters but you can change that to whatever granu

[R] Solution: Available environment variables

2007-12-21 Thread Thompson, David (MNR)
Thank you Duncan (additional comments below), >-Original Message- >From: Duncan Murdoch [mailto:[EMAIL PROTECTED] >Sent: December 20, 2007 02:36 PM >To: Thompson, David (MNR) >Cc: r-help@r-project.org >Subject: Re: [R] Available environment variables > >On 20/12/2007 2:13 PM, Thompson, Da

Re: [R] Sorting factors

2007-12-21 Thread jim holtman
By default, factors are characters and sorted in alphabetical order. It looks like somehow you numeric data was converted to factors. Therefore '10' comes before '2'. If you want the factors in numeric order you have to convert them back. Look in the FAQs. > x <- c(1,2,3,4,5,10,11,20,21,22,30) >

[R] Sorry for "off-topic"

2007-12-21 Thread Michal Kneifl
Mary Christmas and happy new year to all contributors:-) Michael __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal

Re: [R] "gam()" in "gam" package

2007-12-21 Thread Thomas Lumley
On Thu, 20 Dec 2007, Kunio takezawa wrote: > R-users > E-mail: r-help@r-project.org > > I found the answer myself. > '.Fortran("baklo",' in lo.wam() and .Fortran("bakfit",in > s.wam() may carry out backfitting. But I cannot > create an R code which gives the same results as those of > "bakfit".

Re: [R] install R on Unix without Admin Privilege

2007-12-21 Thread Wensui Liu
HI, Gabor, what program do i need to compile R on unix? local dir? do you mean dir under my user name with my profile? On Dec 21, 2007 10:42 AM, Gabor Csardi <[EMAIL PROTECTED]> wrote: > Well i don't know which Unix you have, but it shouldn't matter anyway. > If you're able to compile programs on

[R] Sorting factors

2007-12-21 Thread tom sgouros
Hello all: I'm sorry to take up bandwidth with easy questions, but as an R beginner, I continue to be surprised by factors. If someone can shed some light for me, I'd be grateful: > d Total z02801 z02802 z02804 z02806 z02807 z02808 z02809 z02812 z02813 54813 29 51169 236810

Re: [R] install R on Unix without Admin Privilege

2007-12-21 Thread Gabor Csardi
Well i don't know which Unix you have, but it shouldn't matter anyway. If you're able to compile programs on the machine then download the R source, compile it and install it into your local directory. Gabor On Fri, Dec 21, 2007 at 10:25:36AM -0500, Wensui Liu wrote: > Good morning, Dear Lister,

Re: [R] (no subject)

2007-12-21 Thread John Kane
Hi, I don't think you have given enough information for anyone to offer advice. A small working (or almost working ) exmple that reproduces the problem would help. --- Jiaming Zuo <[EMAIL PROTECTED]> wrote: > Dear R Users, > > I am working for the United Nations to construct > a complet

Re: [R] Finding overlaps in vector

2007-12-21 Thread jim holtman
Here is a modification of the algorithm to use a specified value for the overlap: > vector <- c(0,0.45,1,2,3,3.25,3.33,3.75,4.1,5,6,6.45,7,7.1,8) > # following add 0.5 as the overlap detection -- can be changed > x <- rbind(cbind(value=vector, oper=1, id=seq_along(vector)), +cbind(valu

[R] install R on Unix without Admin Privilege

2007-12-21 Thread Wensui Liu
Good morning, Dear Lister, I really like to use R on our unix server. However, I am not an admin of the system. I am wondering if it is possible to install R on unix without admin previlege. Thank you so much! Have a nice Xmas season! __ R-help@r-projec

[R] Plotting multiple jpegs on a plot?

2007-12-21 Thread jchernev
Hello, I have a set of data which I currently represent graphically on a plot. To do that, I use the embedded functions from the graphics package (symbols, rect, etc.) However, I would like to use static jpeg images to represent the data if possible. Here's my progress so far: 1. Used the rimage

Re: [R] Efficient way to find consecutive integers in vector?

2007-12-21 Thread Martin Maechler
> "MS" == Marc Schwartz <[EMAIL PROTECTED]> > on Thu, 20 Dec 2007 16:33:54 -0600 writes: MS> On Thu, 2007-12-20 at 22:43 +0100, Johannes Graumann wrote: >> Hi all, >> >> Does anybody have a magic trick handy to isolate directly consecutive >> integers from somethin

Re: [R] NaN as a parameter in NLMINB optimization

2007-12-21 Thread Duncan Murdoch
On 21/12/2007 12:07 AM, Rebecca Sela wrote: > I am trying to optimize a likelihood function using NLMINB. After running > without a problem for quite a few iterations (enough that my intermediate > output extends further than I can scroll back), it tries a vector of > parameter values NaN. Thi

[R] RE : RE: using apply to loop [SEC=UNCLASSIFIED]

2007-12-21 Thread Louis Martin
Thanks Joe, It works well. How can I expand the output of z to include all other classes from 1 to nclass Right now I the output is: > z truth classified 4 5 8 10 1 1 1 1 0 4 6 10 8 0 10 1 0 1 1 > Is it possible to have something like: [,1]

[R] Within and Between matrix

2007-12-21 Thread cgenolin
Hi the list... I am working on cluster analysis, more precisely on Calinski and Harabasz criterion ( C(g)=Trace(between)/Trace(Within)*(n-g)/(g-1)) ) Is there a package for calculating within and between matrix? Morever, when we are using matrix with missing value, how do they handle it? On a mo

Re: [R] RMySQL installation problem - partially solved

2007-12-21 Thread Knut Krueger
First question: is this the right mailing list to discuss installation problems? if yes: (Windows XP professional environment ) Installing a second packege after the first causes an error: ---

Re: [R] array addition

2007-12-21 Thread Robin Hankin
[snip snip snip] > suppose I have two arrays x1,x2 of dimensions a1,b1,c1 and > a2,b2,c2 respectively. > > I want x = x1 "+" x2 with dimensions c(max(a1,a2), max(b1,b2),max > (c1,c2)) [snip snip snip] perhaps it wouldn't be too much to ask for you to check the most recent version of the "

[R] Odp: creating a factor from dates by subject?

2007-12-21 Thread Petr PIKAL
[EMAIL PROTECTED] napsal dne 20.12.2007 16:33:05: > Dear R-help, > > I have a data set consisting of measurements made on multiple > subjects. Measurement sessions are repeated for each subject on > multiple dates. Not all subjects have the same number of > sessions. To create a factor that re

Re: [R] Multicore computation in Windows network: How to set up Rmpi

2007-12-21 Thread Samu Mäntyniemi
Some progress in my problem: Samu Mäntyniemi kirjoitti: > With MPICH2 I managed to connect my computers so that I was able to > remotely launch Rgui on both machines but R hanged when calling > "library(Rmpi)". If only one Rgui was launched on the localhost, > "library(Rmpi)" worked without error

[R] Finding overlaps in vector

2007-12-21 Thread Johannes Graumann
Dear all, I'm trying to solve the problem, of how to find clusters of values in a vector that are closer than a given value. Illustrated this might look as follows: vector <- c(0,0.45,1,2,3,3.25,3.33,3.75,4.1,5,6,6.45,7,7.1,8) When using '0.5' as the proximity requirement, the following groups