[R] how to apply the function cut( ) to many columns in a data.frame?

2007-02-28 Thread ahimsa campos-arceiz
Dear useRs, In a data.frame (df) I have several columns (x1, x2, x3xn) containing data as a continuous numerical response: df var x1x2 x3 1143 147 137 2 9393 117 316439 101 4123 11897 5 63 125 97 612983 124

Re: [R] problem with help.start and ?somefunction

2007-02-28 Thread Berwin A Turlach
G'day Jenny, On Wed, 28 Feb 2007 16:41:18 -0600 (CST) <[EMAIL PROTECTED]> wrote: > I am going to be teaching a workshop next week using R and > Bioconductor in one of our university's computer labs. They have > recently installed R 2.4.1 for me, and I'm checking all my scripts. I > just noticed t

Re: [R] prop.test or chisq.test ..?

2007-02-28 Thread Dylan Beaudette
On Wednesday 28 February 2007 01:07, Christoph Buser wrote: > Hi > > Some comments are inside. > > Dylan Beaudette writes: > > Hi everyone, > > > > Suppose I have a count the occurrences of positive results, and the > > total number of occurrences: > > > > > > pos <- 14 > > total <- 15 > >

Re: [R] What is a expression good for?

2007-02-28 Thread Gabor Grothendieck
You can evaluate it, differentiate it, pick apart its components, use it as a title or legend in a plot, use it as a function body and probably other things too: e <- expression(x+y) eval(e, list(x = 1, y = 2)) # 3 D(e, "x") e[[1]][[1]] # + e[[1]][[2]] # x e[[1]][[3]] # y plot(1:10, main = e)

Re: [R] Removing directory?

2007-02-28 Thread Henrik Bengtsson
Thanks! ...as the help indeed says: If 'recursive = FALSE' directories are not deleted, not even empty ones. It works; > dir.create("foo") > file.info("foo") size isdir mode mtime ctime atime foo0 TRUE 777 2007-02-28 15:50:33 2007-02

Re: [R] PROC TABULATE with R

2007-02-28 Thread Laurent Valdes
Hi Hadley, It seems aggregate does some trick !!! I will have a look at this ! Many many thanks ! Laurent 2007/2/28, hadley wickham <[EMAIL PROTECTED]>: > > You might want to have a look at the reshape package, > http://had.co.nz/reshape, which can create multi-dimensional cubes and > perform

[R] does function deriv take a vector of expressions?

2007-02-28 Thread toby909
Hi All I was just wondering if the function deriv also takes a vector of expressions If someone knows I would appreciate letting me know. Thanks Toby __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Efficient way to repeat rows (or columns) of a matrix?

2007-02-28 Thread Talbot Katz
Brilliant, thanks! -- TMK -- 212-460-5430home 917-656-5351cell >From: "James W. MacDonald" <[EMAIL PROTECTED]> >To: Talbot Katz <[EMAIL PROTECTED]> >CC: r-help@stat.math.ethz.ch >Subject: Re: [R] Efficient way to repeat rows (or columns) of a matrix? >Date: Wed, 28 Feb 2007 18:20:32 -

Re: [R] Packages in R for least median squares regression a nd computing outliers (thompson tau technique etc. )

2007-02-28 Thread D L McArthur
lalitha viswanath yahoo.com> writes: > > Hi > I am looking for suitable packages in R that do > regression analyses using least median squares method > (or better). Additionally, I am also looking for > packages that implement algorithms/methods for > detecting outliers that can be discarded befo

Re: [R] Efficient way to repeat rows (or columns) of a matrix?

