[R] Anova with Scheffe Tests

Hi R-people, I am wanting to run Factorial ANOVA followed by Scheffe tests on some spatial subjective data. I'm comparing X-Y independent coordinates against x-y dependent coordinates. There are only four independent spatial coordinates that form a square. I am wondering whether I am doing the

Re: [R] Sweave and \input or \include LaTeX commands

On Tue, 01 Mar 2005 07:53:47 +0100, Gregor GORJANC (GG) wrote: Hi! What is wrong if there would be the same command? 2 reasons: 1) I try to keep the Sweave and LaTeX namespaces disjunct, i.e., I currently do not overload regular LaTeX commands with a different meaning, and I want

[R] Beginner - simple simulation

Hallo to everybody. I am new to the list and would appreciate some help in a basic first demo of how to use R for simulating a simple game; I would like my students to use R and this may stimulate their interest. The problem is simply: Two players (A and B) play the following game. Each player

Re: [R] Beginner - simple simulation

a simple approach could be: n - seq(10, 5000, 10) means - numeric(length(n)) for(i in seq(along=n)){ mat - sample(1:6, 2*n[i], TRUE); dim(mat) - c(n[i], 2) means[i] - mean(ifelse(!mat[,1]%%2, mat[,1], -mat[,2])) } means plot(n, means, xlab=sample size, ylab=mean) I hope it helps. Best,

Re: [R] Journal Quality R Graphs?

Werner Wernersen writes: the graphs look nice on the screen but when printed in black and white every color apart from black doesn't look very nice. My advice is: If you want a black-and-white or grayscale printout, don't plot in colors. -- Bjørn-Helge Mevik

[R] Variogram with time series

Hello, I'm using a dataset with unequally spaced time series and I'd want to know if there is in R some functions in order to calculate the autocorrelation function, because acf() in stats package cannot calculate it, because I have many missing data. And if so, is it necessary to organize

Re: [R] Sweave and \input or \include LaTeX commands

Friedrich.Leisch at tuwien.ac.at writes: I currently also plan only a \SweaveInput{}, I agree that its clearer to keep the two systems distinct rather than having subtle changes in functionality depending on the extension of the file in \input. While we are discussing I just wanted to bring

[R] part of name to Date

hi everybody, i try to extract a part of name to a date : like : VGT1_VGT2_CONTR_B020030401.H0V0.MIRto20030401 but the beginning of the files changes i have a list of files: [,1] [1,] VGT1_CONTR_B020020301.H0V0.MIR [2,]

Re: [R] part of name to Date

Please check string manipulation functions: - substr - strsplit Regards, Carlos Ortega On Tue, 01 Mar 2005 13:36:00 +0100, Jonathan Charrier [EMAIL PROTECTED] wrote: hi everybody, i try to extract a part of name to a date : like : VGT1_VGT2_CONTR_B020030401.H0V0.MIRto20030401

[R] Negative intercept in glm poisson model

Dear list, I'm trying to fit a glm model using family=poisson(link = identity). The problem is that the glm function fits a model with a negative intercept, which sounds like a nonsense to me, being the response a Poisson variable. From a previous discussion on this list I've understood that

Re: [R] Sweave and \input or \include LaTeX commands

On Tue, 1 Mar 2005 12:08:25 + (UTC), Gabor Grothendieck (GG) wrote: Friedrich.Leisch at tuwien.ac.at writes: I currently also plan only a \SweaveInput{}, I agree that its clearer to keep the two systems distinct rather than having subtle changes in functionality depending on

[R] Variogram with time series

Hello, I'm using a dataset with unequally spaced time series and I'd want to know if there is in R some functions in order to calculate the autocorrelation function, because acf() in stats package cannot calculate it, because I have many missing data. And if so, is it necessary to organize

Re: [R] draw random samples from empirical distribution

bogdan romocea wrote: Dear useRs, I have an empirical distribution (not normal etc) and I want to draw random samples from it. One solution I can think of is to compute let's say 100 quantiles, then use runif() to draw a random number Q between 1 and 100, and finally run runif() again to pull a

Re: [R] Negative intercept in glm poisson model

Federico Gherardini wrote: I'm trying to fit a glm model using family=poisson(link = identity). The problem is that the glm function fits a model with a negative intercept, which sounds like a nonsense to me, being the response a Poisson variable. Not really. The negative intercept is on the

