Re: [R] Vector Manipulation

2006-06-19 Thread Petr Pikal
Hi or slightly quicker y-(x != N)*1 HTH Petr On 16 Jun 2006 at 22:57, Dimitrios Rizopoulos wrote: Date sent: Fri, 16 Jun 2006 22:57:48 +0200 From: Dimitrios Rizopoulos [EMAIL PROTECTED] To: Davis, Jacob B. [EMAIL PROTECTED] Copies to:

Re: [R] more  of a statistics/display question

2006-06-19 Thread Petr Pikal
Hi you could probably use image/contour/persp plots if your data are not completely chaotic. HTH Petr On 18 Jun 2006 at 17:56, [EMAIL PROTECTED] wrote: Date sent: Sun, 18 Jun 2006 17:56:35 -0500 (CDT) From: [EMAIL PROTECTED] To: [EMAIL

Re: [R] How to change the margin widths in png-plots?

2006-06-19 Thread Atte Tenkanen
I found one solution. It may be unofficial but it works. I put the par(mar=c(5.1, 7.1, 4.1, 2.1)) command between png() and plot() -commands. Atte Hello, I have tried to change the margin widths so that mtext (here sd of consecutive pc intervals, look at the picture) and

[R] Border line width?

2006-06-19 Thread Atte Tenkanen
Is there some way to change the line widths of plot borders? I couldn't find any parameters for that purpose. Atte Tenkanen University of Turku, Finland __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Border line width?

2006-06-19 Thread Prof Brian Ripley
On Mon, 19 Jun 2006, Atte Tenkanen wrote: Is there some way to change the line widths of plot borders? I couldn't find any parameters for that purpose. Do you mean the (optional) box? Use bty=n and then box(lwd=). plot(1:10, bty=n) box(lwd=3) works for me. -- Brian D. Ripley,

Re: [R] How to change the margin widths in png-plots?

2006-06-19 Thread Petr Pikal
Hi Perfectly official solution. First open graphic device (ps, png, jpeg, ...) and then change parameters according to your wish. HTH Petr On 19 Jun 2006 at 10:15, Atte Tenkanen wrote: Date sent: Mon, 19 Jun 2006 10:15:28 +0300 From: Atte Tenkanen [EMAIL

[R] combining tables

2006-06-19 Thread Robin Hankin
Hi Suppose I have two tables of counts of different animals and I wish to pool them so the total of the sum is the sum of the total. Toy example follows (the real ones are ~10^3 species and ~10^6 individuals). x and y are zoos that have merged and I need a combined inventory including animals

Re: [R] Border line width?

2006-06-19 Thread Petr Pikal
Hi plot(1,1) box(, lwd=5) HTH Petr On 19 Jun 2006 at 10:22, Atte Tenkanen wrote: Date sent: Mon, 19 Jun 2006 10:22:03 +0300 From: Atte Tenkanen [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Priority: normal Subject:

Re: [R] Border line width?

2006-06-19 Thread Dimitris Rizopoulos
maybe you're looking for: plot(0:1, 0:1) box(lwd = 3) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015

Re: [R] Border line width?

2006-06-19 Thread Atte Tenkanen
TNX Brian! Atte - Original Message - From: Prof Brian Ripley [EMAIL PROTECTED] Date: Monday, June 19, 2006 10:33 am Subject: Re: [R] Border line width? On Mon, 19 Jun 2006, Atte Tenkanen wrote: Is there some way to change the line widths of plot borders? I couldn't find any

Re: [R] Fitting Distributions Directly From a Histogram

2006-06-19 Thread Prof Brian Ripley
On Sun, 18 Jun 2006, Spencer Graves wrote: Won't 'interval censoring using 'survreg' in the 'survival' package handle this? (http://finzi.psych.upenn.edu/R/library/survival/html/survreg.html) Yes, but only for distributions it supports (and I think even with a user-supplied

Re: [R] combining tables

2006-06-19 Thread Dimitris Rizopoulos
one approach is the following: x - as.table(c(cats=3,squid=7,pigs=2,dogs=1,slugs=0)) y - as.table(c(cats=4,dogs=5,slugs=3,crabs=0)) out - merge(as.data.frame(x), as.data.frame(y), by = Var1, all = TRUE) desired.output - rowSums(out[c(Freq.x, Freq.y)], na.rm = TRUE)

[R] Bayesian Networks with deal

2006-06-19 Thread Carsten Steinhoff
Hello, I want to use R to model a bayesian belief network of frequencies for system failiures in various departments of a company. For the nodes I want to use a poisson-distribution parameterized with expert knowledge (e.g. with a gamma prior). Then I want to fill in learning-data to improve

Re: [R] Question concerning mle

2006-06-19 Thread Rainer M Krug
Thanks a lot Spencer for your tips - I'll look into all of them. Rainer Spencer Graves wrote: It's difficult to say much at this level of generality, but I have four suggestions: 1. Have you tried creating a reasonable grid of starting values using expand.grid and then

[R] Generating Random Numbers after Momentfitting

2006-06-19 Thread Jole, M. van
L.S. After moment fitting, which I programmed as follows: #Function which, given the first two moment, computates the density, #according to momentfitting momentfit-function(ES,c){ res-(1:5) if(c=1){ k-floor(1/c^2)

[R] Foometrics in R

2006-06-19 Thread Jan de Leeuw
One of the outcomes of useR! 2006 is that JSS is planning to publish a series of special volumes. They will be guest edited (I have guest editors already, although somewhat tentative in some cases). Each volume will have 5-10 issues (articles) of the usual JSS format. Psychometrics in R

[R] can I call user-created functions without source() ?

2006-06-19 Thread Rob Campbell
Hi, I have to R fairly recently from Matlab, where I have been used to organising my own custom functions into directories which are adding to the Matlab search path. This enables me to call them as I would one of Matlab's built-in functions. I have not found a way of doing the same thing in R. I

[R] function call

2006-06-19 Thread vincent
Dear R-users, When a function f1() is called, is there a way to know, from inside the function f1(), if f1() is directly called from the R console, or indirectly from another function f2() ? Of course, I may add an argument to f1(..., callbyf2=FALSE) only used by f2() giving explicitely this

Re: [R] function call

2006-06-19 Thread Peter Dalgaard
[EMAIL PROTECTED] writes: Dear R-users, When a function f1() is called, is there a way to know, from inside the function f1(), if f1() is directly called from the R console, or indirectly from another function f2() ? Of course, I may add an argument to f1(..., callbyf2=FALSE) only used

[R] how to do this sum?

2006-06-19 Thread GuangXing
Hi, Everybody! I have a big table which named table_x, and all the elements in the table is very large! Now I want to do the summay on the talbe_x[,3]. Unfornately, I can't get the right result! And the R give the warning messages as follow: sum(table_x[,3]) [1] NA Warning message: interger

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Ted Harding
On 19-Jun-06 Rob Campbell wrote: Hi, I have to R fairly recently from Matlab, where I have been used to organising my own custom functions into directories which are adding to the Matlab search path. This enables me to call them as I would one of Matlab's built-in functions. I have not

Re: [R] function call

2006-06-19 Thread vincent
Peter Dalgaard a écrit : [EMAIL PROTECTED] writes: Dear R-users, When a function f1() is called, is there a way to know, from inside the function f1(), if f1() is directly called from the R console, or indirectly from another function f2() ? Of course, I may add an argument to f1(...,

Re: [R] how to do this sum?

2006-06-19 Thread GuangXing
I am sorry for my ignorance. table_x is defined as follow: table_x-read.table(stat.txt) in the stat.txt,there are 5 cols,there are so many big intergers in the every cols. Now, I want to to that: sum(table_x[,3]) But it does not work. How can I? Pls help me! --

[R] useR! Thanks

2006-06-19 Thread Frank E Harrell Jr
After attending my first useR! conference I want to thank the organizers for doing a wonderful job and the presenters for their high quality presentations and stimulating ideas. The conference venue was excellent and of course Vienna is one of the greatest cities in the world to visit.

[R] MLE maximum number of parameters

2006-06-19 Thread Federico Calboli
Hi All, I would like to know, is there a *ballpark* figure for how many parameters the minimisation routines can cope with? I'm asking because I was asked if I knew. Cheers, Federico -- Federico C. F. Calboli Department of Epidemiology and Public Health Imperial College, St. Mary's Campus

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Duncan Murdoch
On 6/19/2006 5:36 AM, Rob Campbell wrote: Hi, I have to R fairly recently from Matlab, where I have been used to organising my own custom functions into directories which are adding to the Matlab search path. This enables me to call them as I would one of Matlab's built-in functions. I have

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Duncan Murdoch
On 6/19/2006 6:25 AM, (Ted Harding) wrote: On 19-Jun-06 Rob Campbell wrote: Hi, I have to R fairly recently from Matlab, where I have been used to organising my own custom functions into directories which are adding to the Matlab search path. This enables me to call them as I would one of

[R] Qurey : How to add trendline( st. line) in Graph

2006-06-19 Thread priti desai
How to add trendline (i.e. straight line passing through maximum points) in graph. I have worked on the data given below. Please tell me how to add trendline in the graph. The script is as follows === start #

Re: [R] Qurey : How to add trendline( st. line) in Graph

2006-06-19 Thread Petr Pikal
Hi maybe abline(0,1) HTH Petr On 19 Jun 2006 at 16:22, priti desai wrote: Date sent: Mon, 19 Jun 2006 16:22:24 +0530 From: priti desai [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject:[R] Qurey : How to add trendline( st.

Re: [R] how to do this sum?

2006-06-19 Thread Petr Pikal
Hi sum(as.numeric(table_x[,3]), na.rm=T) shall help (maybe) HTH Petr On 19 Jun 2006 at 19:05, GuangXing wrote: Date sent: Mon, 19 Jun 2006 19:05:51 +0800 From: GuangXing [EMAIL PROTECTED] To: Dr L. Y Hin [EMAIL PROTECTED] Copies to:

Re: [R] combining tables

2006-06-19 Thread Gabor Grothendieck
Try this: both - c(x,y) as.table(tapply(both, names(both), sum)) On 6/19/06, Robin Hankin [EMAIL PROTECTED] wrote: Hi Suppose I have two tables of counts of different animals and I wish to pool them so the total of the sum is the sum of the total. Toy example follows (the real ones are

Re: [R] function call

2006-06-19 Thread Gabor Grothendieck
Perhaps you what you want to do is to return an object that has a print method like this: f1 - function(x) structure(c(x, x^2), class = f1) print.f1 - function(x) cat(x is, x[1], x squared is, x[2], \n) Now you could do this: f1(3) # gives fancy output y - f1(3) # returns structure unclass(y)

Re: [R] Qurey : How to add trendline( st. line) in Graph

2006-06-19 Thread Gabor Grothendieck
Try: library(quantreg) abline(rq(g ~ f, tau = .999)) On 6/19/06, priti desai [EMAIL PROTECTED] wrote: How to add trendline (i.e. straight line passing through maximum points) in graph. I have worked on the data given below. Please tell me how to add trendline in the graph. The script is as

[R] saving rounded numbers as a new variable in a dataframe

2006-06-19 Thread Bob Green
A basic question, but one that eludes me. I have created a new variable $numurder, which I have rounded off. I want to save the rounded off version of this variable to an existing datafile called 'ngri.csv' . numurder -c((murder*no.of.cases)/100) [[1]] [1] 48.952 112.073 182.160

Re: [R] Qurey : How to add trendline( st. line) in Graph

2006-06-19 Thread Sean O'Riordain
lines(predict(lm(g~f))~f,lty=2) or if its a curvy line you're looking for... lines(predict(loess(g~f))~f,lty=3) On 19/06/06, priti desai [EMAIL PROTECTED] wrote: How to add trendline (i.e. straight line passing through maximum points) in graph. I have worked on the data given below. Please

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Joerg van den Hoff
Duncan Murdoch wrote: On 6/19/2006 6:25 AM, (Ted Harding) wrote: On 19-Jun-06 Rob Campbell wrote: Hi, I have to R fairly recently from Matlab, where I have been used to organising my own custom functions into directories which are adding to the Matlab search path. This enables me to call

Re: [R] saving rounded numbers as a new variable in a dataframe

2006-06-19 Thread Neuro LeSuperHéros
Hi 1-Please read this: R Data Import/Export cran.r-project.org/doc/manuals/R-data.pdf 2-This is the code: rounded -round(numurder, 0) write.csv(rounded,c:/ngri.csv) Neuro From: Bob Green [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject: [R] saving rounded numbers as a new variable in a

Re: [R] Qurey : How to add trendline( st. line) in Graph

2006-06-19 Thread Gabor Grothendieck
One other item. You had asked for the trendline passing through the maximum points, which is what my prior answer gives, but others interpreted your question as just as ordinary trendline so if that is what you are looking for try this: abline(lm(g ~ f)) On 6/19/06, Gabor Grothendieck [EMAIL

Re: [R] MLE maximum number of parameters

2006-06-19 Thread Ben Bolker
Federico Calboli f.calboli at imperial.ac.uk writes: Hi All, I would like to know, is there a *ballpark* figure for how many parameters the minimisation routines can cope with? I think I would make a distinction between theoretical and practical limits. A lot depends on how fast

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Duncan Murdoch
Just a few comments below on alternative ways to do the same things: On 6/19/2006 8:19 AM, Joerg van den Hoff wrote: for short term usage of some specialized functions I have added some lines to the `.Rprofile' in my home(!) directory as follows (probably there are smarter solutions, but at

Re: [R] lmer and mixed effects logistic regression

2006-06-19 Thread Rick Bilonick
On Sat, 2006-06-17 at 09:46 -0700, Spencer Graves wrote: 'lmer' RETURNS AN ERROR WHEN SAS NLMIXED RETURNS AN ANSWER Like you, I would expect lmer to return an answer when SAS NLMIXED does, and I'm concerned that it returns an error message instead. Your example is not

Re: [R] saving rounded numbers as a new variable in a dataframe

2006-06-19 Thread Petr Pikal
Hi On 19 Jun 2006 at 8:20, Neuro LeSuperHéros wrote: From: Neuro LeSuperHéros [EMAIL PROTECTED] To: [EMAIL PROTECTED], r-help@stat.math.ethz.ch Date sent: Mon, 19 Jun 2006 08:20:16 -0400 Subject:Re: [R] saving rounded numbers

[R] Function hints

2006-06-19 Thread hadley wickham
One of the recurring themes in the recent UserR conference was that many people find it difficult to find the functions they need for a particular task. Sandy Weisberg suggested a small idea he would like to see: a hints function that given an object, lists likely operations. I've done my best

Re: [R] MLE maximum number of parameters

2006-06-19 Thread Roger D. Peng
It really depends on how well-behaved your objective function is, but I've been able to fit a few models with 10--15 parameters. But I felt like I was stretching the limit there. -roger Federico Calboli wrote: Hi All, I would like to know, is there a *ballpark* figure for how many

[R] multivariate splits

2006-06-19 Thread Bálint Czúcz
Dear R users! Does someone know about any algorithms / packages in R, that perform classification / regression / decision trees using multivariate splits? I have done some research, but I found nothing. Packages tree and rpart seem only to be able to do CART with univariate splits. Thank you

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Joerg van den Hoff
Duncan Murdoch wrote: Just a few comments below on alternative ways to do the same things: On 6/19/2006 8:19 AM, Joerg van den Hoff wrote: for short term usage of some specialized functions I have added some lines to the `.Rprofile' in my home(!) directory as follows (probably there are

[R] help on nlm (gradient) (fwd)

2006-06-19 Thread Luz Maria Palacios Derflingher
Hello No worries anymore. Figured it out. Thanks for everything. Luz -- Forwarded message -- Date: Sun, 18 Jun 2006 23:12:25 -0600 (MDT) From: Luz Maria Palacios Derflingher [EMAIL PROTECTED] To: [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Subject: help on nlm (gradient) Hello. I

Re: [R] lmer and mixed effects logistic regression

2006-06-19 Thread Rick Bilonick
On Sun, 2006-06-18 at 13:58 +0200, Douglas Bates wrote: If I understand correctly Rick it trying to fit a model with random effects on a binary response when there are either 1 or 2 observations per group. I think that is very optimistic because there is so little information available per

Re: [R] function call

2006-06-19 Thread vincent
Gabor Grothendieck a écrit : Perhaps you what you want to do is to return an object that has a print method like this: f1 - function(x) structure(c(x, x^2), class = f1) print.f1 - function(x) cat(x is, x[1], x squared is, x[2], \n) Thank you Gabor for this idea. For the moment interactive()

[R] eacf

2006-06-19 Thread Michel Helcias
I am Brazilian and I don't know how to speak English, for that I apologize for my writing. I'd like of informations about the extended (sample) autocorrelation function (*EACF*) for the time series. where to acquire some related command. thak you. [[alternative HTML version deleted]]

Re: [R] Fitting Distributions Directly From a Histogram

2006-06-19 Thread Myers, Brent
I'm desperately uninformed about the methods you have been discussing, but for a non-parametric solution, how about splines? Carl DeBoor (1978) outlines a method of area preserving quadratic splines to fit histograms. All you need to know is the value each bar represents and its interval. I am

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Duncan Murdoch
On 6/19/2006 10:19 AM, Joerg van den Hoff wrote: Duncan Murdoch wrote: Just a few comments below on alternative ways to do the same things: On 6/19/2006 8:19 AM, Joerg van den Hoff wrote: for short term usage of some specialized functions I have added some lines to the `.Rprofile' in my

Re: [R] multivariate splits

2006-06-19 Thread Torsten Hothorn
On Mon, 19 Jun 2006, Bálint Czúcz wrote: Dear R users! Does someone know about any algorithms / packages in R, that perform classification / regression / decision trees using multivariate splits? I have done some research, but I found nothing. Packages tree and rpart seem only to be able to

[R] Graph and function

2006-06-19 Thread Ahmed Elhabti
I have a program which make calculate them of 4 estimator, B, BE, wB and wBE. Now I want to varied n example: n-seq(1,100,1), and to make a graph of B and Be according to n, to make a comparison. Tank you very much. n-50 a-0.2 theta-rexp(n,a)

Re: [R] nested mixed-effect model: variance components

2006-06-19 Thread Christoph Buser
Dear John I've put your mail back to the R list since I have no explanation for the lmer result. Maybe someone else has an idea. I adapted it to show some else that I do not understand. ## Creation of the data (habitat is nested in lagoon): set.seed(1) dat - data.frame(y = rnorm(100), lagoon =

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Joerg van den Hoff
Duncan Murdoch wrote: On 6/19/2006 10:19 AM, Joerg van den Hoff wrote: Duncan Murdoch wrote: Just a few comments below on alternative ways to do the same things: On 6/19/2006 8:19 AM, Joerg van den Hoff wrote: for short term usage of some specialized functions I have added some lines to

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Michael H. Prager
Rob, I accomplish what you ask by putting my functions into a dedicated directory, starting R there and sourcing them, and loading the resulting workspace in .Rprofile with this: attach(d:/R/MHP/MHPmisc/.RData) I find that procedure simpler than learning the package mechanism. It is easy to

Re: [R] function call

2006-06-19 Thread Martin Maechler
vincent == vincent [EMAIL PROTECTED] on Mon, 19 Jun 2006 16:30:54 +0200 writes: vincent Gabor Grothendieck a écrit : Perhaps you what you want to do is to return an object that has a print method like this: f1 - function(x) structure(c(x, x^2), class = f1)

Re: [R] MLE maximum number of parameters

2006-06-19 Thread Spencer Graves
Applications with lots of parameters also tend to have parameters in a relatively small number of families, and each of these few families could be considered to have a distribution. Splines, for example, have lots of parameters -- sometimes more parameters than observations (as do

Re: [R] Graph and function

2006-06-19 Thread Petr Pikal
Hi well I am a bit puzzled a-0.2 t-rexp(1,a) xn-rpois(1,t) B-(xn+1)/(1+a) so B has nothing to do with n b-mean(x) BE-(1+xn)*b/(1+b) so as Be you can use n - 1:100 and n as an input vector for some vectorized function but I wonder where to put such vector in your pieces of

Re: [R] Function hints

2006-06-19 Thread Joerg van den Hoff
hadley wickham wrote: One of the recurring themes in the recent UserR conference was that many people find it difficult to find the functions they need for a particular task. Sandy Weisberg suggested a small idea he would like to see: a hints function that given an object, lists likely

Re: [R] multivariate splits

2006-06-19 Thread Liaw, Andy
If by multivariate split the OP meant splitting on combinations of covariates (instead of multi-way split on a single covariate at a time), there aren't that many methods published (AFAIK). All the ones I know about are in the cheminformatics area: RP-SA and generalizations/extensions of it. No

Re: [R] Problem in Projection pursuit regression . [Broadcast]

2006-06-19 Thread Liaw, Andy
ppr() estimates the ridge functions by either the Friedman super smoother or smoothing splines. (See ?supsmu and ?smooth.spline.spline.) Andy From: anil kumar rohilla Hi list, Is any body has implemented the projection pursuit regression (ppr) method in MASS library.Running

[R] analyze amino acid sequence (composition)of proteins

2006-06-19 Thread stubben
Have you checked the package seqinR at CRAN? It may help. Also, if you split a string into single characters using strsplit, then use table to count characters. seq-ATGAAC table(strsplit(seq, )) A C G T 3 1 1 1 3. based on AAABBB,how can i get some statistics of this string such as

[R] objects

2006-06-19 Thread H. Paul Benton
Sorry I'm still learning R so here my question. I have a subroutine that I have made. I want to be able to return two items from it. One is an object (the nls object) the other is matrix/data.frame. I would like to pass both of these to another subroutine and also be able to use it outside as

[R] Display Conditional Probabilities on y-label

2006-06-19 Thread Brian Phillip Weaver
Hello, I am trying to display a conditional probability for the y-label on a beta distribution plot. Here is the code I have been using: z = expression(f(x|) g = paste(z,expression(alpha),,,expression(beta),)) plot(x,y,ylim=c(0,3),xlab=x,ylab=g,lwd=3,type=l,cex.lab=1.3) The output on the plot

[R] Compute in C a loop with the lm() function and back to R

2006-06-19 Thread Javier López de Lacalle Beltrán de Heredia
I want to do the computations from a loop by means of a C routine that will be called from R. The lm() function is entailed in that loop. As a guidance, I have taken the following example from http://www.math.mcgill.ca/steele/ (As far as my question is concerned, I get similar information from

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Duncan Murdoch
On 6/19/2006 11:38 AM, Michael H. Prager wrote: Rob, I accomplish what you ask by putting my functions into a dedicated directory, starting R there and sourcing them, and loading the resulting workspace in .Rprofile with this: attach(d:/R/MHP/MHPmisc/.RData) This seems like a nice

[R] Plotting Upper triangle of Matrix with diagonal as the Base

2006-06-19 Thread Cal Stats
Hi.. I a have a symmetric matrix to plot . I would like to plot only the Upper triangle but with the diagonal as the Base of the rectangle. Is there an easy way to do it. Thanks. Harsh - [[alternative HTML

Re: [R] Display Conditional Probabilities on y-label

2006-06-19 Thread Charles Annis, P.E.
The trick is not to paste expressions but to make an expression of the pastes: composite.expression - expression(paste(f(x | , alpha, , , beta, ))) dev.off() par(mar=c(5,5,1,1)+0.1) plot(NA,xlim=c(0,1), ylim=c(0,3),xlab=x,ylab=composite.expression,lwd=3,type=l,cex.lab=1.3) Charles Annis,

Re: [R] Function hints

2006-06-19 Thread Duncan Murdoch
On 6/19/2006 12:14 PM, Joerg van den Hoff wrote: just a feedback: that's a useful function, thank you. but the problem is probably more general: frequently I do not really want to know what I generally can do with a data frame, for instance, but rather I would like to use `help.search' as

[R] frechet distance

2006-06-19 Thread Rajarshi Guha
Hi, is there any package (or source code snippet) that will evaluate the Frechet distance for curves represented as sets of points? Searching around only threw up references to a Frechet distribution. Thanks, --- Rajarshi Guha

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Spencer Graves
I'd like to echo Duncan's comment about creating 'man' pages for functions. I find myself writing man pages for functions I may never include in an R package, just because it helps me think through what I'm trying to do. I think I get better code faster. It's like preparing a map

Re: [R] Function hints

2006-06-19 Thread Jonathan Baron
On 06/19/06 13:13, Duncan Murdoch wrote: `help.search' does not allow full text search in the manpages (I can imagine why (1000 hits...), but without such a thing google, for instance, would probably not be half as useful as it is, right?) and there is no sorting by relevance in the

Re: [R] Foometrics in R

2006-06-19 Thread Albyn Jones
On Mon, Jun 19, 2006 at 02:34:17AM -0700, Jan de Leeuw wrote: One of the outcomes of useR! 2006 is that JSS is planning to publish a series of special volumes. They will be guest edited (I have guest editors already, although somewhat tentative in some cases). Each volume will have 5-10

[R] lattice xyplot - aligning date labels so that they align with the grid lines in panel.grid

2006-06-19 Thread john.gavin
Hi, I have a basic question about aligning date labels for the x-axis in an xyplot so that they align with the grid lines from the panel.grid argument. For example, with x - data.frame( date = seq(as.Date(2005/01/01), as.Date(2006/06/01), length.out = 20), value = runif(20)) xyplot(value ~

Re: [R] MLE maximum number of parameters

2006-06-19 Thread Patrick Burns
Seagulls have a very different perspective to ballparks than ants. Nonetheless, there is something that can be said. There are several variables in addition to the number of parameters that are important. These include: * The complexity of the likelihood * The number of observations in the

[R] strange digit switching

2006-06-19 Thread Hans-Peter
Hi, I execute the following code snippet.What I don't understand is, that in the first two cases the numbers shown don't correnspondant to the input (though I have set the options to show 22 digits). When I assign the third number which has one place, the display of the numbers is suddenly

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Sean O'Riordain
Indeed, some folk say that the documentation should be written before the code... Sean On 19/06/06, Duncan Murdoch [EMAIL PROTECTED] wrote: 2. If you use the package system, you will be encouraged to create man pages for your functions. This is work, but I think it pays off in the end. I

Re: [R] strange digit switching

2006-06-19 Thread Peter Dalgaard
Hans-Peter [EMAIL PROTECTED] writes: Hi, I execute the following code snippet.What I don't understand is, that in the first two cases the numbers shown don't correnspondant to the input (though I have set the options to show 22 digits). When I assign the third number which has one place,

Re: [R] strange digit switching

2006-06-19 Thread Thomas Lumley
On Mon, 19 Jun 2006, Hans-Peter wrote: Hi, I execute the following code snippet.What I don't understand is, that in the first two cases the numbers shown don't correnspondant to the input (though I have set the options to show 22 digits). When I assign the third number which has one place,

Re: [R] lmer binomial model overestimating data?

2006-06-19 Thread Thomas Lumley
On Wed, 14 Jun 2006, Martin Henry H. Stevens wrote: Hi folks, Warning: I don't know if the result I am getting makes sense, so this may be a statistics question. The fitted values from my binomial lmer mixed model seem to consistently overestimate the cell means, and I don't know why. I

Re: [R] MLE maximum number of parameters

2006-06-19 Thread Albyn Jones
I regularly optimize functions of over 1000 parameters for posterior mode computations using a variant of newton-raphson. I have some favorable conditions: the prior is pretty good, the posterior is smooth, and I can compute the gradient and hessian. albyn On Mon, Jun 19, 2006 at 06:53:00PM

Re: [R] number of iteration s exceeded maximum of 50

2006-06-19 Thread Leaf Sun
Thanks to Douglas and all others who responded. I applied nls(y ~ a*x^b, start = list(a = a1, b = b1), control = list(maxiter = 500), trace=TRUE) to increase the number of iterations, found it successful. The suggestion Douglas raised in plotting the data and then tracing the

Re: [R] lattice xyplot - aligning date labels so that they align with the grid lines in panel.grid

2006-06-19 Thread Deepayan Sarkar
On 6/19/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Hi, I have a basic question about aligning date labels for the x-axis in an xyplot so that they align with the grid lines from the panel.grid argument. For example, with x - data.frame( date = seq(as.Date(2005/01/01),

Re: [R] analyze amino acid sequence (composition)of proteins

2006-06-19 Thread Jean lobry
i am a newer to R and i am doing some protein category classification based on the amino acid sequence.while i have some questions urgently. I'm not sure to understand exactly what you are trying to do, but if this is a kind of clustering from the amino-acid composition of proteins beware that

Re: [R] multivariate splits

2006-06-19 Thread Bálint Czúcz
Yes, I am interested in tree-building algorithms that consider some (linear) combination of the covariates during the splitting. Sorry if I have not been clear enough. There seem to be a large variety of methods (e.g. see http://www.tigr.org/~salzberg/murthy_thesis/survey/node11.html), I hoped

Re: [R] can I call user-created functions without source() ?

2006-06-19 Thread Michael H. Prager
on 6/19/2006 1:32 PM Spencer Graves said the following: I'd like to echo Duncan's comment about creating 'man' pages for functions. I find myself writing man pages for functions I may never include in an R package, just because it helps me think through what I'm trying to do. I

[R] Nested variance-covariance matrix in Multilevel model

2006-06-19 Thread Tobias Guennel
Dear R community, I have trouble implementing a nested variance-covariance matrix in the lme function. The model has two fixed effects called End and logpgc, the response variable is the logarithm to base 2 of Intensity ( log2(Intensity) ) and the random effects are called Probe and ProbeNo.

[R] Nested variance-covariance matrix in Multilevel model

2006-06-19 Thread Tobias Guennel
I completely forgot to supply the R code I tried: vov1i2-read.table(VOV1_INHIBITED6-16-2006-13h35min33sec.txt,header=TRUE) test.lme-lme(fixed=log2(Intensity)~End+logpgc,random=list(Probe=pdBlocked(list(~1,pdCompSymm(~1 ),End=~1,ProbeNo=pdCompSymm(~1)),data=vov1i2) It doesn't look right to me

Re: [R] Change the range of a 'ternaryplot'

2006-06-19 Thread Greg Snow
Try the function triangle.plot from the ade4 package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Martínez Ovando Juan Carlos

[R] binary for Rmath standalone library for Visual C++ 8.0?

2006-06-19 Thread Faheem Mitha
Hi, Does anyone happen to have compiled the Rmath standalone library for Visual C++ 8.0, and can offer it for me to download? I'm not that familiar with Microsoft or Visual C++. I mostly use Linux. However, my boss needs to use this library on Windows to compile some code I wrote, and I

Re: [R] binary for Rmath standalone library for Visual C++ 8.0?

2006-06-19 Thread Faheem Mitha
On Mon, 19 Jun 2006, Faheem Mitha wrote: Hi, Does anyone happen to have compiled the Rmath standalone library for Visual C++ 8.0, and can offer it for me to download? I'm not that familiar with Microsoft or Visual C++. I mostly use Linux. However, my boss needs to use this library on

Re: [R] Post Stratification

2006-06-19 Thread Thomas Lumley
On Sun, 18 Jun 2006, Mark Hempelmann wrote: Dear WizaRds, having met some of you in person in Vienna, I think even more fondly of this community and hope to continue on this route. It was great talking with you and learning from you. Thank you. I am trying to work through an

Re: [R] function call

2006-06-19 Thread Gabor Grothendieck
On 6/19/06, Martin Maechler [EMAIL PROTECTED] wrote: vincent == vincent [EMAIL PROTECTED] on Mon, 19 Jun 2006 16:30:54 +0200 writes: vincent Gabor Grothendieck a écrit : Perhaps you what you want to do is to return an object that has a print method like this: f1 -

Re: [R] Custom Command to Generate SQL

2006-06-19 Thread Alex Restrepo
Hi Gabor: Thanks for the great example. I am an R newbie, so please forgive my question, but could you describe what the sub() function is doing in your example? Why is there an and in the first argument to the sub() function? Many Thanks: Alex From: Gabor Grothendieck [EMAIL PROTECTED]

[R] how to put the results of loop into a dataframe

2006-06-19 Thread zhijie zhang
Dear friends, suppose i want to do the following caulation for 100 times, how to put the results of x , y and z into the same dataframe/dataset? x-runif(1) y-x+1 z-x+y thanks in advance! -- Kind Regards, Zhi Jie,Zhang ,PHD Department of Epidemiology School of Public Health Fudan University

[R] how to put the results of loop into a dataframe

2006-06-19 Thread zhijie zhang
Dear friends, suppose i want to do the following caulation for 100 times, how to put the results of x , y and z into the same dataframe/dataset? x-runif(1) y-x+1 z-x+y thanks in advance! -- Kind Regards, Zhi Jie,Zhang ,PHD Department of Epidemiology School of Public Health Fudan University

Re: [R] Custom Command to Generate SQL

2006-06-19 Thread Gabor Grothendieck
Hi, Here is an example of paste: paste(and, 1:3, collapse = ) [1] and 1 and 2 and 3 in which and was prepended to each of the vector entries and then everything was collapsed together. In our situation we don't actually want the first 'and' so we substitute for it, using sub, the first

Re: [R] Effect size in mixed models

2006-06-19 Thread Spencer Graves
You have asked a great question: It would indeed be useful to compare the relative magnitude of fixed and random effects, e.g. to prioritize efforts to better understand and possibly manage processes being studied. I will offer some thoughts on this, and I hope if there are errors

  1   2   >