2007-02-28 Thread James W. MacDonald
Hi Talbot, ma <- matrix(1:6, nrow = 2) rma <- ma[rep(1:2, c(2,3)),] Best, Jim Talbot Katz wrote: > Hi. > > If I have a vector, v_1, and another vector of positive integers, i_1, the > same length as v_1, then rep(v_1,i_1) will repeat v_i[j] exactly i_1[j] > times, like so: > > >>rep(c(1

Re: [R] Removing directory?

2007-02-28 Thread James W. MacDonald
I think you want to add a recursive = TRUE to your call to unlink(). Best, Jim Henrik Bengtsson wrote: > Hi, > > I'm trying to remove/delete a directory usingR. I've tried the > following with no success: > > % Rterm --vanilla > >>getwd() > > [1] "C:/Documents and Settings/hb/braju.com.R/ar

[R] Efficient way to repeat rows (or columns) of a matrix?

2007-02-28 Thread Talbot Katz
Hi. If I have a vector, v_1, and another vector of positive integers, i_1, the same length as v_1, then rep(v_1,i_1) will repeat v_i[j] exactly i_1[j] times, like so: >rep(c(1,2,3),c(3,2,1)) [1] 1 1 1 2 2 3 > I'd like to do the same sort of thing where I replace v_1 with a matrix, and the jth

[R] Removing directory?

2007-02-28 Thread Henrik Bengtsson
Hi, I'm trying to remove/delete a directory usingR. I've tried the following with no success: % Rterm --vanilla > getwd() [1] "C:/Documents and Settings/hb/braju.com.R/aroma.affymetrix/test" > dir.create("foo") > file.info("foo") size isdir mode mtime ctime

[R] problem with help.start and ?somefunction

2007-02-28 Thread drnevich
Hi all, I am going to be teaching a workshop next week using R and Bioconductor in one of our university's computer labs. They have recently installed R 2.4.1 for me, and I'm checking all my scripts. I just noticed that using the ?somefunction call to access the documentation for that function

Re: [R] What is a expression good for?

2007-02-28 Thread Alberto Vieira Ferreira Monteiro
Bert Gunter wrote: > > An expression object is the output of parse(), and so is R's representation > of a parsed expression. It is a type of list -- a parse tree for the > expression. This means that you can actually find the sorts of things you > mention by taking it apart as a list: > >> ex <- pa

Re: [R] What is a expression good for?

2007-02-28 Thread Bert Gunter
See V&R's S PROGRAMMING, esp. section 3.5; and section 6.1 and subsequent of the "R Language Definition." An expression object is the output of parse(), and so is R's representation of a parsed expression. It is a type of list -- a parse tree for the expression. This means that you can actually fi

Re: [R] question regression trees

2007-02-28 Thread David Farrar
Have a look at the function arguments. I think the function may not split a set of fewer than 20 objects, which is a default setting that can be changed. In addition to rpart, you might want to look at tree. regards, Farrar Wensui Liu <[EMAIL PROTECTED]> wrote: with seeing more

Re: [R] use tapply in a list which contains many data.frame

2007-02-28 Thread jim holtman
This should work on your data using the 'reshape' package; > nmr_19_model_aa<-read.table(url(" http://omega.psi.iastate.edu/data/nmr_19_model_aa.txt";),sep ="") > name.col<-c("pr","mod","mod_No","aa","aa_po","omega") > colnames(nmr_19

[R] What is a expression good for?

2007-02-28 Thread Alberto Monteiro
I mean, I can generate a expression, for example, with: z <- expression(x+y) But then how can I _use_ it? Is it possible to retrieve information from it, for example, that z is a sum, its first argument is x (or expression(x)) and its second argument is y? Alberto Monteiro _

Re: [R] map population density

2007-02-28 Thread Roger Bivand
On Wed, 28 Feb 2007, Barbara Diane-Spillmann wrote: > dear all, > > i?m trying to plot a map of population density in german communities by > using shape files. > to plot the polygons is not the problem but as soon as I want to add the > population data to the map I get the following message: >

Re: [R] Packages in R for least median squares regression and computing outliers (thompson tau technique etc.)

2007-02-28 Thread roger koenker
It's not often one gets needs to correct Gabor, but no, least median of squares is not the same as least absolute error regression. Take a look at the package robust if you want the lms. url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Depar

Re: [R] Packages in R for least median squares regression and computing outliers (thompson tau technique etc.)

2007-02-28 Thread Gabor Grothendieck
Try rq in quantreg using the default value for tau. On 2/28/07, lalitha viswanath <[EMAIL PROTECTED]> wrote: > Hi > I am looking for suitable packages in R that do > regression analyses using least median squares method > (or better). Additionally, I am also looking for > packages that implement a

[R] map population density

2007-02-28 Thread Barbara Diane-Spillmann
dear all, i?m trying to plot a map of population density in german communities by using shape files. to plot the polygons is not the problem but as soon as I want to add the population data to the map I get the following message: Fehler in if (attr(theMap$Shapes[[ii]], "nParts") == 1) { : Argum

Re: [R] question regression trees

2007-02-28 Thread Darin A. England
Claudia, Can you send us the actual call you are making to rpart()? The call to plot() doesn't really help. Darin On Wed, Feb 28, 2007 at 11:53:45AM -0500, Wensui Liu wrote: > with seeing more code and output, i guess your tree fails to grow. > > On 2/28/07, Claudia Romero <[EMAIL PROTECTED]> wro

[R] use tapply in a list which contains many data.frame

2007-02-28 Thread Aimin Yan
In order to make this question be clear, Here is my code >source("http://omega.psi.iastate.edu/data/nmr_19.R";) > class(temp) [1] "list" > names(temp) [1] "1A24" "1A57" "1A5J" "1A6X" "1AB7" "1AF8" "1AFI" "1AGG" "1AH9" "1AHL" [11] "1AJ3" "1AJW" "1AK7" "1ALG" "1AOY" "1APF" "1AUZ" "1AZJ" "1AZK"

Re: [R] Packages in R for least median squares regression and computingoutliers (thompson tau technique etc.)

2007-02-28 Thread Bert Gunter
Packages MASS and robustbase both have this functionality. There may also be others. Bert Gunter Genentech Nonclinical Statistics South San Francisco, CA 94404 650-467-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of lalitha viswanath Sent: Wednesday

Re: [R] matrix manipulations

2007-02-28 Thread Petr Klasterecky
Don't know what's wrong - it works: > data100 y X0 X1 X2 X3 ind [1,] 11 1 2.79581511 -23.5477 -33.6123061 1 [2,] 8 1 21.43289242 21.52826214 3.8415209 2 [3,] 6 1 6.18688631 21.51057247 -50.5547410 3 [4,] 12 1 -5.95172686

Re: [R] Package RGtk2, rattle, libatk-1.0.0.dll Errors

2007-02-28 Thread Liaw, Andy
The way I have that problem resolved is by installing the rggobi package using the command shown on http://www.ggobi.org/rggobi/, which is source("http://www.ggobi.org/download/install.r";) That will install all the things that Ggobi needs. Since rattle depends on rggobi, it's probably a good

[R] configure error on Mandriva 2007

2007-02-28 Thread Bricklemyer, Ross S
All, I am having difficulty installing R-2.4.1 on Mandriva Linux 2007 Discovery. Thanks to the help of Doug Bates I got further along the path (Thanks Doug!!). I now cannot figure out which line of code needs modified in order for configure to complete. Please see the attached output script

[R] Packages in R for least median squares regression and computing outliers (thompson tau technique etc.)

2007-02-28 Thread lalitha viswanath
Hi I am looking for suitable packages in R that do regression analyses using least median squares method (or better). Additionally, I am also looking for packages that implement algorithms/methods for detecting outliers that can be discarded before doing the regression analyses. Although some webs

Re: [R] matrix manipulations

2007-02-28 Thread Alberto Monteiro
Anup Menon Nandialath wrote: > > I have a basic question with R. I'm generating a set > of random variables and then combining them using the > cbind statement. The code for that is given below. > > for (i in 1:100) > { > y <- rpois(i,lambda=10) > X0 <- seq(1,1,length=i) > X1 <- rn

Re: [R] matrix manipulations

2007-02-28 Thread Peter Dalgaard
Anup Nandialath wrote: > Dear friends, > > I have a basic question with R. I'm generating a set > of random variables and then combining them using the > cbind statement. The code for that is given below. > > for (i in 1:100) > { > y <- rpois(i,lambda=10) > X0 <- seq(1,1,length=i) > X

[R] matrix manipulations

2007-02-28 Thread Anup Nandialath
Dear friends, I have a basic question with R. I'm generating a set of random variables and then combining them using the cbind statement. The code for that is given below. for (i in 1:100) { y <- rpois(i,lambda=10) X0 <- seq(1,1,length=i) X1 <- rnorm(i,mean=5,sd=10) X2 <- rnorm(

Re: [R] legend question

2007-02-28 Thread Sebastian P. Luque
On Wed, 28 Feb 2007 16:52:05 + (GMT), Jenny Barnes <[EMAIL PROTECTED]> wrote: > Hi folks, Do you mind if I ask a related question that I have been > having trouble with - how do you put the legend outside of the plot area > (to the bottom of the area - below the x-axis title)? Could anybody sh

[R] use tapply in a list which contains many data.frame

2007-02-28 Thread Aimin Yan
I have a list like the following: > class(temp) [1] "list" > names(temp) [1] "1A24" "1A57" "1A5J" "1A6X" "1AB7" "1AF8" "1AFI" "1AGG" "1AH9" "1AHL" [11] "1AJ3" "1AJW" "1AK7" "1ALG" "1AOY" "1APF" "1AUZ" "1AZJ" "1AZK" > head(temp[[1]]) pr modelaa omega 1 1A24 MODEL_1 1_ALA 84.47 2 1

Re: [R] legend question

2007-02-28 Thread Jenny Barnes
Hi folks, Do you mind if I ask a related question that I have been having trouble with - how do you put the legend outside of the plot area (to the bottom of the area - below the x-axis title)? Could anybody show me using the example given below: x <- seq(-pi, pi, len = 65) plot(x, sin(x), type

Re: [R] question regression trees

2007-02-28 Thread Wensui Liu
with seeing more code and output, i guess your tree fails to grow. On 2/28/07, Claudia Romero <[EMAIL PROTECTED]> wrote: > Hello, > This is my first time addressing such a big audience so apologies in > advance in case I fail to formulate this question. > > I am working with 13 species of trees, a

Re: [R] legend question

2007-02-28 Thread jim holtman
try: y<-c(1960, 1965, 1970, 1975) z<-c(1, 2, 3, 4) plot(y, z, type="l", col = 2) legend("topleft", "legend text", pch = 1, xjust = 0.5) On 2/28/07, Emili Tortosa-Ausina <[EMAIL PROTECTED]> wrote: > > Hi to all, > > I'm sorry for posting this question, I am sure I am missing something > importan

[R] question regression trees

2007-02-28 Thread Claudia Romero
Hello, This is my first time addressing such a big audience so apologies in advance in case I fail to formulate this question. I am working with 13 species of trees, and the data I have are: 1 continuous (phenolic concentration in xylem and in phloem) and 2 categorical variables: lineage (3 su

Re: [R] legend question

2007-02-28 Thread Stephen Tucker
Hi Emili, Even though you are calling your horizontal coordinate y, and vertical coordinate z, the first and second arguments to legend(), namely x and y, should be the horizontal and vertical coordinates, respectively; and they are given in user coordinates (e.g., legend()'s x should be between 1

Re: [R] legend question

2007-02-28 Thread Ranjan Maitra
On Wed, 28 Feb 2007 17:06:18 +0100 Emili Tortosa-Ausina <[EMAIL PROTECTED]> wrote: > y<-c(1960, 1965, 1970, 1975) > z<-c(1, 2, 3, 4) > plot(y, z, type="l", col = 2) > legend(x = -3, y = .9, "legend text", pch = 1, xjust = 0.5) your x and y are outside the plotting area. try using a different set

Re: [R] get more than get

2007-02-28 Thread Barry Rowlingson
Ingmar Visser wrote: > it's unclear what you want ... > but > get("x")^2 > does not give an error Neither does get("x^2") if you actually have an object called "x^2": > "x^2"=4.01 > get("x^2") [1] 4.01 but that would be a perverse thing to do. Barry __

Re: [R] Help on GAM

2007-02-28 Thread Wensui Liu
which library are you using for gam? On 2/28/07, Dacha Atienza <[EMAIL PROTECTED]> wrote: > 1) I have a semiparametric model, like > *Y~x1+s(x2)+s(x3)* > When I rum gam package I only obtained the estimates and the statistics of > the nonparametric part. How can I get the parametric part? Please c

[R] legend question

2007-02-28 Thread Emili Tortosa-Ausina
Hi to all, I'm sorry for posting this question, I am sure I am missing something important but after reading the documentation I cannot find where the problem is. I want to add a legend to a figure. If I use a simple example drawn from the R Reference Manual such as, for instance: x <- seq(-p

Re: [R] PROC TABULATE with R

2007-02-28 Thread hadley wickham
You might want to have a look at the reshape package, http://had.co.nz/reshape, which can create multi-dimensional cubes and perform some OLAP type functionality (although it's likely a lot slower than a fully-fledged OLAP, it will probably do the job) Hadley On 2/28/07, Laurent Valdes <[EMAIL PR

[R] Fwd: How to do GAM

2007-02-28 Thread Dacha Atienza
I have to evaluate how 4 different environmental variables could be use to explain the variability on the abundance of an specific species of zooplankton. I know that Generalized additive models are useful tools to evaluate this kind of relationships. 1) I want to know the difference to use the "g

Re: [R] survival analysis using rpart

2007-02-28 Thread Walter345
Thanks a lot for the reply, Terry! Concerning my question #3, I was thinking about the following scenario: Suppose that you have a data set of survival time data. We use rpart on two different subsets of the data set, where the columns contain the covariates. Say, for instance, we build the mod

[R] Help on GAM

2007-02-28 Thread Dacha Atienza
1) I have a semiparametric model, like *Y~x1+s(x2)+s(x3)* When I rum gam package I only obtained the estimates and the statistics of the nonparametric part. How can I get the parametric part? Please could you give me the complete comand to do it. 2) How are the negative coefficients identified. I

Re: [R] get more than get

2007-02-28 Thread Ingmar Visser
it's unclear what you want ... but get("x")^2 does not give an error hth, Ingmar On 28 Feb 2007, at 15:35, Alberto Monteiro wrote: > This must be a stupid question, but is there any "extension" of get? > > For example: > x <- 10 > get("x") # gives me 10 > get("x^2") # gives me an error > > Albert

Re: [R] get more than get

2007-02-28 Thread Uwe Ligges
Alberto Monteiro wrote: > This must be a stupid question, but is there any "extension" of get? > > For example: > x <- 10 > get("x") # gives me 10 > get("x^2") # gives me an error There is no object called "x^2" in the environments within your search path. If you want to calculate x^2, then t

Re: [R] get more than get

2007-02-28 Thread Barry Rowlingson
Alberto Monteiro wrote: > This must be a stupid question, but is there any "extension" of get? > > For example: > x <- 10 > get("x") # gives me 10 > get("x^2") # gives me an error 'get' really only gets R objects - you want to evaluate an expression - like this: > x=2 > eval(parse(text="x^2

Re: [R] topTable function from LIMMA

2007-02-28 Thread M Perez
Dear Roger, I think there is no a default value for the FDR in the toptable funtion. You will get a p-values corrected by the method you choose(fdr, bh). Briefly the shorted list of p-values will give you the FDR. It is up to you choose what is the estimated number of false positives you want to

[R] get more than get

2007-02-28 Thread Alberto Monteiro
This must be a stupid question, but is there any "extension" of get? For example: x <- 10 get("x") # gives me 10 get("x^2") # gives me an error Alberto Monteiro __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEAS

[R] rpart minimum sample size

2007-02-28 Thread Terry Therneau
Look at rpart.control. Rpart has two "advisory" parameters that control the tree size at the smallest nodes: minsplit (default 20): a node with less than this many subjects will not be worth splitting minbucket (default 7) : don't create any final nodes with <7

Re: [R] topTable function from LIMMA

2007-02-28 Thread James W. MacDonald
Hi Roger, Vallejo, Roger wrote: > Dear R-Help, > > > > I am using the function "topTable" from the LIMMA package. To estimate > adjusted P-values there are several options (adjust="fdr" , adjust="BH") > as shown below: > > > > topTable(fit, number = 10, adjust = "BH", fit$Name) > >

Re: [R] SEM - standardized path coefficients?

2007-02-28 Thread John Fox
Dear Tim, See ?standardized.coefficients (after loading the sem package). Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox >

[R] PROC TABULATE with R

2007-02-28 Thread Laurent Valdes
Hi ! with apply or tapply-like functions, is it possible to create multidimensional cubes with R ? Like with SAS and its function PROC TABULATE or OLAP ? Is there some functions or modules to access OLAP databases with R ? My idea is to create a package for that, since XMLA and JOLAP specificatio

[R] 2 Courses: (1) R/Splus Advanced Programming - San Francisco - (2) Data Mining: Practical Tools and Techniques in R/Splus - Salt Lake City

2007-02-28 Thread Sue Turner
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce our*** R/S Advanced Programming ***course in San Francisco on March 15-16 and *** Data Mining: Practical Tools and Techniques in R/Splus *** course in Salt Lake City on March 26-27 Ask for group discount and reserve your sea

[R] Piecewise linear regression with RMA

2007-02-28 Thread amsegura
Hi, I'm working in crustacean maturity, using breakpoint analysis. I used total length and Carapace length (linearly related) to find a break point, wich is actually the size at maturity. Both sizes are measured with error, so Modell II aplies. Piecewise regression uses model I to find breakpoint

[R] SEM - standardized path coefficients?

2007-02-28 Thread Tim Holland
Hello - Does anybody know how to get the SEM package in R to return standardized path coefficients instead of unstandardized ones? Does this involve changing the covariance matrix, or is there an argument in the SEM itself that can be changed? Thank you, Tim [[alternative HTML version d

[R] topTable function from LIMMA

2007-02-28 Thread Vallejo, Roger
Dear R-Help, I am using the function "topTable" from the LIMMA package. To estimate adjusted P-values there are several options (adjust="fdr" , adjust="BH") as shown below: topTable(fit, number = 10, adjust = "BH", fit$Name) I guess any of these options (fdr, BH, etc.) is using a de

[R] D() and deriv() accessing components of returned expression

2007-02-28 Thread toby909
Hi All I like the D() function since it directly gives me the derivative. D() however does not take multiple arguments like c("x1", "x2"), but deriv() does, but deriv() I am having trouble accessing the actual derivative function. So strange, its sitting direct in front of me but I cant access

Re: [R] LME without convergence

2007-02-28 Thread Douglas Bates
On 2/28/07, Ralf Sommerfeld <[EMAIL PROTECTED]> wrote: > Dear R-help list readers, > > I am fitting a mixed model using the lme function (R V 2.3.1 for > Windows). This is an example: > > dep<-c(25,40,33.33,60,70.83,72,71.43,50,40,53.33,64,54.17,60,53.57) > yes<-c(0,1,2,3,4,5,6,0,1,2,3,4,5,6) > tre

Re: [R] bootstrap

2007-02-28 Thread Frank E Harrell Jr
Indermaur Lukas wrote: > Hi, > I would like to evaluate the frequency of the variables within the best > selected model by AIC among a set of 12 competing models (I fit them with > GLM) with a bootstrap procedure to get unbiased results. So I would ike to > do the ranking of the 12-model-set 10

[R] delete selecting rows and columns

2007-02-28 Thread jastar
Hi, I'm working with a big square matrix (15k x 15k) and I have some trouble. I want to deleting selecting rows and columns. I use something like this: > sel_r=c(15,34,384,985,4302,6213) > sel_c=c(3,151,324,3384,7985,14302) > matrix=matrix[-sel_r,-sel_c] but it works very slow. Does anybody know

Re: [R] sort of OT: bootstrap tutorial

2007-02-28 Thread John Kane
Thank you Patrick. That looks very useful. --- Patrick Burns <[EMAIL PROTECTED]> wrote: > There is now a tutorial on bootstrapping and other > resampling > methods at: > > http://www.burns-stat.com/pages/Tutor/bootstrap_resampling.html > > Corrections and other suggestions are welcome. > > The

[R] LME without convergence

2007-02-28 Thread Ralf Sommerfeld
Dear R-help list readers, I am fitting a mixed model using the lme function (R V 2.3.1 for Windows). This is an example: dep<-c(25,40,33.33,60,70.83,72,71.43,50,40,53.33,64,54.17,60,53.57) yes<-c(0,1,2,3,4,5,6,0,1,2,3,4,5,6) treat<-c(1,1,1,1,1,1,1,0,0,0,0,0,0,0) #factor If I now fit a model wit

[R] bootstrap

2007-02-28 Thread Indermaur Lukas
Hi, I would like to evaluate the frequency of the variables within the best selected model by AIC among a set of 12 competing models (I fit them with GLM) with a bootstrap procedure to get unbiased results. So I would ike to do the ranking of the 12-model-set 10'000 times separately and calcula

Re: [R] adding graphics after dev.copy

2007-02-28 Thread Giampiero Salvi
Hi, sorry for the multiple posting ;) Here is the way to do it: pdf(file="file1.pdf") dev.control(displaylist='enable') plot(1:10) dev.copy(pdf,file="file2.pdf") points(2,3) dev.off() dev.off() Giampiero __ R-help@stat.math.ethz.ch mailing list https:/

Re: [R] adding graphics after dev.copy

2007-02-28 Thread Giampiero Salvi
Hi, now I see the problem is that it is not possible to run dev.copy from non-interactive devices (see the post in r-devel: "[Rd] feature enhancement request & patch: dev.control(displaylist='en (PR#3424)" Tue 08 Jul 2003 - 03:50:38 EST). I'll move this question to the r-devel list. Giampiero O

Re: [R] .C HoltWinters

2007-02-28 Thread Roger Bivand
On Tue, 27 Feb 2007, sj wrote: > Hello, > > I would like to look at the compiled C code behind HoltWinters from the > stats package. Is that possible? If so where do I find it? Please see the R help desk note in R News 6(4): 43-45, October 2006, by Uwe Ligges and Duncan Murdoch: "Accessing the

Re: [R] Package RGtk2, rattle, libatk-1.0.0.dll Errors

2007-02-28 Thread Hans W. Borchers
j.joshua thomas gmail.com> writes: > > Dear Group, > > I have followed the instructions from the link > http://datamining.togaware.com/survivor/Installing_GTK.html > However i couldn't fix the libatk01.0.0.dll application error > It is important that libatk-1.0.0.dll is in your path, that is

[R] sort of OT: bootstrap tutorial

2007-02-28 Thread Patrick Burns
There is now a tutorial on bootstrapping and other resampling methods at: http://www.burns-stat.com/pages/Tutor/bootstrap_resampling.html Corrections and other suggestions are welcome. The project started because a novice asked me about bootstrapping. My response was, "How dare you bug me while

Re: [R] compiling issues with Mandriva Linux 2007 Discovery

2007-02-28 Thread Paul Bivand
You may find that Mandriva Discovery is missing gfortran which is necessary for compiling R - having recently compiled it on a mandriva free edition. You may need to add appropriate package repositories to the package management system - http://easyurpmi.zarb.org is helpful for configuring Mandriva

Re: [R] factor documentation issue

2007-02-28 Thread Heinz Tuechler
At 09:41 28.02.2007 +1030, Geoff Russell wrote: >There is a warning in the documentation for ?factor (R version 2.3.0) >as follows: > >" The interpretation of a factor depends on both the codes and the > '"levels"' attribute. Be careful only to compare factors with the > same set of levels (in

[R] adding graphics after dev.copy

2007-02-28 Thread Giampiero Salvi
Hi, I am trying to create a plot and save it as pdf to the file file1.pdf, then add some details and save the new plot as file2.pdf. I would like to avoid repeating all the instructions needed to create file1.pdf when I create file2.pdf. This is what I have tried: pdf(file="file1.pdf") plot(1:10)

Re: [R] R in linux

2007-02-28 Thread Robert Lischke
hi, there are a couple of alternatives, i'm sure someone is going to bring up vi/emacs ;) i'd suggest you take a look at rkward (http://rkward.sourceforge.net/) and eclipse with the statet-plugin (http://www.eclipse.org/, http://www.walware.de/goto/statet). depending on which desktop your workin

Re: [R] prop.test or chisq.test ..?

2007-02-28 Thread Christoph Buser
Hi Some comments are inside. Dylan Beaudette writes: > Hi everyone, > > Suppose I have a count the occurrences of positive results, and the total > number of occurrences: > > > pos <- 14 > total <- 15 > > testing that the proportion of positive occurrences is greater than 0.5 > g

Re: [R] .C HoltWinters

2007-02-28 Thread Gregor Gorjanc
sj gmail.com> writes: > Hello, > > I would like to look at the compiled C code behind HoltWinters from the > stats package. Is that possible? If so where do I find it? Get the R source tarball from CRAN unpack it and search for the C code. Gregor __

Re: [R] Package RGtk2, rattle, libatk-1.0.0.dll Errors

2007-02-28 Thread Daniel Nordlund
Joshua, I have a few questions about the problems you are having which I will intersperse in your email below. I recall having some of the same difficulties you are having, which for me, were at least in part due to not following the installation directions exactly, in the correct order. > -

Re: [R] cluster analysis under contiguity constraints with R ?

2007-02-28 Thread Gavin Simpson
On Fri, 2007-02-16 at 13:10 +0100, Bellanger Lise wrote: > Hello, > > I would like to know if there is a function in an R library that > allows to do cluster analysis under contiguity constraints ? > > > Thank you very much for your answer ! > > Lise Bellanger Lise, One option

[R] no df to test the effect of an interaccion on a lmer mixed model

2007-02-28 Thread ahimsa campos-arceiz
Dear useRs, I am fitting a mixed model using the function lmer from the package lme4, but I have some problems when I try to test the effect of my factors of interest. First let me explain the structure of the model: I'm measuring animal movements. Explicitly, I am interested in displacement (str