Re: [R] Negative intercept in glm poisson model

Douglas Bates [EMAIL PROTECTED] writes: Federico Gherardini wrote: I'm trying to fit a glm model using family=poisson(link = identity). The problem is that the glm function fits a model with a negative intercept, which sounds like a nonsense to me, being the response a Poisson

Re: [R] Negative intercept in glm poisson model

On Tuesday 01 March 2005 14:31, Douglas Bates wrote: Not really. The negative intercept is on the scale of the linear predictor. The expected response, which is the exponential of the linear predictor, is always positive. Thank you very much for the quick response. Actually I used link =

[R] Help : delete at random

Hello I would like to delete some values at random in a data frame. Does anyone know how I could do? With best regards Caroline __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] Negative intercept in glm poisson model

Peter Dalgaard wrote: Douglas Bates [EMAIL PROTECTED] writes: Federico Gherardini wrote: I'm trying to fit a glm model using family=poisson(link = identity). The problem is that the glm function fits a model with a negative intercept, which sounds like a nonsense to me, being the response a

Re: [R] Help : delete at random

Caroline, You probably want to look at ?sample. Use sample to choose the rows for deletion then use: df.new = df[-sampled,] Sean On Mar 1, 2005, at 9:04 AM, Caroline TRUNTZER wrote: Hello I would like to delete some values at random in a data frame. Does anyone know how I could do? With best

Re: [R] Negative intercept in glm poisson model

On Tue, 1 Mar 2005 15:10:35 +0100, Federico Gherardini [EMAIL PROTECTED] wrote : Dear list, I'm trying to fit a glm model using family=poisson(link = identity). The problem is that the glm function fits a model with a negative intercept, which sounds like a nonsense to me, being the response a

Re: [R] Help : delete at random

On Tue, 2005-03-01 at 15:04 +0100, Caroline TRUNTZER wrote: Hello I would like to delete some values at random in a data frame. Does anyone know how I could do? With best regards Caroline The basic process is to randomly select row indices from the possible number of rows. If your data

Re: [R] Help : delete at random

Caroline TRUNTZER wrote: Hello I would like to delete some values at random in a data frame. Does anyone know how I could do? What about sample()-ing (if I understand at random correctly) a certain number of values from 1:nrow(data) and using the result as negative index the data.frame? Uwe

[R] GARCH

Hi, everyone! Is there a function to do single-variable GARCH in R? If yes, what library is it in? Thanks! Toby -- ** When Thomas Edison invented the light bulb he tried over 2000 experiments before he got it to

Re: [R] Help : delete at random

