Re: [R] Screen settings for point of view in lattice and misc3d

2010-03-05 Thread Deepayan Sarkar
On Wed, Mar 3, 2010 at 1:22 PM, Waichler, Scott R wrote: > I'm making some 3D plots with contour3d from misc3d and wireframe from > lattice.  I want to view them from below; i.e. the negative z-axis.  I can't > figure out how to do so.  I would like my point of view looking up from > below, wit

Re: [R] How to create a line and bar panel chart with two different axes?

2010-03-05 Thread Deepayan Sarkar
On Wed, Mar 3, 2010 at 12:35 PM, DougNiu wrote: > > I need to create a line and bar panel chart with two different axes. I tried > in lattice but couldn't get it worked. Here is my code: > > data(barley) > barchart(yield ~ variety | site, data = barley, >              groups = year, layout = c(1,6

Re: [R] Robust SE for lrm object

2010-03-05 Thread Achim Zeileis
On Sat, 6 Mar 2010, David Winsemius wrote: On Mar 5, 2010, at 11:54 PM, Patrick Shea wrote: I'm trying to obtain the robust standard errors for a multinomial ordered logit model: mod6 <- lrm(wdlshea ~ initdesch + concap + capasst + qualrat + terrain,data=full2) The model is fine but w

Re: [R] scientific (statistical) foundation for Y-RANDOMIZATION in regression analysis

