[R] how to use biplot to show PC1 and PC3 or else?

hello: I am trying to use biplot show relationship between PC1 and PC3 of princomp. There is example showing how to do it with PC1 and PC2. But I am totally out of idea for PC1 and PC3. Thanks very much for your kind help. Best Wishes, yours, xin __

RE: [R] how to use biplot to show PC1 and PC3 or else?

xin, data(USArrests) biplot(princomp(USArrests)) # This plots the default 1st and 2nd components biplot(princomp(USArrests),choices=c(1,3)) # use choices to plot other components...in this case 1st and 3rd Bill -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]

Re: [R] graph algorithms in R

You found graph, but there is also an R API for the Boost.Graph libraries, RBGL, there. best, -tony Andre Skusa [EMAIL PROTECTED] writes: Dear R users, just a quick question: Is there a reliable and good graph library for R, eg. with shortest path algorithms on adjacency matrixes? I

[R] Hankel Matrix

Hello is there a function in R to generate Hankel Matrix? thanks __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

RE: [R] Mixed-effects model for nested design data

Dear Han, I am using nlme for data from nested design. That is, tows are nested within trip, trips nested within vessel, and vessels nested within season. I also have several covariates, say tow_time, latitude and depth My model is y = season + tow_time + latitude + depth +

Re: [R] Stichprobe bilden

Thomas Lumley wrote: On Thu, 29 Apr 2004, Thomas Petzoldt wrote: If you want to make a subsample, you can select cases (rows) using sub(), see ?sub for details, e.g. somethink like: I think you mean subset(). sub() is a regular expression substitution function. Absolutely sure! I should not try

[R] absolute value

Hello could you please tell me what is the function to get the absolute value of the real or complex number. most of other languages it is abs(x) , what is it in r? I did few searhces in the help docs for no avail. thanks F.J __ [EMAIL PROTECTED]

[R] randomization test

Hello all, I am conducting a randomization test on a given dataset. One of the covariates, gender, is randomly assigned 1000 times to create a randomization dataset (rand.data). To these 1000 datasets, I fit a full model (see below) and the aim is to generate a distribution of LRT statistic

Re: [R] Rtemp directories accumulating over time

On Wed, 28 Apr 2004 [EMAIL PROTECTED] wrote: On 28 Apr 2004 at 13:53, Prof Brian Ripley wrote: On Tue, 27 Apr 2004 [EMAIL PROTECTED] wrote: . . I don't see why: package base is never unloaded so that hook function is never run. (Indeed, no package/namespace is unloaded except by

Re: [R] Fedora 1 RPM Packages

On Wed, 28 Apr 2004, Ian Wallace wrote: That makes sense since most people will build with the defaults and that's what we'd like packages for the main stream. I'm already working around it ... but thought I would mention that some of us do like the shared library. I hadn't even given any

[R] I'm trying to use package ts (decompose). How do you set up the data/ See attached. thanks

InDATA -read.table(C:/Data/May 2004/season.txt,header=T) X - decompose(InDATA) print(X) Period Connections Q1 67519 Q2 69713 Q3 68920 Q4 69452 Q1 70015 Q2 59273 Q3 57063 Q4 65596 Q1 73527 Q2 58586 Q3

Re: [R] R-crash using read.shape (maptools)

On Wed, 28 Apr 2004, Danny Heuman wrote: Hi Herry, On Thu, 29 Apr 2004 12:20:44 +1000, you wrote: Hi List, I am trying to read a large shapefile (~37,000 polys) using read.shape [winxp, 1gig ram, dellbox). I receive the following error: AppName: rgui.exe AppVer: 1.90.30412.0

[R] problem with pca

Hello I have some kind of problem with R. I want to do a pca (I work on R and ade4 library) but R doesn't wan to realize my pca because of my data: I have some binary variable (presence/absence of birds, hole, vegetals...) and metric variable (cliff 's size...) and each time I have an error

Re: [R] outer

K Fung [EMAIL PROTECTED] writes: outer(m,n,cord.mag) this gives me a matrix of zeroes outer(1,2,cord.mag) this gives the right value on the other hand Your cord.mag function probably doesn't vectorize. Try Cord.mag - function(m,n) mapply(cord.mag, m, n) outer(m, n, Cord.mag) --

[R] Sweave in makefile under Win XP

Hi all. To automate the production of a document I want to have something like: test.tex : test.Rnw Sweave(test.Rnw) | R --no-restore --no-save in a makefile. This however gives an error as R states: Sweave(test.Rnw) : not found. I use R-1.9.0pat (040424 I think) on Win XP. Any

[R] accessing information in lists