On Tue, 01 Mar 2005 15:04:17 +0100, Caroline TRUNTZER [EMAIL PROTECTED] wrote : Hello I would like to delete some values at random in a data frame. Does anyone know how I could do? Assuming your data.frame is named df: To delete index i, use df[-i, ] (i.e. the row selection is the negative of

[R] altering legend with plot(density(..))

I'm having a hard time modifying the legend that is emitted by default when calling plot(density(...)) and then legend(). I've looked at the docs and I'm not sure what I'm doing wrong, if anything. Any advice appreciated. Thanks, Brian __

[R] lattice -- panel order display

R-help, I'm using 'xyplot' in lattice package which plots length frecuencies by year (10). The order I get is not logical and the 'index.cond' argument to 'xyplot' is a bit cumbersome when it comes to plot a great deal of (in my case years). I have tried sorting the conditioning variable but

Re: [R] Help : delete at random

Hello, d - data.frame(a=c(2,3,4), b=c(2,4,1), c=c(3,5,6)) ## one NA s.r - sample(dim(d)[1], 1) s.c - sample(dim(d)[2], 1) d.na - d d.na[s.r, s.c] - NA d.na # Here a matrix is more comfortable by using sample. For multiple NA, you should write a loop, but to choose e.g. exact 4 values, it

Re: [R] part of name to Date

Here is one way x - c( VGT1_CONTR_B020020301.H0V0.MIR, VGT1_VGT2_CONTR_B020020611.H0V0.MIR, VGT1_VGT2_CONTR_B020030401.H0V0.MIR, VGT1_VGT2_CONTR_B020030711.H0V0.MIR, VGT1_VGT2_CONTR_B020031211.H0V0.MIR ) ( tmp1 - sapply( strsplit(x, split=_), tail, n=1 ) ) [1]

[R] repost my question of cda in case

Dear R-helpers: I sent this question 3 days ago but I didn't get any reply. In case this question was somewhat not seen by people who happpened to know the answer, I repost it here. Sorry for bother but I am kind of needing some help. BTW, if the question itself was not well expressed, please let

[R] write a library under 2.0.1 version

Hi there, I had written a library under R 1.9.0 and now I would like to import that library under R 2.0.1 Apparently it does not work; I can install the package, but when I try to read it the error is the following: Error in library(compvar) : 'compvar' is not a valid package -- installed 2.0.0

Re: [R] altering legend with plot(density(..))

Brian Hook wrote: I'm having a hard time modifying the legend that is emitted by default when calling plot(density(...)) and then legend(). I've looked at the docs and I'm not sure what I'm doing wrong, if anything. Any advice appreciated. So, can you tell us what is wrong with

Re: [R] lattice -- panel order display

Luis Ridao Cruz wrote: R-help, I'm using 'xyplot' in lattice package which plots length frecuencies by year (10). The order I get is not logical and the 'index.cond' argument to 'xyplot' is a bit cumbersome when it comes to plot a great deal of (in my case years). I have tried sorting the

Re: [R] Using mutiply imputed data in NLME

For example Will NLME work with mitools written by Thomas Lumley? Shige On Mon, 28 Feb 2005 13:02:54 -0500, Doran, Harold [EMAIL PROTECTED] wrote: Well, I may be missing something but why don't you call each dataframe from within a loop and save the results in a new data frame each time?

Re: [R] lattice -- panel order display

Here is a 'subset' of the ' test ' data . xyplot ( number ~ cm | as.factor(test$year) , data=test) test year cm number 34 1995 72 34 35 1995 73 37 36 1995 74 31 37 1995 75 13 38 1995 76 18 39 1995 77 4 40 1995 78 10 41 1995 79 6 42 1995 80 4

Re: [R] Help : delete at random

Might be slightly more interesting. If we want to generate values which are completely missing at random, then we can just simply sample all available index of a 2-d array. # simulate data # set.seed(1) # for reproducibility m - matrix( rnorm(12), nr=4, nc=3 ) m [,1] [,2]

Re: [R] lattice -- panel order display

Luis Ridao Cruz allegedly said on 3/1/2005 8:50 AM: R-help, I'm using 'xyplot' in lattice package which plots length frecuencies by year (10). The order I get is not logical and the 'index.cond' argument to 'xyplot' is a bit cumbersome when it comes to plot a great deal of (in my case years). I

Re: [R] lattice -- panel order display

On Tuesday 01 March 2005 08:50, Luis Ridao Cruz wrote: R-help, I'm using 'xyplot' in lattice package which plots length frecuencies by year (10). The order I get is not logical and the 'index.cond' argument to 'xyplot' is a bit cumbersome when it comes to plot a great deal of (in my case

[R] Two problems building a package

Hello all: I have written a few R scripts and am trying to turn them into a package for submission to CRAN. All of these scripts are R code only, no C or C++ or anything else. I'm working with R 2.0.1 running on a Windows XP machine. So far running rcmd install --build --docs=normal mypkge seems

Re: [R] lattice -- panel order display

Luis Ridao Cruz wrote: Here is a 'subset' of the ' test ' data . xyplot ( number ~ cm | as.factor(test$year) , data=test) test year cm number 34 1995 72 34 35 1995 73 37 [SNIP, was NTW a) not easily reproducible and b) too much waste of bandwith] I still don't see any problem? The

Re: [R] altering legend with plot(density(..))

On Tue, 01 Mar 2005 16:38:42 +0100, Uwe Ligges wrote: So, can you tell us what is wrong with set.seed(123) plot(density(rnorm(100))) and where the legend is? Maybe you don't like the defaults for xlab or main? Ah, yes, I was confusing xlab with the legend. Sorry for that. In addition,

Re: [R] write a library under 2.0.1 version

[EMAIL PROTECTED] wrote: Hi there, I had written a library You are talking about a *package*, I guess. under R 1.9.0 and now I would like to import that library under R 2.0.1 You have to (re-)install it from the sources. Apparently it does not work; I can install the package, but when I try to

Re: [R] altering legend with plot(density(..))

