Re: [R] Identify and plotting symbols.

On Thu, 20 Sep 2007, Rolf Turner wrote: > I have been trying, unsuccessfully, to use identify() to (simply) return > a list of the indices of points clicked on and overplot (with say a > solid dot) each clicked-on point so that I can see where I've been. > I.e. I don't want to see the indices p

Re: [R] SHLIB problem under Vista

This is a known problem with your compiler installation, discussed in the R-admin manual. From the version in R-patched: There are known problems with one of the compilers sets in this toolset on Windows Vista: see the workaround below. (With that workaround it has been used on both 32-

Re: [R] SEM - standardized path coefficients? - significant style

Dear John Being an R user is like having Christmas every day and having a direct line to Santa for present suggestions. How about if path.diagram.sem had the ability to print non-significant paths in, for example dotted style? I had a go myself but couldn't see how to get (summary(model)$coeff)$"P

[R] Ambiguities in vector

Hello all you helpful people out there! I am stil R Beginner using R 2.5.1 on a Apple Power Book G4 with Mac OS X 10.4.10 . Perhaps you haven´t understood my question in the mail yesterday. So I will try to describe my problem in a different way You see the tables. I would like to test the var

Re: [R] RE : Must be easy, but haven't found the function (numerical integration)

Hello Ravi, I was also trying to write down such a calculation but I did not manage (I didn't use diff(x)). Thank you for the solution. Ptit Bleu. Ravi Varadhan wrote: > > Here is a simple trapezoidal rule integrator: > > x <- time > y <- value > > area <- sum(diff(x)*(y[-1]+y[-length(y)]))

[R] problem with generalized singular value decomposition using LAPACK

Hi All, I'm trying to run generalized singular value decomposition (GSVD) function from LAPACK library. Basically my problem is that I can not run it for large matrices, I get a memory error. I'm using R 2.5.1. I tried this on intel centos5 machines with 2 GB memory and 8 GB memory. I have u

Re: [R] Smooth line in graph

Both worked, thanks! N. Fernandez > Katharine Mullen: > > require(splines) > x<-1:5 > y <- c(0.31, 0.45, 0.84, 0.43, 0.25) > yy <-predict(interpSpline(x, y)) > plot(x, y) > lines(yy) > > Greg Snow: > lines(spline(x,y, method='n', n=250)) __ R-help@

[R] R

I am going to fit my model with R-2.5.1.tar my model is y=fixed effect+b(x) and I have three column in my data file (y a x) and about 20 observation(y) but I had some problems. I was wondering If you could help me the error massage is: > res <- read.table("pro.sdf") > fm <- glm(y=A+x,data=res)

Re: [R] R

On 20-Sep-07 08:45:43, mohammad sahebhonar wrote: > I am going to fit my model with R-2.5.1.tar > my model is y=fixed effect+b(x) and I have three column in my > data file (y a x) and about 20 observation(y) > but I had some problems. I was wondering If you could help me > the error massage is

[R] combine mathematical expressions with plain text

Hello, I would like to create a mtext. The argument 'text' should combine the mathematical expressions 'bar' with the plain text ' = 3.07'. Is that possible? mtext(text=expression(bar(x)) ...) ' = 3.07' Thanks for any ideas Tobias -- View this message in context: http://www.nabble.com/combine

Re: [R] Smooth line in graph

Sorry, I answered too quickly. It worked with the "simplified" example I provided but not with non-regular intervals in x: x<-c(-45,67,131,259,347) y <- c(0.31, 0.45, 0.84, 0.43, 0.25) plot(x,y) lines(spline(x,y, method='n', n=250)) #or: lines(predict(interpSpline(x, y))) Produce the same decre

Re: [R] problem with generalized singular value decomposition using LAPACK

> I'm trying to run generalized singular value decomposition (GSVD) > function from LAPACK library. Basically my problem is that I can not run > it for large matrices, I get a memory error. I'm using R 2.5.1. I tried > this on intel centos5 machines with 2 GB memory and 8 GB memory. I have >

Re: [R] combine mathematical expressions with plain text

plot(0,0,"n") mtext( expression( bar(x) == 3.07) ) but you can also simply 'paste' things: text(0,0,labels=expression(paste(bar(omega), " = 1")),srt=90) __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] Need help on "date"

This mail is continuation of my previous one. i have some raw data from Excel which was carried to R: data = read.delim(file="clipboard", header=T) > data Date Price 1 09/01/05 365 2 09/02/05 360 3 09/03/05 360 4 09/05/05 370 5 09/06/05 370 6 09/08/05 365 7 09/09/05 365

Re: [R] combine mathematical expressions with plain text

Hi, have you tried: mtext(expression(bar(x) == 3.07), ...) or mtext(bquote(bar(x) == 3.07), ...) ? Fernando Mayer. squall44 escreveu: > Hello, > > I would like to create a mtext. The argument 'text' should combine the > mathematical expressions 'bar' with the plain text ' = 3.07'. Is that >

Re: [R] Need help on "date"

Since this is a time series you might want to look at the zoo package: Lines <- "Date Price 09/01/05 365 09/02/05 360 09/03/05 360 09/05/05 370 09/06/05 370 09/08/05 365 09/09/05 365 09/10/05 365 09/12/05 365 09/13/05 360 09/14/05 360 09/15/05 360 " library(zoo) # replace

[R] referencing packages?

Hi I know how to referenc R in a scientific paper - but is there a standardised way to reference packages? Thanks Rainer -- NEW EMAIL ADDRESS AND ADDRESS: [EMAIL PROTECTED] [EMAIL PROTECTED] WILL BE DISCONTINUED END OF MARCH Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation Biology (

Re: [R] referencing packages?

use citation("pkg-name") On Thu, 20 Sep 2007, Rainer M. Krug wrote: > Hi > > I know how to referenc R in a scientific paper - but is there a > standardised way to reference packages? > > Thanks > > Rainer > > -- > NEW EMAIL ADDRESS AND ADDRESS: > > [EMAIL PROTECTED] > > [EMAIL PROTECTED] WILL BE

Re: [R] referencing packages?

Thanks a million Rainer Katharine Mullen wrote: > use citation("pkg-name") > > On Thu, 20 Sep 2007, Rainer M. Krug wrote: > >> Hi >> >> I know how to referenc R in a scientific paper - but is there a >> standardised way to reference packages? >> >> Thanks >> >> Rainer >> >> -- >> NEW EMAIL ADDR

[R] Cutting & pasting help examples into script window

Hi All, When I cut & paste help file examples into a script window, about half the time it pastes as a single long line. The steps I follow are: 1. Open a help file e.g. ?data.frame. 2. Select the examples at the bottom. 3. Choose File: Copy. 4. Return to the console. 5. Choose File: New script

[R] SEM - singularity error

Good morning, I am trying to develop a structural equation model of snake abundance using habitat variables. In attempting to estimate the model using the "sem" package in R version 2.4.0, I receive the following error message: "Error in solve.default(C) : system is computationally singular: reci

[R] How to find "p"(proportion) in binomial(n, p)?

Hi, I got a problem. I am trying to find "p" in binomial. X~bin(n, p) I want to find value "p", so that Pr(X <= k) <= alpha. Here, n, k are known. Thank you for helping me with this! Catherine _ [[replacing trailing spam]] ___

[R] The Survey package again

Good afternoon! 1 I have a bit of a problem with the Survey package, but one that is more theoretical than practical. If for example on wishes to use some additional information which comprises of 4 variables with 7, 5,5 and 5 categories (resulting in 875 post-strata if the whole crossin

Re: [R] Identify and plotting symbols.

Rolf Turner auckland.ac.nz> writes: ... > I tried > > ind <- identify(x,y,labels=rep("\021",length(x)),offset=0) Have you tried "\225" or "\244" or "\370" (which can be previewed by typing > "\225" etc), and clicking just above the plotted datapoints? D L McArthur, UCLA Sch of Medicine

[R] ggplot and xlim/ylim

Hello everyone, I am (happily) using ggplot2 for all my plotting now and I wondered is there is an easy way to specify xlim and ylim somewhere when using the ggplot syntax, as opposed to the qplot syntax. Eg. qplot(data=mtcars,y=wt, x=qsec,xlim=c(0,30)) <-> ggplot(mtcars, aes(y=wt, x=qsec

Re: [R] Cutting & pasting help examples into script window

On Thursday 20 September 2007 15:08:45 Muenchen, Robert A (Bob) wrote: MR > Hi All, MR > MR > When I cut & paste help file examples into a script window, about half MR > the time it pastes as a single long line. MR > MR > Any ideas? MR > MR > Thanks, MR > Bob MR > MR > P.S. What would really be sli

[R] Comprehensive New Learning Resource for R

Hi Folks, I have created a website with over 70 R tutorials on a wide range of subjects. The site is aimed at data analysts who have an idea about what they would like to do, but don't know how to do it in R. The site should be particularly useful for users of statistical packages like SAS, SPS

[R] Package Survey

Hello, How I use the function as.svrepdesign without memory.size problems? desenho_npc_JK <- as.svrepdesign(desenho_npc,type="JKn") Error: cannot allocate vector of size 161.3 Mb In addition: Warning messages: 1: Reached total allocation of 1022Mb: see help(memory.size) 2: Reached tota

Re: [R] Must be easy, but haven't found the function (numerical integration)

Ptit Bleu. try the 'adapt' function for multi-dimensional numerical integration. library(adapt) Adaptive Numerical Integration in 2–20 Dimensions John Checked by AVG Free Edition. 11:53 __ R-help@r-project.org mailing list https://stat.ethz.ch/ma

Re: [R] Cutting & pasting help examples into script window

Stephan Grosse replied: > > What I do not understand is why you not just type example(yourcommand)? > > Stefan That's a good question. I want to play around with variations of the examples rather than run them exactly as they are. Thanks, Bob __ R-h

[R] comparing survival curves(weighted log rank test)

hie I know I can use the log rank test for comparing survival curves but my problem is that the G(rho) by fleming & harrington(1982) which is the one in R seems to be inadequate are they any other test ie fleming & harrington (1991) as there are some some differences on the survival curves

Re: [R] ggplot and xlim/ylim

Hi Jiho, To figure this out, you need to think about exactly what you are doing when you change the limits - what part of the plot are you changing? You are changing the scales, and if you have a look at http://had.co.nz/ggplot2/scale_continuous.html, you'll see you need something like: + scale_y

Re: [R] Package Survey

On Thu, 20 Sep 2007, Rita Sousa wrote: How I use the function as.svrepdesign without memory.size problems? desenho_npc_JK <- as.svrepdesign(desenho_npc,type="JKn") Error: cannot allocate vector of size 161.3 Mb There is currently no easy way to affect the amount of memory that this uses, un

Re: [R] Comprehensive New Learning Resource for R

On 20.09.2007, at 09:24, Rob Kabacoff wrote: > I hope that you find this useful. Please feel free to link to the > site. Very nicely done from what I saw by browsing the site a little. Good colors, good content, easy to navigate. Thanks! You might want to add self-organizing maps as a cross-

Re: [R] Comprehensive New Learning Resource for R

Also use of T and F is bad form since its possible to define variables T and F. Use TRUE and FALSE. On 9/20/07, Ted Harding <[EMAIL PROTECTED]> wrote: > Excellent! Thanks. Just what I need for referring > friends to. > > However, your very first feedback: > > Under http://www.statmethods.net/inpu

Re: [R] warnings although R CMD check runs without any problem

Can you give examples? Uwe Christian Kohler wrote: > Dear list members, > > as I am new to R, I would kindly like to ask for your help. > > I checked a package with the 'R CMD check' command and all worked out > well, (only 'OK' status). Loading its vignette afterwards with > 'source("package.

Re: [R] Comprehensive New Learning Resource for R

Excellent! Thanks. Just what I need for referring friends to. However, your very first feedback: Under http://www.statmethods.net/input/datatypes.html you have: ##Vectors a <- c(1,2,5.3,6,-2,4) # numberic vector b <- c("one","two","three") # character vector c <- c(T,T,T,F,T,F) #logical vector

[R] warnings although R CMD check runs without any problem

Dear list members, as I am new to R, I would kindly like to ask for your help. I checked a package with the 'R CMD check' command and all worked out well, (only 'OK' status). Loading its vignette afterwards with 'source("package.Rcheck/inst/doc/vignette.R")' yields a message saying that several (

Re: [R] SEM - standardized path coefficients? - significant style

Dear Steve, I'll file this suggestion for when I next revise the sem package. I'm unlikely to get to it soon since I rarely use the graph-drawing program dot and would have to figure out how to do what you want. As well, I think that I'd more likely allow the user to specify the line type for each

[R] Problem with lars

Hi! Ok, probably my last mail was too long. Questions are: * Did someone try lars on data with much more components than data points? Did you observe bad overfitting, too? Any other observations? * Why might the lars (least angle regression from lars package) method produce beta values that

[R] It works but I don't know how (nice scientific notation with Lattice)

Hi, I had problems to display my results in a nice way (according to me). Thanks to posts of people of this forum, I solved my problems (ticks with a scientific notation and how to write nicely "V/m2" in the label of the x-axis) but ... I don't understand the script. Can somebody explain me how

[R] SUPER PARAMAGNETIC CLUSTERING

Dear all, we would like to know if there exist some routines (software) to build models based on the super paramagnetic clustering. Best regards, Silvia Figini -- Silvia Figini Department of Statistics "L. Lenti" University of Pavia Via Strada Nuova, 65 - 27100 - Pavia, Italy tel.+39- 0382-9846

Re: [R] warnings although R CMD check runs without any problem

Uwe Ligges wrote: > Can you give examples? > > Uwe > > > Christian Kohler wrote: > >> Dear list members, >> >> as I am new to R, I would kindly like to ask for your help. >> >> I checked a package with the 'R CMD check' command and all worked out >> well, (only 'OK' status). Loading its vignette

Re: [R] warnings although R CMD check runs without any problem

Hi, Christian Kohler wrote: > Uwe Ligges wrote: >> Can you give examples? >> >> Uwe >> >> >> Christian Kohler wrote: >> >>> Dear list members, >>> >>> as I am new to R, I would kindly like to ask for your help. >>> >>> I checked a package with the 'R CMD check' command and all worked out >>> we

Re: [R] warnings although R CMD check runs without any problem

On Thu, 20 Sep 2007, Christian Kohler wrote: > Dear list members, > > as I am new to R, I would kindly like to ask for your help. > > I checked a package with the 'R CMD check' command and all worked out > well, (only 'OK' status). Loading its vignette afterwards with > 'source("package.Rcheck/ins

[R] Plotmath issue superscript "-"

Dear List, I'm trying to typeset some chemical ions in axis labels. These have both super and subscript components, and for some, I need a superscript "-". In LaTeX I might use $NO_3^-$ to do the typesetting, but I'm having a problem getting the correct invocation for expression: > expression(NO^

Re: [R] warnings although R CMD check runs without any problem

Prof Brian Ripley wrote: > On Thu, 20 Sep 2007, Christian Kohler wrote: > > >> Dear list members, >> >> as I am new to R, I would kindly like to ask for your help. >> >> I checked a package with the 'R CMD check' command and all worked out >> well, (only 'OK' status). Loading its vignette afterw

Re: [R] Plotmath issue superscript "-"

Try: plot(0,0,"n") text(0,0,expression( {NO[3]}^'-')) __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

Re: [R] Plotmath issue superscript "-"

Gavin Simpson wrote: > Dear List, > > I'm trying to typeset some chemical ions in axis labels. These have both > super and subscript components, and for some, I need a superscript "-". > In LaTeX I might use $NO_3^-$ to do the typesetting, but I'm having a > problem getting the correct invocation f

Re: [R] Smooth line in graph

I don't know what SigmaPlot and Excel are doing for you, but I would guess that they are not doing cubic splines (as a general rule, when R and Excel differ, it is safest to assume that R is not the one doing something wrong) Often differences between packages are due to differences in assumptio

[R] Portfolio Optimization

Hello, I would like to solve a portfolio optimization problem in R. As far as I searched, I found the example of "solve.QP" &"portfolio.optim". In my understanding, both of them are based on given expected return, finding the minimum variance. Is there a way of doing this in an opposite way?i.e m

Re: [R] SEM - singularity error

[EMAIL PROTECTED] wrote: > Good morning, > > I am trying to develop a structural equation model of snake abundance using > habitat variables. In attempting to estimate the model using the "sem" package > in R version 2.4.0, I receive the following error message: > > "Error in solve.default(C) : s

[R] help with making a function of scatter plot with multiple variables

Dear list, I have done a scatter plot of multiple variables in the same graph, with different col and pch. I managed to do it with the following code but not know how to make a function of these so that next time if I want to do similar graph but with new variables, I dont have to copy th

Re: [R] Plotmath issue superscript "-"

If you need a subscript as well, I like plot(0, main=quote({NO^'\x96'}[3])) Peter Ehlers Peter Dalgaard wrote: > Gavin Simpson wrote: >> Dear List, >> >> I'm trying to typeset some chemical ions in axis labels. These have both >> super and subscript components, and for some, I need a superscri

Re: [R] SEM - singularity error

Good afternoon Chuck, I really appreciate your help. I just ran the factor analysis and did not receive the singularity error (please see results below). Do you happen to have any additional ideas or suggestions? Thank you so much. Nicki FACTANAL RESULTS: Call: factanal(x = thsi, factors = 1)

[R] how can I attach a variable stored in

Hi - Any help would be greatly appreciated. I'm loading a list variable that's stored in an .RData file and would like attach it. I've used attach( ), but that only lets me see the variable that's stored in the file. As the variable name is of the form "comp.x.x", I've tried using attach(

[R] Superimposing vector polygons over raster grid in a plot

Hello: I would like to superimpose vector polygons (state outlines) from a Shape file on top of a satellite image, imported into a SpatialGridDataFrame from GEOTIFF via gdal_translate and readGDAL. When I plot polygon and point shape files in R, into SpatialPointDataFrame and SpatialPolygonD

Re: [R] Plotmath issue superscript "-"

Hallo, > If you need a subscript as well, I like > > plot(0, main=quote({NO^'\x96'}[3])) I tried this but I get: > plot(0, main=quote({NO^'\x96'}[3])) Errore in title(...) : stringa multibyte non valida ('invalid multibyte string') My R version is: platform i686-redhat-linux-gnu versi

Re: [R] Cutting & pasting help examples into script window

Does this look like a bug? If so, is there a different way to report it? Thanks, Bob > -Original Message- > From: Dirk Eddelbuettel [mailto:[EMAIL PROTECTED] > Sent: Thursday, September 20, 2007 10:17 AM > To: Muenchen, Robert A (Bob) > Subject: Re: [R] Cutting & pasting help examples into

[R] Time series graphs

I'm fairly new to S-Plus and I need to get this done quickly. Suppose I have the following fake data below: There are two companies, call them Bob and Tom. Each have two variables, call them A and B, that have observations. Bob

Re: [R] SEM - singularity error

Dear Nicolette and Chuck, I apologize for not seeing Nicolette's original posting. The model does appear to be correctly specified (though I didn't see the path diagram) and over-identified. I assume that all of the variables with lower-case names are observed variables and that "Moist" and "Hab"

[R] Non-metric multidimensional scaling

Hello everyone, I'm working with R 2.4.1 on a PC running with XP. Trying to run isoMDS as follows: Gquad.mat <- Gquads[4:10] # extracts only the metric data variables Gquad.dist <- dist(Gquad.mat) Gquad.mds <- isoMDS(Gquad.dist) plot(Gquad.mds$points, type = "n") text(Gquad.mds$p

Re: [R] Strange behaviour of lars method

Thanks for your answer! On 20.09.2007, at 12:48, Tim Hesterberg wrote: > Might you be having numerical problems? That might very well be what's happening. Did hope that it only happens late in the procedure so I can ignore it, but those Cp values indicate I maybe can't. > If you're using th

Re: [R] Time series graphs

Using plot.zoo in the zoo package try this: Lines <- "Bob.A Bob.BTom.ATom.B Jan 84 9 8 Feb 7 5 4 7 Mar 6 8 4 4 Apr 3 7

[R] Conditional Logit and Mixed Logit

Hello, Could anybody provide me with codes (procedure) how to obtain Conditional Logit (McFadden) and Mixed Logit (say, assuming normal distribution) estimates in R? Thanks, David U. -- View this message in context: http://www.nabble.com/Conditional-Logit-and-Mixed-Logit-tf4489238.html#a128029

Re: [R] Ambiguities in vector

Birgit Lemcke wrote: > > > Perhaps you haven´t understood my question in the mail yesterday. So > I will try to describe my problem in a different way > > You see the tables. I would like to test the variables between the > tables. > I'm afraid that even before we start to deal with th

[R] convert data from

Hi, I am trying to convert a data frame from: X1X2 A12 B34 C56 to: A01 1 A02 2 B013 B024 C015 C026. How can I do it in R? I appreciate your help. Zheng - [[alternative HTML version deleted]] ___

Re: [R] how can I attach a variable stored in

I don't think I understand your question but John Fox has written a very nice documentat about scoping and environments on his website. It's probably easy to find the site by googling "John Fox" but, if you can't find it, let me know. As I said, I don't think that I understand your question but,

Re: [R] Non-metric multidimensional scaling

>Error in isoMDS(Gquad.dist) : zero or negative distance between objects >179 and 180. >How can I handle this, is it valid to add 0.5 to every element in the >distance matrix or is someother alternative more appropriate? This means that your objects 179 and 180 are identical so just remove one

Re: [R] random Forests

1) See the help pages for save() and load(). We do this all the time. 2) Prediction (or what you call scoring) given the matrix structure is quite simple. Take a look at the underlying C source code. You can write that in any language you want without much problem. Andy From: David Montgomer

Re: [R] SEM - singularity error

Dear John, Thank you so much for your assistance, it is greatly appreciated. You are correct in your interpretation of the variables (Moist and Hab are latent, the other, lower case variables are observed). I've attempted to run the model again adding Moist <-> Hab (please see code below) and sti

Re: [R] Ambiguities in vector

First thanks for your answer. Now I try to explain better: I have species in the rows and morphological attributes in the columns coded by numbers (qualitative variables; nominal and ordinal). In one table for the male plants of every species and in the other table for the female plants of eve

Re: [R] Time series graphs

In pure R (without other packages) is IMHO simpler to understand (at least if data are so simple: months 1 till 9 with no missing values...). All you need is: dummybob <- ' month A B Jan 8 4 Feb 7 5 Mar 6 8 Apr 3 7 May 5 1 Jun 6 4 Jul 2 8 Aug 1 2 Sep 4 3 ' bob <- read.t

Re: [R] Superimposing vector polygons over raster grid in a plot

On Thursday 20 September 2007, Rick Reeves wrote: > Hello: > > I would like to superimpose vector polygons (state outlines) from a > Shape file on top of a satellite image, > imported into a SpatialGridDataFrame from GEOTIFF via gdal_translate and > readGDAL. > > When I plot polygon and point shape

Re: [R] Plotmath issue superscript "-"

Yes, sorry, I should have said that I was on Windows. In a UTF-8 locale, you could try \u2013 in place of \x96. The character is an endash. Peter Ehlers Scionforbai wrote: > Hallo, > >> If you need a subscript as well, I like >> >> plot(0, main=quote({NO^'\x96'}[3])) > > > I tried this but I

Re: [R] how can I attach a variable stored in

Hi Mark - Thanks for the reply. Sorry I didn't really clarify too well what I'm trying to do. The issue is not that I can't see the variable that gets loaded. The issue is that the variable is a list variable, and I'd like to write a function that will take the .RData filename and attach th

Re: [R] Conditional Logit and Mixed Logit

On Thu, 20 Sep 2007, David U. wrote: > > Hello, > > Could anybody provide me with codes (procedure) how to obtain Conditional > Logit (McFadden) and Mixed Logit (say, assuming normal distribution) > estimates in R? RSiteSearch("conditional logit") as the posting guide tells you! > > Th

[R] error in manova

I work with Windows, R 2.4.1. I'm a beginner with R! After doing a Discriminant Function Analysis, I am trying to run manova to get a measure of significance of my lda results. I want to predict groups 1 through 4 using 78 variables (bad group/var ratio, I know, but I'm just exploring the possib

Re: [R] how can I attach a variable stored in

try: attach(get(ls(pat="^foo"))) On 9/20/07, Peter Waltman <[EMAIL PROTECTED]> wrote: > Hi Mark - > > Thanks for the reply. Sorry I didn't really clarify too well what I'm > trying to do. The issue is not that I can't see the variable that gets > loaded. > > The issue is that the variable is a

Re: [R] convert data from

This should get you close to what you want: > x <- read.table(textConnection(" X1X2 + A12 + B34 + C56"), header=TRUE) > # add rownames to the dataframe > x$name <- row.names(x) > require(reshape) # use reshape package [1] TRUE > melt(x, id='name') name variable val

Re: [R] Identify and plotting symbols. [SEC=UNCLASSIFIED]

On 20/09/2007, at 4:35 PM, Crombie, Joe wrote: > Or maybe: > >> while(length(ind <- identify(x,y,n = 1, plot = F))) >> points(x[ind], y[ind], pch = 19) > > (highlights each point as you select it, until you click _stop_) Bewdy!!! That's eggs-actly what I was looking for. Thanks very much.

Re: [R] help with making a function of scatter plot with multiple variables

The simple way is to enclose it in a 'function' and pass parameters. Assuming that you have the same number of parameters, then the following will do: my.func <- function(x,y,d1,v1,s1,t1,s2,t2,s3,t3,s4,t4,s5,t5) { op <- par(bg = "grey97") par(mfrow=c(1,2)) plot(d1,v1, pch="v", col="ora

Re: [R] how can I attach a variable stored in

Here's a function that does what I think you want to do: > attach.firstvar <- function(file) { + tmpenv <- new.env() + vars <- load(file, envir=tmpenv) + x <- get(vars[1], envir=tmpenv, inherits=FALSE) + if (is.list(x)) + attach(x, name=vars[1]) + return(vars) + } > x

Re: [R] Robust or Sandwich estimates in lmer2

Abdus Sattar <[EMAIL PROTECTED]> asked: > > I am trying to find the robust (or sandwich) estimates of the standard > error of fixed effects parameter estimates using the package "lmer2". > Others have already pointed out that this is not implemented in lmer2. You could try a delete-n jackknife, w

Re: [R] Superimposing vector polygons over raster grid in a plot

On Thursday 20 September 2007, Rick Reeves wrote: > Thanks for this, Dylan. Great example. One question: 'elev.pred' and > 'elev.var' are columns in the result 'p' from predict? > thanks, RR Hi, elev.pred and elev.var are the 'zcol' attribute, or rather the value at each cell within the spatial

Re: [R] Time series graphs, question about using zoo

Hi, Can you tell me what is the meaning for "tail, 1" in "aggregate"? I also want to get some similar graph, but the data is not time series data. Suppose here is my data one, I want a graph with x-axis is just the index(1:9). The graph plot all the variable A, B,C,D. So there should be 4 lines f

Re: [R] SEM - singularity error

Dear Nicki, > -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > Sent: Thursday, September 20, 2007 3:06 PM > To: John Fox > Cc: Chuck Cleland; [EMAIL PROTECTED] > Subject: Re: [R] SEM - singularity error > > Dear John, > > Thank you so much for your assistance, it

Re: [R] Cutting & pasting help examples into script window

On 20/09/2007 1:49 PM, Muenchen, Robert A (Bob) wrote: > Does this look like a bug? If so, is there a different way to report it? It sounds like a bug, but I can't reproduce it. You said it is intermittent on your system. Can you try to work out the conditions that reliably trigger it? It mig

Re: [R] Time series graphs, question about using zoo

Try this: Lines <- "A BCD 8 4 9 8 7 5 4 7 6 8 4 4 3 7 6 2 5 1 8 5 6 4 7

Re: [R] Time series graphs, question about using zoo [SEC=UNCLASSIFIED]

Hi Fang, An easy way of doing this is by: > matplot(one, pch = LETTERS[1:4], type = 'b') Cheers Joe Joe Crombie Information and Risk Sciences Bureau of Rural Science Canberra Australia p: +61 2 6272 5906 e: [EMAIL PROTECTED] -Original Message- From: [EMAIL PROTECTED] [mailto:[

[R] Help create a loopto conduct multiple pairwise operations

#Hello, #I have three data frames, X,Y and Z with two columns each and different numbers of rows. # creation of data frame X X.alleles <- c(1,5,6,7,8) X.Freq<- c(0.35, 0.15, 0.05 , 0.10, 0.35) Loc1 <- cbind( X.alleles,X.Freq) X <- data.frame(Loc1) #creation

Re: [R] Help create a loopto conduct multiple pairwise operations

You can put the dataframes in a list and then use 'combn' to determine the possible pairwise combinations and then use this in an lapply to compute the products which is also in a list: > > X.alleles <- c(1,5,6,7,8) > X.Freq<- c(0.35, 0.15, 0.05 , 0.10, 0.35) > Loc1 <- cbind(

Re: [R] Cutting & pasting help examples into script window

Now I'm working in 2.5.1 on a home machine also running XP. It has the same problem, and I think I finally figured it out. I've noticed that if the cursor is directly over the text, it becomes an I-beam. When hovering over the blank space around the text, the cursor becomes an arrow. Selections v

[R] BRugs package question

Hi there, I installed the latest OpenBUGS version (3.0.3) in program files folder but it seems that bugs use its own version of OpenBUGS. here is part of the message bugs returns . Welcome to BRugs running on OpenBUGS version 2.2.0 beta how can i configure BRugs to use the mi

[R] Line Graph - Greater than 2 variables on plot

Hello all, I was wondering if anyone knew how to construct a multiple line graph on R, where there are 2 (or more) sets of data points plotted against some x axis of data, and you can draw a line on the graph connecting each set of data points. For example: A B C

Re: [R] Line Graph - Greater than 2 variables on plot

On 21/09/2007, at 1:57 PM, Wayne Aldo Gavioli wrote: > > > Hello all, > > I was wondering if anyone knew how to construct a multiple line > graph on R, > where there are 2 (or more) sets of data points plotted against > some x axis of > data, and you can draw a line on the graph connecting ea

[R] RJDBC connection help needed

Hi everyone, I'm obviously missing something simple here... Trying to connect to an external db with RJDBC. (I can connect OK with ROLAP and with other java apps, eg dbVisualizer) JDBC call seems to work ok: > driverClass <- "net.sourceforge.jtds.jdbc.Driver" > classPath <- "C:\\Dwns\\jtds\\jtd

Re: [R] Line Graph - Greater than 2 variables on plot

This should do it for you: > x <- read.table(textConnection("A B C D + 0.65662.11851.23205 + 0.647 2.08651.232510 + 0.65322.10601.228715 + 0.64872.12901.231320 + 0.65

[R] Estimate correlation with bootstrap

Hello. I would like to estimate the correlation coefficient from two samples with Bootstrapping using the R-function sample(). The problem is, that I have to sample pairwise. For example if I have got two time series and I draw from the first series the value from 1912 I need the value from 1912