I've created a dataframe containing multiple ACF lists through the command rev.acf-apply(rev.matrix, 2, acf, na.action=na.contiguous, lag.max=12, plot=FALSE) where rev.matrix is an n by t matrix containing n time series in columns. I'd now like to pull out only the ACF information and store it

Re: [R] Sweave in makefile under Win XP

On Thu, 29 Apr 2004, Ladelund, Steen wrote: Hi all. To automate the production of a document I want to have something like: test.tex : test.Rnw Sweave(test.Rnw) | R --no-restore --no-save in a makefile. This however gives an error as R states: Sweave(test.Rnw) : not found.

Re: [R] Sweave in makefile under Win XP

Ladelund, Steen [EMAIL PROTECTED] writes: Hi all. To automate the production of a document I want to have something like: test.tex : test.Rnw Sweave(test.Rnw) | R --no-restore --no-save in a makefile. This however gives an error as R states: Sweave(test.Rnw) : not found.

Re: [R] accessing information in lists

Try: *= rev.matrix-matrix(rnorm(40),10,4) rev.acf-apply(rev.matrix, 2, acf, na.action=na.contiguous, lag.max=12, plot=FALSE) x.acf-lapply(rev.acf,function(x) x[[acf]][,,1]) matrix(unlist(x.acf),ncol=4) output-start Thu Apr 29 09:56:22 2004 [,1] [,2][,3][,4]

Re: [R] Matrix efficiency in 1.9.0 vs 1.8.1

Stephen == Stephen Ellner [EMAIL PROTECTED] on Wed, 28 Apr 2004 17:19:38 -0400 writes: Stephen I'm seeking some advice on effectively using the Stephen new Matrix library in R1.9.0 for operations with Stephen large dense matrices. I'm working on integral Stephen operator

[R] line transect method

Hello, I'm a final-year student in statistics and I deal with as subject of thesis rare and elusive populations. I want make a simulation, in which I want generate a population and estimate its density with line transect method using R language. Is there someone that can send me the program

Re: [R] line transect method

It might be worth looking at: http://www.ruwpa.st-and.ac.uk/estimating.abundance/ best, Simon On Thu, 29 Apr 2004, Villirillo wrote: Hello, I'm a final-year student in statistics and I deal with as subject of thesis rare and elusive populations. I want make a simulation, in which I want

Re: [R] I'm trying to use package ts (decompose). How do you set up the data/ See attached. thanks

Package ts in R contains no functionality, it having been moved to package stats. The function ts() was not in package ts in older versions either. I have no idea what your data refer to (quarters, probably, but of which years?). Assumming you know the time base, see the help page ?ts for how

Re: [R] label separators in boxplots with conditioning

Webb Sprague wwsprague at ucdavis.edu writes: : : Hi R-helpers, : : I have a data.frame with three columns (lots more reps though in each), : like so: : : 'FOO''BAR''RESULT' : 1 .0175 : 1 .0512 : 1.1.01100 :

[R] GUI checkbox linked to entry