Brian Hook wrote: On Tue, 01 Mar 2005 16:38:42 +0100, Uwe Ligges wrote: So, can you tell us what is wrong with set.seed(123) plot(density(rnorm(100))) and where the legend is? Maybe you don't like the defaults for xlab or main? Ah, yes, I was confusing xlab with the legend. Sorry for that. In

Re: [R] altering legend with plot(density(..))

Brian == Brian Hook [EMAIL PROTECTED] on Tue, 1 Mar 2005 09:42:46 -0500 writes: Brian I'm having a hard time modifying the legend that is Brian emitted by default when calling plot(density(...)) Brian and then legend(). I've looked at the docs and I'm Brian not sure what

Re: [R] GARCH

Yes. Probably the best approach is if you go to the Rmetrics website http://www.rmetrics.org/ Patrick Burns Burns Statistics [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User) Tobias Muhlhofer wrote: Hi, everyone! Is there a

RE: [R] Temporal Analysis of variable x; How to select the outlier threshold in R?

I'm not sure I understand. You have financial data and want to throw away some outliers?? Why would you ever do this? First of all, I'd suggest you pay close attention to what the data is trying to say. Maybe your distribution is not normal after all (see tests for normality etc). Maybe you

Re: [R] write a **package** under 2.0.1 version

marta == marta [EMAIL PROTECTED] on Tue, 1 Mar 2005 16:23:17 +0100 (CET) writes: marta Hi there, marta I had written a library under R 1.9.0 and now I would like to import that marta library under R 2.0.1 package, package, package -- a library is something really different!

Re: [R] write a library under 2.0.1 version

On Tue, 1 Mar 2005 16:23:17 +0100 (CET), [EMAIL PROTECTED] wrote : Hi there, I had written a library under R 1.9.0 and now I would like to import that library under R 2.0.1 Apparently it does not work; I can install the package, but when I try to read it the error is the following: Error in

Re: [R] lattice -- panel order display