2010-03-05 Thread Greg Snow
In the stats literature these are more often called permutation tests. Looking up that term should give you some results (if not, I have some references, but they are at work and I am not, I could probably get them for you on Monday if you have not found anything before then). -- Gregory (Gre

[R] complex--scatter plot

2010-03-05 Thread OEM Configuration (temporary user)
My input "ID" "Label" "*Stype*" "Ntype" "Stype_No""*log*" "S1" "xxx" "A/A" 1 6 2.8 "S1" "xxx" "A/G" 2 2 3 "S1" "xxx" "G/G" 3 1 4 "S2" "yyy" "A/A" 1 1 6.8 "S2" "yyy" "A/G" 2 2 7 "S2" "yyy" "G/G" 3 6

Re: [R] conditioning variable in panel.xyplot?

2010-03-05 Thread Deepayan Sarkar
On Fri, Mar 5, 2010 at 12:45 PM, Seth W Bigelow wrote: > I'm stumped after an hour or so reading about subscripts in panel.xyplot. > Apparently the panel function is executed for each subset of data in the > main dataset (specified by the conditioning variable, 'site' in my > example), and the 'su

Re: [R] Robust SE for lrm object

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 11:54 PM, Patrick Shea wrote: I'm trying to obtain the robust standard errors for a multinomial ordered logit model: mod6 <- lrm(wdlshea ~ initdesch + concap + capasst + qualrat + terrain,data=full2) The model is fine but when I try to get the RSE I get an error.

[R] Robust SE for lrm object

2010-03-05 Thread Patrick Shea
I'm trying to obtain the robust standard errors for a multinomial ordered logit model: mod6 <- lrm(wdlshea ~ initdesch + concap + capasst + qualrat + terrain,data=full2) The model is fine but when I try to get the RSE I get an error. coeftest(mod6, vcov = vcovHAC(mod6)) Error in match.arg(

Re: [R] converting multiple lines of text to a data frame

2010-03-05 Thread Phil Spector
Andrew- Maybe something like this: dd = read.table(filename) unstack(dd,V2~V1) A. B. C. 1 1 2 10.0 2 34 20 6.7 3 2 78 35.0 - Phil Spector Statistical Computing Facility

Re: [R] Writing own simulation function in C

2010-03-05 Thread Sharpie
TheSavageSam wrote: > > I am wishing to write my own random distribution simulation function using > C programmin language(for speed) via R. I am familiar with R programming > but somewhat new to C programming. I was trying to understand "Writing R > extensions" -guide and its part 6.16, but I f

Re: [R] transposing data

2010-03-05 Thread Phil Spector
Frank - I think you need to create a composite time variable to do what you want to do: exp1.r5$key = with(exp1.r5,paste(CannonAngle,CannonOriB, CannonOriR,nRedPelelts,TargetColor,tbearing,sep='.')) exp1.r5.use = subset(exp1.r5,select=-c(CannonAngle,CannonOriB,

[R] converting multiple lines of text to a data frame

2010-03-05 Thread Andrew Yee
I'm trying to find a way for converting multiple lines of text into a table.  I'm not sure if there's a way where you can use read.delim() to read in multiple lines of text and create the following data frame with something akin to rehape()?. Apologies if there is an obvious way to do this. A: 1

[R] transposing data

2010-03-05 Thread Frank Tamborello
Hi. I have repeated measures data of the form where each observation is a trial, and trials are grouped by subject, and variables encode whatever level of a factor was present during that trial, and the dependent variable is response time (RT). I want to transpose the data to a form suitabl

Re: [R] how to make this sequence: 1,2,3,4,5,4,3,2,1

2010-03-05 Thread kMan
c(x,(x<-1:5)[4:1]) -Original Message- From: baptiste auguie [mailto:baptiste.aug...@googlemail.com] Sent: Friday, March 05, 2010 1:08 AM To: kensuguro Cc: r-help@r-project.org Subject: Re: [R] how to make this sequence: 1,2,3,4,5,4,3,2,1 c(x <- 1:5, rev(x[-length(x)])) On 5 March 2010

[R] memory error in for loop

2010-03-05 Thread Peter
hi, I have been attempting to run this script and am getting some strange results. The script connects to a database and retrieves a series of tables, using sequential sql statements. I have tested all of the sql statements in the PostGreSQL terminal and they all return the desired results. I

Re: [R] Three most useful R package

2010-03-05 Thread kMan
(1) - nlme, lattice, stats (2) - a usable large-file/out of memory regression package that abstracts "all" the details of connections & etc from the user, accept perhaps the initial function call, so I don't have to actually know anything about the file I'm opening, how big it is, how many lines of

[R] Plot interaction in multilevel model

2010-03-05 Thread dadrivr
I am trying to plot an interaction in a multilevel model. Here is some sample data. In the following example, it is longitudinal (i.e., repeated measures), so the outcome, score (at each of the three time points), is nested within the individual. I am interested in the interaction between gende

Re: [R] Writing own simulation function in C

2010-03-05 Thread stephen sefick
Look in the sources of runif for the C code? On Fri, Mar 5, 2010 at 3:54 PM, TheSavageSam wrote: > > I am wishing to write my own random distribution simulation function using C > programmin language(for speed) via R. I am familiar with R programming but > somewhat new to C programming. I was try

[R] Plotting Comparisons with Missing Data

2010-03-05 Thread Alastair
Hi, I'm new to R and I've run into a problem that I'm not really sure how to express properly in the language. I've got a data table that I've read from a file containing some simple information about the performance of 4 algorithms. The columns are the name of the algorithm, the problem instanc

Re: [R] test question, variance

2010-03-05 Thread Phil Spector
Contact your instructor, and let them know that you are having problems with an assignment. Most instructors have office hours, where you can discuss problems like this. - Phil Spector Statistical Computing Facilit

[R] scientific (statistical) foundation for Y-RANDOMIZATION in regression analysis

2010-03-05 Thread Damjan Krstajic
Dear all, I am a statistician doing research in QSAR, building regression models where the dependent variable is a numerical expression of some chemical activity and input variables are chemical descriptors, e.g. molecular weight, number of carbon atoms, etc. I am building regression models a

[R] Bootstrap standard errors

2010-03-05 Thread Rman
I have this coding to help me work out the bootstrap standard errord for the following array/matrix dat <- matrix(c(8,23,14,13,9,11,25,25,14,11,4,25,19,7,13,1,8,16,2,6,3,4,6,5,7), nrow = 5, ncol=5, byrow=TRUE) the function is: mybootstrap <- function(dat, nr, nc, m) { # dat is a two dimentiona

[R] test question, variance

2010-03-05 Thread SquareAce
Can someone get me going in the right direction with this test question? We have went as far as t.test in our class, but not to power or anova yet. I have pasted the question and a summary of the dataset "gpa". 5. Given that all the GPA data in this dataset is comprised of means, would you expect

Re: [R] How to parse the arguments from a function call and evaluate them in a dataframe?

2010-03-05 Thread Thomas Lumley
On Fri, 5 Mar 2010, Ravi Varadhan wrote: Hi, I would like to write a function which has the following syntax: myfn <- function(formula, ftime, fstatus, data) { # step 1: obtain terms in `formula' from dataframe `data' # step 2: obtain ftime from `data' # step 3: obtain fstatus from `data' # s

Re: [R] How to assign week numbers to a time-series

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 4:46 PM, David Winsemius wrote: On Mar 5, 2010, at 4:23 PM, Hosack, Michael wrote: Hello everyone, My progress has stalled on finding a way of creating a somewhat complicated variable to add to my existing dataframe and I am hoping one of you could help me out. The da

[R] How to parse the arguments from a function call and evaluate them in a dataframe?

2010-03-05 Thread Ravi Varadhan
Hi, I would like to write a function which has the following syntax: myfn <- function(formula, ftime, fstatus, data) { # step 1: obtain terms in `formula' from dataframe `data' # step 2: obtain ftime from `data' # step 3: obtain fstatus from `data' # step 4: do model estimation # step 5: return

Re: [R] How to assign week numbers to a time-series

2010-03-05 Thread Achim Zeileis
On Fri, 5 Mar 2010, Hosack, Michael wrote: Hello everyone, My progress has stalled on finding a way of creating a somewhat complicated variable to add to my existing dataframe and I am hoping one of you could help me out. The dataframe below contains only a fraction of the data of my complete

[R] R algorithm for maximum curvatures measures of nonlinear models

2010-03-05 Thread Walmes Marques Zeviani
Hi all, I'm looking for the algorithm to calculate maximum intrinsic and parameter effects curvature of Bates & Watts (1980). I have Bates & Watts (1980) original article, Bates et al (1983) article, Seber & Wild (1989) and Ratkowsky (1983) books. Those sources show steps and algorithm to get

[R] Redhat Linux Install

2010-03-05 Thread Ryan Garner
I just installed R on Redhat Linux at work for the first time and have two questions. 1. I tried to install R to have png and cairo capabilities and was unsuccessful. Before running make, I ran ./configure --with-libpng=yes --with-x=no --with-cairo=yes --with-readline-yes . R installed fine, but

[R] How to assign week numbers to a time-series

2010-03-05 Thread Hosack, Michael
Hello everyone, My progress has stalled on finding a way of creating a somewhat complicated variable to add to my existing dataframe and I am hoping one of you could help me out. The dataframe below contains only a fraction of the data of my complete dataframe, but all of the variables. What I

[R] Writing own simulation function in C

2010-03-05 Thread TheSavageSam
I am wishing to write my own random distribution simulation function using C programmin language(for speed) via R. I am familiar with R programming but somewhat new to C programming. I was trying to understand "Writing R extensions" -guide and its part 6.16, but I found it hard to understand(http:

Re: [R] Sweave and optional document sections

2010-03-05 Thread Aleksey Naumov
Max, Thank you very much! I'll give LaTeX a go (haven't used it in a while). Is there a recommended distribution of LaTeX for Windows? Aleksey On Fri, Mar 5, 2010 at 1:50 PM, Max Kuhn wrote: > That isn't hard to do with Sweave. > > Within you if branch, you can have your code write out any mar

Re: [R] Update RMySQL and ... it's no more running

2010-03-05 Thread bernese
Yep this worked. I encountered this problem only after installing the ODBC connector plugin. Not sure if there's any correlation. I installed the MySQL 5.1.33 binary from http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL and that will fixed 'RMySQL was compiled with MySQL 5.0.67 but loading My

[R] convert rpart tree to phylo (in ape)

2010-03-05 Thread Chris Hane
Hello, Has anyone written a conversion from rpart tree class to phylo or hclust trees? The plot.phylo method allows much more readable layouts. Thanks! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.eth

Re: [R] how to make this sequence: 1,2,3,4,5,4,3,2,1

2010-03-05 Thread j verzani
kensuguro gmail.com> writes: > > > so basically, it's impossible to do with just seq() and rep().. Doesn't seem > like a good question for chapter 1... > > Also, problem 1.13 is even more crazy.. it asks you to build the fibonacci > sequence. Now I'm a programmer, and so went way ahead in t

Re: [R] Equation for model generated by auto.arima

2010-03-05 Thread Stephan Kolassa
Hi, In that case, I'd recommend reading a good book on time series analysis. "Forecasting: Methods and Applications" by Makridakis, Wheelwright and Hyndman is very accessible. Alternatively, there are probably tons of webpages on ARIMA, so google around. Best, Stephan testuser schrieb: Tha

Re: [R] How to assign week numbers to a time-series

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 4:23 PM, Hosack, Michael wrote: Hello everyone, My progress has stalled on finding a way of creating a somewhat complicated variable to add to my existing dataframe and I am hoping one of you could help me out. The dataframe below contains only a fraction of the data o

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Dylan Beaudette
On Friday 05 March 2010, Ted Harding wrote: > Again, I have the same attitude as Gabor. I don't need to know the > real identity of a poster, if they ask a sensible question (which > may well be elementary). > > As a contribution to the discussion, may I point out that many people > use different e

[R] How to assign week numbers to a time-series

2010-03-05 Thread Hosack, Michael
Hello everyone, My progress has stalled on finding a way of creating a somewhat complicated variable to add to my existing dataframe and I am hoping one of you could help me out. The dataframe below contains only a fraction of the data of my complete dataframe, but all of the variables. What I

Re: [R] Sweave and optional document sections

2010-03-05 Thread Max Kuhn
I haven't used windows for a while... when I did I used MikTex and TeXnicCenter. Max On Fri, Mar 5, 2010 at 2:27 PM, Aleksey Naumov wrote: > Max, > > Thank you very much! I'll give LaTeX a go (haven't used it in a while). Is > there a recommended distribution of LaTeX for Windows? > > Aleksey >

Re: [R] conditioning variable in panel.xyplot?

2010-03-05 Thread Seth W Bigelow
I'm stumped after an hour or so reading about subscripts in panel.xyplot. Apparently the panel function is executed for each subset of data in the main dataset (specified by the conditioning variable, 'site' in my example), and the 'subscripts' keyword passes a vector of the corresponding row n

Re: [R] zoo Package Question

2010-03-05 Thread Gabor Grothendieck
If your series is nearly regular, that is it is regularly spaced except for a few missing points then as.ts(z) will convert it to a regular ts series and as.zoo(as.ts(z)) does that and converts it back. You may need to add back the class. For example see: https://stat.ethz.ch/pipermail/r-help/201

Re: [R] Is it possible to recursively update a function?

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 11:06 AM, Seeker wrote: Thanks for your suggestion, Carl. Actually I am looking for a series of functions like these exp(-x)*.5^x exp(-x)*.5^x*(1-.4^x) exp(-x)*.5^x*(1-.4^x)*(1-.3^x) where(.5,.4,.3...) are from the coming results. Those are not really functions but rather

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Douglas Bates
On Fri, Mar 5, 2010 at 12:16 PM, Gabor Grothendieck wrote: > On Fri, Mar 5, 2010 at 12:58 PM, Matthew Dowle wrote: >> As far as I know you are wrong that there is no moderator.  There are in >> fact an uncountable number of people who are empowered to moderate i.e. all >> of us. In other words it

[R] REvolutions blog: February Roundup

2010-03-05 Thread David M Smith
I write about R every weekday at the Revolutions blog: http://blog.revolution-computing.com and every month I post a summary of articles from the previous month of particular interest to readers of r-help. http://bit.ly/9GoTVd announced the availability on YouTube of "What is R", a 4-part video b

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Ben Bolker
blue sky gmail.com> writes: I almost certainly shouldn't feed the trolls, but: 1. ?kruskal.test (listed in "see also" in ?wilcox.test) 2. One of the disadvantages of nonparametric tests is that it is in general difficult to generalize them to analogues of arbitrarily complex linear mixed mod

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Robert A LaBudde
A search on "bluesky...@gmail.com" shows the user is in Norfolk, VA, USA. At 01:26 PM 3/5/2010, John Sorkin wrote: The sad part of this interchanges is that Blue Sky does not seem to be amiable to suggestion. He, or she, has not taken note, or responded to the fact that a number of people belie

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Frank E Harrell Jr
Matthew Dowle wrote: John, So you want BlueSky to change their name to "Paul Smith" at "New York University", just to give a totally random, false name, example, and then you will be happy ? I just picked a popular, real name at a real, big place. Are you, or is anyone else, going to c

[R] how can I look at .Internal(model.matrix(t, data))?

2010-03-05 Thread Werner W.
Hi, I would like to see how model.matrix expands factor column to a set of dummy columns. I think that is done int .Internal(model.matrix(t, data)) which is called from model.matrix.default. But I have not idea how I can look at this function. How can I get to such internal functions? Thanks s

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Ted Harding
Again, I have the same attitude as Gabor. I don't need to know the real identity of a poster, if they ask a sensible question (which may well be elementary). As a contribution to the discussion, may I point out that many people use different email addresses for different purposes, so that any traf

[R] Assistance with pointers to code for B-spline derivatives (S-plus related).

2010-03-05 Thread Jeffrey Racine
Hi. I have been using the splines package for my work, in particular, the bs() function and associated predict() method. I now find myself in need of the derivatives of this beast. In the man page for predict.bSpline I found `predict(object, x, nseg=50, deriv=0, ...)' but alas this is not impl

Re: [R] End of line marker?

2010-03-05 Thread jonas garcia
Jim, Duncan and David, Thanks to you guys I manage to solve the problem and I have learnt a lot. Best regards J On Fri, Mar 5, 2010 at 4:55 AM, Duncan Murdoch wrote: > On 04/03/2010 11:40 PM, David Winsemius wrote: > >> On Mar 4, 2010, at 10:58 PM, Duncan Murdoch wrote: >> >> On 04/03/2010 10:

Re: [R] I can't find "rpart" help (linux)

2010-03-05 Thread Grześ
Thank you very much! :) -- View this message in context: http://n4.nabble.com/I-can-t-find-rpart-help-linux-tp1579403p1579846.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailm

Re: [R] install rJava in linux

2010-03-05 Thread Grześ
Thank you vary much! =^D -- View this message in context: http://n4.nabble.com/install-rJava-in-linux-tp1579395p1579852.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/lis

Re: [R] Is it possible to recursively update a function?

2010-03-05 Thread Seeker
Thanks for your suggestion, Carl. Actually I am looking for a series of functions like these exp(-x)*.5^x exp(-x)*.5^x*(1-.4^x) exp(-x)*.5^x*(1-.4^x)*(1-.3^x) where(.5,.4,.3...) are from the coming results. I am not sure how to construct such an irregular fucntion vector. On Mar 4, 6:40 pm, Carl

Re: [R] Is it possible to recursively update a function?

2010-03-05 Thread Ming Zhong
I was trying to replicate one CRM simulation. The following code works but seems redundant so I want to create a loop. O'Quigley CRM example 1# f1<-function(x) exp(-x) dose<-c(-1.47,-1.1,-.69,-.42,0.0,.42) p<-c(0.05,0.1,0.2,0.3,0.5,0.7) y<-c(0,0,1,0,1,1,0,0,0,0,0,0,1,0,0,0,0,0,1,0,0,1,0,1,

[R] zoo Package Question

2010-03-05 Thread testuser
I want to use the zoo package to convert an irregular time series to a regular one at every 15 min. interval. I would like to read in the csv file as a zoo object. When trying to do a seq(from,to,by), how can I specify the "from" to be the first time element in the zoo object and "to" to be the la

[R] Breadt-first-search algorithm

2010-03-05 Thread Thomas Jensen
Dear R-list, does anyone of you know whether there exist a breadth-first-search algorithm for R? Best, Thomas Jensen __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-pro

Re: [R] Sweave and optional document sections

2010-03-05 Thread Max Kuhn
That isn't hard to do with Sweave. Within you if branch, you can have your code write out any markup you choose (including \section). For figures, I would suggest writing the image file "manually" (i.e with the code chunk option fig = FALSE) and also write out the \includegraphics statements too.

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Matthew Dowle
John, So you want BlueSky to change their name to "Paul Smith" at "New York University", just to give a totally random, false name, example, and then you will be happy ? I just picked a popular, real name at a real, big place. Are you, or is anyone else, going to check its real ? We wan

Re: [R] About the interaction A:B

2010-03-05 Thread Jeff Laake
You are correct that you need to use ~-1+A:B. I use that all the time and just spaced it out when I was writing the response. Using ~A:B will produce one too many columns. Didn't follow your other question. You can always look at the result of model.matrix to see if it is correct or dig int

Re: [R] How to match vector with a list ?

2010-03-05 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Carlos Petti > Sent: Friday, March 05, 2010 9:43 AM > To: r-help@r-project.org > Subject: [R] How to match vector with a list ? > > Dear list, > > I have a vector of characters

Re: [R] How to match vector with a list ?

2010-03-05 Thread jim holtman
try this: > i <- c("a","a","b","b","b","c","c","d") > > j <- list(j1 = c("a","c"), j2 = c("b","d")) > > # create a matrix for mapping > map <- do.call(rbind, lapply(names(j), function(x) cbind(x, j[[x]]))) > # generate your output > map[match(i, map[,2]),1] [1] "j1" "j1" "j2" "j2" "j2" "j1" "j1"

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread John Sorkin
The sad part of this interchanges is that Blue Sky does not seem to be amiable to suggestion. He, or she, has not taken note, or responded to the fact that a number of people believe it is good manners to give a real name and affiliation. My mother taught me that when two people tell you that yo

Re: [R] [Rd] Changing the Prompt for browser()

2010-03-05 Thread Prof Brian Ripley
On Fri, 5 Mar 2010, Andrew Redd wrote: Is there a way that I can change the prompt for within a browser() call. I often use use code like with(obj1,browser()) Browse[1]> Is there a way that I can set it so that I can get something like with(obj1,browser(prompt="obj1")) obj1[1]> I know t

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Gabor Grothendieck
On Fri, Mar 5, 2010 at 12:58 PM, Matthew Dowle wrote: > As far as I know you are wrong that there is no moderator.  There are in > fact an uncountable number of people who are empowered to moderate i.e. all > of us. In other words its up to the responders to moderate.  The posting I think moderat

Re: [R] About the interaction A:B

2010-03-05 Thread blue sky
On Fri, Mar 5, 2010 at 11:41 AM, Jeff Laake wrote: > On 3/5/2010 9:19 AM, Frank E Harrell Jr wrote: >> >> You neglected to state your name and affiliation, and your question >> demonstrates an allergy to R documentation. >> >> Frank > > I agree with Frank but will try to answer some of your questi

Re: [R] About the interaction A:B

2010-03-05 Thread RICHARD M. HEIBERGER
On Fri, Mar 5, 2010 at 12:41 PM, Jeff Laake wrote: > On 3/5/2010 9:19 AM, Frank E Harrell Jr wrote: > >> You neglected to state your name and affiliation, and your question >> demonstrates an allergy to R documentation. > > A:B indicates a two-level naming scheme. The naming is identical in any

Re: [R] Sweave and optional document sections

2010-03-05 Thread Gabor Grothendieck
Its possible to do that using latex if/then/else macros which you can find in latex package ifthen; however, once you are getting to that level of complexity, Sweave is more of a hindrence than a help and is probably no longer the appropriate tool for the job. On Fri, Mar 5, 2010 at 12:55 PM, Alek

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Matthew Dowle
Frank, I respect your views but I agree with Gabor. The posting guide does not support your views. It is not any of our views that are important but we are following the posting guide. It covers affiliation. It says only that "some" consider it "good manners to include a concise signature spe

[R] Sweave and optional document sections

2010-03-05 Thread Aleksey Naumov
Dear R and Sweave users, Is there a way to have optional sections in a Sweave-generated report document, complete with section header(s), text and code chunks? In other words, I'd like for my report to include or omit certain sections based on the data itself. For example, If I examine the input

[R] How to match vector with a list ?

2010-03-05 Thread Carlos Petti
Dear list, I have a vector of characters and a list of two named elements : i <- c("a","a","b","b","b","c","c","d") j <- list(j1 = c("a","c"), j2 = c("b","d")) I'm looking for a fast way to obtain a vector with names, as follows : [1] "j1" "j1" "j2" "j2" "j2" "j1" "j1" "j2" I used : match <-

Re: [R] About the interaction A:B

2010-03-05 Thread Jeff Laake
On 3/5/2010 9:19 AM, Frank E Harrell Jr wrote: You neglected to state your name and affiliation, and your question demonstrates an allergy to R documentation. Frank I agree with Frank but will try to answer some of your questions as I understand it. First, model.matrix uses the options$contr

Re: [R] Running complete R script from Java

2010-03-05 Thread Romain Francois
You can source the script, e.g run the command : eval( "source( '" + script + "') " ) ; Questions about rJava/JRI are better on the stats-rosuda-devel mailing list: http://mailman.rz.uni-augsburg.de/mailman/listinfo/stats-rosuda-devel Romain On 03/05/2010 03:51 AM, Ralf B wrote: Is it poss

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread blue sky
I don't see which a link has the GENERAL and COMPLETE MATHEMATICAL description of nonparametric ANOVA for ARBITRARY MODEL. Would you please be specific which one does so? On Fri, Mar 5, 2010 at 10:52 AM, Jeremy Miles wrote: > Two links for you which will get your answer much quicker than a mailin

Re: [R] About the interaction A:B

2010-03-05 Thread Frank E Harrell Jr
You neglected to state your name and affiliation, and your question demonstrates an allergy to R documentation. Frank blue sky wrote: The following is the code for the model.matrix. But it still doesn't answer why A:B is interpreted differently in Y~A+B+A:B and Y~A:B. By 'why', I mean how R in

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread blue sky
On Fri, Mar 5, 2010 at 9:42 AM, David Winsemius wrote: > > On Mar 5, 2010, at 10:34 AM, Matthias Gondan wrote: > >>> This is your first of three postings in the last hour and they are all in >>> a category that could well be described as requests for tutoring in basic >>> statistical topics. I am

Re: [R] About the interaction A:B

2010-03-05 Thread blue sky
The following is the code for the model.matrix. But it still doesn't answer why A:B is interpreted differently in Y~A+B+A:B and Y~A:B. By 'why', I mean how R internally does it and what is the rational behind the way of doing it? And it didn't answer why in the model.matrix of Y~A, there are a-1 t

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Frank E Harrell Jr
Gabor Grothendieck wrote: I am happy to answer posts to r-help regardless of the name and email address of the poster but would draw the line at someone excessively posting without a reasonable effort to find the answer first or using it for homework since such requests could flood the list makin

Re: [R] Cross tabulation with fixed dimensions

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 11:43 AM, Tarmo Remmel wrote: Hello, This is a seemingly simple task, but it has been frustrating me for too long, so I am turning to this list for some help. I have two vectors of factors which are quite long; two simple examples are shown here: No, those are not "

Re: [R] Cross tabulation with fixed dimensions

2010-03-05 Thread Gabor Grothendieck
Those are numerical, not factors. If they were factors all levels would be represented. > a <- factor(c(1,2,3,4,5), levels = 1:5) > b <- factor(c(1,2,5,5,6), levels = 1:6) > table(a, b) b a 1 2 3 4 5 6 1 1 0 0 0 0 0 2 0 1 0 0 0 0 3 0 0 0 0 1 0 4 0 0 0 0 1 0 5 0 0 0 0 0 1 On Fri, M

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Jeremy Miles
Two links for you which will get your answer much quicker than a mailing list: http://lmgtfy.com/?q=non-parametric+anova+R or http://www.justfuckinggoogleit.com/search.pl?query=non+parametric+anova+R Jeremy On 5 March 2010 05:19, blue sky wrote: > My interpretation of the relation between 1-

[R] Changing the Prompt for browser()

2010-03-05 Thread Andrew Redd
Is there a way that I can change the prompt for within a browser() call. I often use use code like > with(obj1,browser()) Browse[1]> Is there a way that I can set it so that I can get something like > with(obj1,browser(prompt="obj1")) obj1[1]> I know that prompt is not a valid option for brows

[R] Cross tabulation with fixed dimensions

2010-03-05 Thread Tarmo Remmel
Hello, This is a seemingly simple task, but it has been frustrating me for too long, so I am turning to this list for some help. I have two vectors of factors which are quite long; two simple examples are shown here: > a <- c(1,2,3,4,5) > b <- c(1,2,5,5,6) If I produce a cross-tabulation of the

Re: [R] I can't find "rpart" help (linux)

2010-03-05 Thread Uwe Ligges
On 05.03.2010 10:51, Grześ wrote: Hi I have installed rpart in my Linux (PLD) but I don't know how I may find help conect this package? Here is my instalaction: install.packages("rpart",dependencies=TRUE) --- Please select a CRAN mirror for use in this session --- trying URL 'http://r.me

Re: [R] Is it possible to recursively update a function?

2010-03-05 Thread Uwe Ligges
On 05.03.2010 01:40, Carl Witthoft wrote: My foolish move for this week: I'm going to go way out on a limb and guess what the OP wanted was something like this. i=1, foo = x*exp(-x) i=2, foo= x^2*exp(-x) i=3, foo = x^3*exp(-x) . . . In which case he really should create a vector bar<-rep(na

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Gabor Grothendieck
I am happy to answer posts to r-help regardless of the name and email address of the poster but would draw the line at someone excessively posting without a reasonable effort to find the answer first or using it for homework since such requests could flood the list making it useless for everyone.

[R] Still on poLCA

2010-03-05 Thread 千早 ケンジ
Dear all, I have just sent a message asking about poLCA but I thought of another question I wanted to ask. I get the G^2 statistic in my output and want to test for its significance. I get that the degrees of freedom for the test are (S-1-p) where S is the number of different patterns observed

Re: [R] Running script with double-click

2010-03-05 Thread Gabor Grothendieck
Note that if you use Rterm.bat in http://batchfiles.googlecode.com in place of Rterm.exe then you won't have to adjust the Run action when you upgrade R. Rterm.bat is a single self contained Windows batch file that you simply place on your Windows path. On Thu, Mar 4, 2010 at 8:51 PM, Steve Taylo

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Ravi Varadhan
David, I agree with your sentiments. I also think that it is bad posting etiquette not to sign one's genuine name and affiliation when asking for help, which "blue sky" seems to do a lot. Bert Gunter has already raised this issue, and I completely agree with him. I would also like to urge the

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 10:34 AM, Matthias Gondan wrote: This is your first of three postings in the last hour and they are all in a category that could well be described as requests for tutoring in basic statistical topics. I am of the impression you have been requested not to engage in such

Re: [R] Nonparametric generalization of ANOVA

2010-03-05 Thread Matthias Gondan
> This is your first of three postings in the last hour and they are all in > a category that could well be described as requests for tutoring in basic > statistical topics. I am of the impression you have been requested not to > engage in such behavior on this list. For this question for instan

Re: [R] hier.part

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 9:52 AM, Marco Jorge wrote: Hi everyone, A beguinner question. - How shall i import 9 different ascii (created from a gis layer (arcmap grid)) into R to create a single dataframe for using in hier.part? Should i use read.table, then turn each created object into a single

Re: [R] for help on building a R package with several R function and a bunch of c, c++

2010-03-05 Thread Whit Armstrong
Pick up Rcpp, make your life easier. http://dirk.eddelbuettel.com/code/rcpp.html -Whit On Fri, Mar 5, 2010 at 9:19 AM, wrote: > Hope I can get quick help from here, I have a bunch of c, c++ included main > function and makefile. It works well on both UNIX and windows. I tried to > build R p

Re: [R] install rJava in linux

2010-03-05 Thread Ista Zahn
On Fri, Mar 5, 2010 at 4:41 AM, Grześ wrote: > > Hi > I'm a problem with installing "rJava" in Linux > > I have got in my system : java-sun, java-sun-tools, java-sun-jre > but if I wont install rJava I got a problem like this: > >>  install.packages("rJava", dependencies = TRUE) > --- Please selec

[R] Improved Nelder-Mead algorithm and Matlab's fminsearch. Was: Hi

2010-03-05 Thread Ravi Varadhan
The `fminsearch' in Matlab uses a version of the Nelder-Mead simplex search algorithm, which is a derivative-free search technique. The Nelder-Mead is also the default algorithm in optim(). Therefore, you can simply call optim() to get your job done. However, I would like to mentiion that the

Re: [R] Data frame column

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 8:26 AM, ManInMoon wrote: I have a big data frame and I have extracted a bit by doing: y<-d[1:10,6] y [1] Headings 0-49 -98 -49 -41 -120 -155 -204 -169 92329 Levels: -0 -1 -10 -100 -1000 -1 -10001 -10002 -10003 -10004 -10005 -10

Re: [R] Data frame query

2010-03-05 Thread David Winsemius
On Mar 5, 2010, at 7:40 AM, ManInMoon wrote: I have created a large dataframe (d) by getting data from file using read.table I now have 79 columns and 3 million rows. How can I plot the 6th column? I tried plot(d[,6]) but it doesn't look right. When I try to do just d[,6] the console get

Re: [R] Operating System Variable

2010-03-05 Thread Prof Brian Ripley
See ?R.version ?Sys.info ?.Platform On Windows, ?win.version It all depends what you mean by 'Operating System'. On Fri, 5 Mar 2010, Rob Helpert wrote: Hi. Is there an easy way (analogous to the $^O variable in perl) to find That's a pretty basic description. Even R.version$os often gives

Re: [R] Operating System Variable

2010-03-05 Thread Duncan Murdoch
On 05/03/2010 9:45 AM, Rob Helpert wrote: Hi. Is there an easy way (analogous to the $^O variable in perl) to find out what operating system R is currently using? .Platform is probably what you want. R.version gives related information. I don't think there's anything built in to determine p

  1   2   >