Hi , i am using R 1.8 and tcltk library in order to make a GUI and i would like to know if it is possible to link a check box to one or more entry dynamically (and eventually other widgets) . ( in concrete terms : if the check box is ticked the entry is avaible for input and unavaible when not

Re[2]: [R] simple repeated measures model: dumb user baffled!

On Thursday, April 29, 2004, 10:17:46 AM, Christophe Pallier wrote: CP Hi, CP Here is how I would go: CP suj - gl(40,4) CP order - gl(4,40) CP time - gl(4,1,160) CP cond - factor(c(rep(c('N','E','E','N'),10), CP rep(c('N','E','N','E'),10), CP rep(c('E','N','E','N'),10), CP

RE: [R] Dummies in R

From: Susana Bird Dear all, my problem is following. I know Stata, but currently I have to use R. I have the opposite of your problem: I don't know Stata, but I don't have to use Stata. 8-) Could You please help in finding the analogy to R. (1) creating of City-Dummy. I

[R] Maen squared deviation

Dear list, i'm looking for how to get the mean squared deviation in R. I tried a lot of (complicated) instructions, sometimes the worked, but now i'm looking for an easier oneHope you can help - greetings. __ [EMAIL PROTECTED] mailing list

Re: [R] problem with pca

Hello, if you want some help, you have to be more precise ! What is the error message, could you give an example ? read the posting guide at http://www.R-project.org/posting-guide.html (I am not sure that this question is devoted to R-help, it would be perhaps better to ask your question to the

Re:[R] p-values

On Thu, 29 Apr 2004, John Maindonald wrote: This is, of course, not strictly about R. But if there should be a decision to pursue such matters on this list, then we'd need another list to which such discussion might be diverted. Ted Harding started such a list (stats-discuss) quite some

[R] RPART drawing the tree

Hello, I am using the RPART library to find patterns in HIV mutations regarding drug-resistancy. My data consists of aminoacid at certain locations and two classes resistant and susceptible. The classification and pruning work fine with Rpart. however there is a problem with displaying the

Re: [R] Dummies in R

On Thu, 29 Apr 2004, [iso-8859-1] Susana Bird wrote: Dear all, my problem is following. I know Stata, but currently I have to use R. Could You please help in finding the analogy to R. (1) creating of City-Dummy. (2) Create a Time-Dummy-Variable Andy Liaw has described how to do this, but

Re: [R] RPART drawing the tree

On Thu, 29 Apr 2004, Rob Kamstra wrote: I am using the RPART library to find patterns in HIV mutations regarding drug-resistancy. My data consists of aminoacid at certain locations and two classes resistant and susceptible. The classification and pruning work fine with Rpart. however

Re: [R] Dummies in R

On Thursday 29 of April 2004 07:49, Susana Bird wrote: my problem is following. I know Stata, but currently I have to use R. Could You please help in finding the analogy to R. I am too clueless about Stata, but I have another hint (check Introduction to R whether I have guessed correctly for

[R] Problems in plot

Hello, I have R1.9.0 under Windows XP. My program plots several plots using x11() par(cex = 0.75) .. x11() par(cex = 0.75) .. x11() par(cex = 0.75) .. x11() par(cex = 0.75) .. Sometimes, one of them generates a small frame only with title area R graphics: Device X

[R] geoR y-scale problem

I am attempting to use the Bayesian kriging function in geoR by Ribeiro and I am encountering a problem with the image command. The program is not responsive to the ylim command, such that: image(temp.bayes,loc=pred.grid,col=rainbow(25, start=.6, end=.1),xlab=Longitude,ylab=Latitude,

[R] memory problems with lm

Hello list, I've seen the recent discussions documenting problems with lm. I have encountered the following problem. I use WinXP Pro with service pack 1, and R 1.9.0, on a XEON 2GHz, with 1GB of RAM. eff.fro std.dev mean NSTRDSP 7.403749e-01 1.215686e-01 CPFGEP

[R] Entering times around the start of daylight savings time

I'm having problems entering dates and times around when daylight savings time starts. If I type (on R 1.8.1 on Gentoo Linux) ISOdatetime(2004,4,4,0:4,0,0,GMT) [1] 2004-04-03 19:00:00 EST 2004-04-03 20:00:00 EST [3] 2004-04-03 22:00:00 EST 2004-04-03 22:00:00 EST [5] 2004-04-03 23:00:00 EST

RE: [R] memory problems with lm

Can you show us the output of str(eff.fro)? Do you have other things in the global environment or the search path that's taking up memory? What does gc() say? Andy From: Adrian Dragulescu Hello list, I've seen the recent discussions documenting problems with lm. I have encountered

Re: [R] Problems in plot

The primary graphics device under Windows is called *windows* not *x11*. Something in your Windows setup is sometimes failing to choose a reasonable window size. I have never seen that, and suspect it is nothing to do with it, but please use windows() and see if the problem vanishes. On Thu, 29

[R] graph algorithms in R

Dear R users, just a quick question: Is there a reliable and good graph library for R, eg. with shortest path algorithms on adjacency matrixes? I already found a graph package as part of Bioconductor. Are there more? If anyone knows and would tell me, I appreciate very much! Cheers, Andre --

Re: [R] memory problems with lm

This may or may not be the same problem (which is already solved). But please read the section on BUGS in the R FAQ and set up a reproducible example. Then try out the current version of r-patched (one dated tomorrow or later, to be safe) and see if the problem recurs. If it does, please

RE: [R] memory problems with lm

If I enforce the variables to be numeric it works fine. str(eff.fro) `data.frame': 6 obs. of 2 variables: $ std.dev: num [, 1:6] 7.40e-01 9.06e+00 4.70e+05 1.02e+06 1.37e+07 ... ..- attr(*, dimnames)=List of 1 .. ..$ : chr NSTRDSP CPFGEP WSWOLF NPILGRIM ... $ mean : num 1.22e-01

RE: [R] memory problems with lm

I believe Prof. Ripley is right. The problem is $ std.dev: num [, 1:6] 7.40e-01 9.06e+00 4.70e+05 1.02e+06 1.37e+07 ... ..- attr(*, dimnames)=List of 1 which looks like an array, rather than a vector. Andy From: Adrian Dragulescu [mailto:[EMAIL PROTECTED] If I enforce the

[R] R plot to HTML image map

Hi, I automatically create some R matplots and now I would like to automatically create HMTL image maps and make the plots clickable. Do you have any suggestions how to do this? regards, R. Alberts __ [EMAIL PROTECTED] mailing list

Re: [R] Entering times around the start of daylight savings time

You might do better using seq.POSIXt. Your basic problem is that your sequence includes a time that doesn't exist, 02:00 on the transition day. seq.POSIXt(ISOdatetime(2004,4,4,0,0,0),by='hour',len=5) [1] 2004-04-04 00:00:00 PST 2004-04-04 01:00:00 PST 2004-04-04 03:00:00 PDT 2004-04-04

Re: [R] R plot to HTML image map

Rudi Alberts wrote: Hi, I automatically create some R matplots and now I would like to automatically create HMTL image maps and make the plots clickable. Do you have any suggestions how to do this? What part of the plots do you want clickable? Lines are normally a bit thin for clickable

Re:[R] p-values

On 29-Apr-04 Thomas Lumley wrote: On Thu, 29 Apr 2004, John Maindonald wrote: This is, of course, not strictly about R. But if there should be a decision to pursue such matters on this list, then we'd need another list to which such discussion might be diverted. Ted Harding started such

[R] Mixed-effects model for nested design data

Hi, I am using nlme for data from nested design. That is, tows are nested within trip, trips nested within vessel, and vessels nested within season. I also have several covariates, say tow_time, latitude and depth My model is y = season + tow_time + latitude + depth + vessel(season) +

RE: [R] p-values

One might begin by considering _conditional_ p-values as elaborated by Hubbard and Bayarri and especially Sellke, Bayarri, and Berger. Record Number: 1545 @article{ Hubbard2003, Author = {Hubbard, R. and Bayarri, M. J.}, Title = {Confusion over measures of evidence ($p$)'s versus errors

Re: [R] GUI checkbox linked to entry

Julien Glenat wrote: Hi , i am using R 1.8 and tcltk library in order to make a GUI and i would like to know if it is possible to link a check box to one or more entry dynamically (and eventually other widgets) . ( in concrete terms : if the check box is ticked the entry is avaible for input

[R] Probability(Markov chain transition matrix)

Hello, My name is Maria, MBA student in Sanfransisco, USA. In my credit scoring class, I have hard time making transition matrix, which explains probability especially in relation to Markov chain model. It is regarding people's monthly credit payment behavior. Does R have function to caculate

[R] Plot.lme error

Dear All; Attempting to reproduce Figure 4.15 of MEMSS (p. 171) in R using plot(Wafer, outer = ~ Wafer) yields an error: plot(Wafer, outer = ~ Wafer) Error in order(na.last, decreasing, ...) : Argument lengths differ The plot() produces the figure without

Re: [R] Plot.lme error

On Thursday 29 April 2004 15:45, Yang, Richard wrote: Dear All; Attempting to reproduce Figure 4.15 of MEMSS (p. 171) in R using plot(Wafer, outer = ~ Wafer) yields an error: plot(Wafer, outer = ~ Wafer) Error in order(na.last, decreasing, ...) : Argument

Re: [R] Probability(Markov chain transition matrix)

On 29 Apr 2004 at 20:27, Maria Gu wrote: Hello, My name is Maria, MBA student in Sanfransisco, USA. In my credit scoring class, I have hard time making transition matrix, which explains probability especially in relation to Markov chain model. It is regarding people's monthly credit payment

Re: [R] Stichprobe bilden

[EMAIL PROTECTED] wrote: Liebe Alle, ich habe ien Problemm. ich habe einen Panel-Datensatz ueber die 185 Staaten (nach Code von IMF) und von 1948-1999. Ich habe auch 12 Variablen in diesem Datensatz, die ich fuer die Regression benoetige. ich muss ein sub-sample bilden fuer die 65 Staaten und

Re: [R] Stichprobe bilden

On Thu, 29 Apr 2004, Thomas Petzoldt wrote: If you want to make a subsample, you can select cases (rows) using sub(), see ?sub for details, e.g. somethink like: I think you mean subset(). sub() is a regular expression substitution function. -thomas

[R] ARIMAX-ARCH

hi everybody i have to fit an ARIMAX-ARCH model, if anybody knows about the existence of a function that can do it or has any suggestion please mail me. thanks a lot Niccolo __ [EMAIL PROTECTED] mailing list

[R] ARIMAX-ARCH

Please try the library(tseries) There are functions for garch in there. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html