Luis Ridao Cruz allegedly said on 3/1/2005 9:53 AM: Here is a 'subset' of the ' test ' data . xyplot ( number ~ cm | as.factor(test$year) , data=test) snip data Luis, So what order are you expecting? Perhaps using `as.table = TRUE' is what you want?? xyplot(number ~ cm | as.factor(year), test,

Re: [R] lattice -- panel order display

I do the following: test$year - factor (test$year , ordered = TRUE) xyplot ( number~cm | year, data = test , type = l ) and still get the same. Luis Sundar Dorai-Raj [EMAIL PROTECTED] 01/03/2005 16:04:34 Luis Ridao Cruz allegedly said on 3/1/2005 8:50 AM: R-help, I'm using 'xyplot' in

[R] constraining initial slope in smoother.spline

Hello. I want to fit a smoother spline (or an equivalent local regression method) to a series of data in which the initial value of the 1st derivative (slope) is constrained to a specific value. Is it possible to do this? If so, how? Bill Shipley [[alternative HTML version

[R] Data in an object = Existing Objects??

Dear R-help... I am rather new so I would appreciate your help.. My question if the following.. I have generated (with R) a data frame DF looking like this.. print(DF) X1 X2 2 C_05_04 C_05_11 3 C_05_04 C_05_17 4 C_05_04 C_06_08 5 C_05_04 C_06_29 6 C_05_04 C_07_16 Where

[R] Fun and games with lowess()

This is not a help request, just an experience with R that I found amusing. I have a machine learning module that I teach which was originally all symbolic, but has a slowly growing numeric/statistical component. Today I taught a two part lecture on instance based methods for learning, with the

RE: [R] memory problem with mac os X

Yes, this came up last week as well. dist uses .C() to call code to compute the upper triangle, which is passed an empty R vector of size N(N-1)/2 to fill in. It returns a list containing the arguments passed in and assigns the result to another vector. I can only guess that because arguments are

[R] Atkinson's score

Hello Dear R Users, I would like to compute the Atkinson's score for a linear model $Y = \beta_0 + \beta_1 X_1$. I could not find a decent literature to do that. Is there a package in R that computes Atkinson's score? Thank you very much. New user, Aram. Many people would sooner die than

Re: [R] write a library under 2.0.1 version

On Tue, 1 Mar 2005 [EMAIL PROTECTED] wrote: Hi there, I had written a library under R 1.9.0 and now I would like to import that library under R 2.0.1 Apparently it does not work; I can install the package, but when I try to read it the error is the following: Error in library(compvar) : 'compvar'

Re: [R] GARCH

I just got 163 hits from www.r-project.org - search - R site search - garch. The first hit mentioned package tseries. The second also mentioned fSeries, fOptions. hope this helps. spencer graves Tobias Muhlhofer wrote: Hi, everyone! Is there a function to do single-variable

Re: [R] repost my question of cda in case

For linear discriminate analysis with only 2 classes, you get one new variable, that being a straight line between the means of the two groups. With 3 classes, the 3 points determine a plane defined by 2 lines or new functions = linear combinations of the original variables.

Re: [R] Using mutiply imputed data in NLME

On Tue, 1 Mar 2005, Shige Song wrote: For example Will NLME work with mitools written by Thomas Lumley? Yes, but since the coefficients are not just returned by the coef() function you will need to use the full MIextract/MIcombine approach rather than the shortcut of just using MIcombine() on

Re: [R] write a **package** under 2.0.1 version

Yet again, yet again, Martin. I told you... On Tue, 1 Mar 2005 17:28:04 +0100, Martin Maechler [EMAIL PROTECTED] wrote: marta == marta [EMAIL PROTECTED] on Tue, 1 Mar 2005 16:23:17 +0100 (CET) writes: marta Hi there, marta I had written a library under R 1.9.0 and now I

RE: [R] Temporal Analysis of variable x; How to select the outlier threshold in R?

--- bogdan romocea [EMAIL PROTECTED] wrote: I'm not sure I understand. You have financial data and want to throw away some outliers?? Why would you ever do this? I would select an outlier threshold, to extract a subset of the data x that had significant difference in financial contributions

[R] almost lower triangular matrices

I have output from a program which produces a distance matrix I want to read into a clustering program in R. The output is a .txt file and is 'almost' lower triangular in the sense that it is just the triangle below the diagonal. So for example a 4-by-4 distance matrix appears as, 1 2

Re: [R] lattice -- panel order display

On Tuesday 01 March 2005 10:41, Luis Ridao Cruz wrote: I do the following: test$year - factor (test$year , ordered = TRUE) xyplot ( number~cm | year, data = test , type = l ) and still get the same. What you haven't told us yet is what's wrong with what you get. -Deepayan

[R] Schedule Published for Data Mining Conference: New York City, March 28-30, Two full-days of Case Study Presentations

Apologies for cross posting - Salford Systems Data Mining 2005 New York, March 28-30, 2005 Focusing on the Contributions of Data Mining to Solving Real World Challenges

RE: [R] Data in an object = Existing Objects??

If you have the name of an object, you can use get to get the object itself. You didn't ask, and I don't know what your intended use is, but this looks like a difficult and inefficient way to keep track of things in R. You may want to look over some of the R introductions (on www.r-project.org)

[R] R with C#

Hello, do you have any examples or R being used with C#? I am currently working on a statistical project implemented in C#. I am researching using R for a third party plug-in to implement the more complicated statistical algorithms. Thanks Kyle Welling [[alternative HTML version

[R] start values for random effects in nlme

I am attempting to run an exponential decay model in nlme with 3 fixed and random effects (an initial value, an asymptote, and a decay rate). I am having trouble getting convergence now that I have added the asymptote to the model, so I was attempting to specify start values for the random

Re: [R] tkRplot under macos x

Luke, Thanks for the pointer to the postings. I still don't have it working. Here is what I have done: % R CMD INSTALL -l ~/Library/R/library/ tkrplot * Installing *source* package 'tkrplot' ... configure: creating ./config.status config.status: creating src/Makevars ** libs gcc

Re: [R] Data in an object = Existing Objects??

achilleas.psomas at wsl.ch writes: : : Dear R-help... : : I am rather new so I would appreciate your help.. : My question if the following.. : : I have generated (with R) a data frame DF looking like this.. : : print(DF) :X1 X2 : 2 C_05_04 C_05_11 : 3 C_05_04 C_05_17 : 4

Re: [R] almost lower triangular matrices

x - 1:6 # In real life x - scan(filename). m - matrix(0,4,4) m[row(m)col(m)] - x m - t(m) m [,1] [,2] [,3] [,4] [1,]0000 [2,]1000 [3,]2300 [4,]4560 The fiddle with the transposing is needed because R puts data into

Re: [R] almost lower triangular matrices

Michael Anyadike-Danes wrote: I have output from a program which produces a distance matrix I want to read into a clustering program in R. The output is a .txt file and is 'almost' lower triangular in the sense that it is just the triangle below the diagonal. So for example a 4-by-4 distance

[R] Users in Ukraine cyrillic support

Hello, everybody: My friends in Ukraine are starting a research lab at a national university and they asked what programs to use. I said R of course, and they then asked me 'what support does it have for Cyrillic'? i've done some snooping in the R website and all the references i find to foreign

[R] Reconstructing Datasets

Hi, Is it possible to recreate smoothed data sets in R, by performing a PCA and then reconstructing a data set from say the first 2/3 EOFs? I've had a look in the help pages and don't seem to find anything relevant. Thanks in advance, Laura Laura Quinn Institute of Atmospheric Science School

[R] AIX 5.1

Hello, I was wondering if someone who has successfully built R (either 2.0.2005 or 1.91) on Platform 3/4 and compiled for AIX 5.1 could send on their compilation flags, tricks/tips ... Thanks, Gavin __ R-help@stat.math.ethz.ch mailing list

Re: [R] Users in Ukraine cyrillic support

As from the next release it has support of almost all human languages. Previews of that are available now. On Tue, 1 Mar 2005, Paul Johnson wrote: Hello, everybody: My friends in Ukraine are starting a research lab at a national university and they asked what programs to use. I said R of course,

Re: [R] Users in Ukraine cyrillic support

On Tue, 1 Mar 2005, Paul Johnson wrote: Hello, everybody: My friends in Ukraine are starting a research lab at a national university and they asked what programs to use. I said R of course, and they then asked me 'what support does it have for Cyrillic'? The development version (to become 2.1.0)

[R] looking for 3D plotting with real-time rotation

I would like to plot points in 3D and then be able to rotate the plot in real time in order to understand the distribution of the points using the visual parallax that results from the rotation. The R-package scatterplot3d will plot in 3D, but I've not found a way to rotate the results in real

Re: [R] AIX 5.1

On Tue, 1 Mar 2005, Gavin La Rowe wrote: I was wondering if someone who has successfully built R (either 2.0.2005 or 1.91) on Platform 3/4 and compiled for AIX 5.1 could send on their compilation flags, tricks/tips ... People have built on AIX 5.1, and their flags and tips are in the R-admin

Re: [R] looking for 3D plotting with real-time rotation

Dear Ben, The scatter3d() function in the Rcmdr package uses rgl to plot rotatable point clouds and regression surfaces. The function is usable with or without the Rcmdr GUI. You'll find an illustrative plot at http://socserv.socsci.mcmaster.ca/jfox/Courses/S-course/index.html. I hope this

Re: [R] na.strings in readLines or is.na?

Dear Ben, Since readLines() returns the lines in the file as character strings, it wouldn't be appropriate for it to substitute for missing values on input. One approach would be to replace -99.99 with NA in the strings, but a simpler method would be to deal with the data frame (say DF): DF[DF

Re: [R] almost lower triangular matrices

Output from which program? If the output is of class dist, then as.matrix should give you what you want: set.seed(1) X - array(rnorm(12), dim=c(4,3)) (dX - dist(X)) 1 2 3 2 1.6598683 3 0.9720025 2.4600033 4 2.2666735 2.2149274

Re: [R] Users in Ukraine cyrillic support

Paul Johnson [EMAIL PROTECTED] writes: Hello, everybody: My friends in Ukraine are starting a research lab at a national university and they asked what programs to use. I said R of course, and they then asked me 'what support does it have for Cyrillic'? i've done some snooping in the R

[R] How to convert a factor to a numeric?

Hi all -- I've got two columns, both of which correspond to three factor levels (e.g., column1 is a, b, or c; column2 is x, y, or z). I'd like to generate a third column, consisting on whether the two factors are correctly aligned for a given case (in this example, a corresponds to x, b to

Re: [R] write a library under 2.0.1 version

On Tue, 1 Mar 2005, Prof Brian Ripley wrote: On Tue, 1 Mar 2005 [EMAIL PROTECTED] wrote: I had written a library under R 1.9.0 and now I would like to import that library under R 2.0.1 Apparently it does not work; I can install the package, but when I try to read it the error is the

RE: [R] How to convert a factor to a numeric?

# If fact1 and fact2 are your factors, let prm be the permutation such that # levels(fact2) corresponds to (aligns to) levels(fact1)[prm] . In your example, the permutation is # apparently the identity, (1:3). #Then levels(fact2)[prm[fact1]]==fact2 ## does what you want. I wouldn't be surprised

Re: [R] How to convert a factor to a numeric?

Chris Bergstresser chris at subtlety.com writes: : : Date: Tue, 01 Mar 2005 16:45:36 -0600 : From: Chris Bergstresser [EMAIL PROTECTED] : To: r-help@stat.math.ethz.ch : Subject: [R] How to convert a factor to a numeric? : : : Hi all -- : : I've got two columns, both of which

Re: [R] Anova with Scheffe Tests

Hi, I don't think there are any packages on CRAN that implement Scheffe's test. If you don't mind using another multiple comparisons procedure, you could look at ?TukeyHSD and/or the multcomp package. Alternatively, you could write your own function to do Scheffe's test. At least one other

[R] Re: R-help Digest, Vol 25, Issue 1

[EMAIL PROTECTED] wrote: Subject: [R] Ask for your help I got following error information when I run the haplo.glm program.Could you tell me which wrong was happened in my program?Many thanks. Shanchun fit.gaus - haplo.glm(y ~ SEX+geno, family = gaussian, + data=my.data, locus.label=label, control

RE: [R] Anova with Scheffe Tests

I stand corrected, although confidence.ellipse is a plotting function, and may not be quite what the questioner had in mind. Cheers, Simon. See confidence.ellipse() in the car() package. (Found from an R site search on Scheffe) -- Bert Gunter Genentech Non-Clinical Statistics South San

Re: [R] A problem about outer()

Thanks for the informative comments and the detailed explanation. - Original Message - From: Adaikalavan Ramasamy [EMAIL PROTECTED] To: Feng Chen [EMAIL PROTECTED] Cc: R-help r-help@stat.math.ethz.ch Sent: Tuesday, March 01, 2005 12:25 AM Subject: Re: [R] A problem about outer() You

[R] How to force specific variables in all the models using regsubsets

Dear friends, I am trying to use R to do Best Subset Regression, I can not force seldom variables in all the subsets models. R is always giving wrong variables when I include force.in option. For example, I like to include X2 X4 in all the models, but I always get X2 X3. Could

Re: [R] Reconstructing Datasets

Laura Quinn a écrit : Hi, Is it possible to recreate smoothed data sets in R, by performing a PCA and then reconstructing a data set from say the first 2/3 EOFs? I've had a look in the help pages and don't seem to find anything relevant. See function reconst in package ade4. Best, Renaud -- Dr

Re: [R] Reconstructing Datasets

On Tue, 2005-03-01 at 20:30 +, Laura Quinn wrote: Hi, Is it possible to recreate smoothed data sets in R, by performing a PCA and then reconstructing a data set from say the first 2/3 EOFs? I've had a look in the help pages and don't seem to find anything relevant. It's not in the

Re: [R] Reconstructing Datasets

On Wed, 2005-03-02 at 08:30 +0200, Jari Oksanen wrote: On Tue, 2005-03-01 at 20:30 +, Laura Quinn wrote: Hi, Is it possible to recreate smoothed data sets in R, by performing a PCA and then reconstructing a data set from say the first 2/3 EOFs? I've had a look in the help pages

Re: [R] return from nested function?

On Feb 25, 2005, at 12:34 PM, [EMAIL PROTECTED] wrote: Is is possible from within a function to cause its caller to return()? This snippet may be of interest: f = function(x) { + print(f) + g(return()) + print(end of f) + } g=function(x) {print(g) + x + print(end of g) + } f(1) [1] f [1] g NULL

Re: [R] write a library under 2.0.1 version

Paul Roebuck wrote: On Tue, 1 Mar 2005, Prof Brian Ripley wrote: On Tue, 1 Mar 2005 [EMAIL PROTECTED] wrote: I had written a library under R 1.9.0 and now I would like to import that library under R 2.0.1 Apparently it does not work; I can install the package, but when I try to read it the error

[R] Leaps regsubsets

Hello I am trying to use all subsets regression on a test dataset consisting of 11 trails and 46 potential predictor variables. I would like to use Mallow's Cp as a selection criterion. The leaps function would provide the required output but does not work with this many variables (see below).