[R] negative axis values in image() and scalebar in cor.plot() {psych}

Dear List, I sent a related message yesterday, but did not receive it through the mailing list; on the R-help archives it reads as "An embedded and charset-unspecified text was scrubbed...". So here it is again, with the (little) progress I have made since then. Any help would be greatly appreciat

Re: [R] Using image/contour with unevenly spaced data...

On Nov 18, 2010, at 10:52 PM, Jason Rupert wrote: Is it possible to plot unevenly spaced data with image/contour function? Below is an example of the type uneven data that I'm trying to plot with image/contour functions. For example, I would like to have the x- distance on the x-axis and

[R] Question on overdispersion

I have a few questions relating to overdispersion in a sex ratio data set that I am working with (note that I already have an analysis with GLMMs for fixed effects, this is just to estimate dispersion). The response variable is binomial because nestlings can only be male or female. I have samples

[R] glmnet_1.5.1 uploaded to CRAN

In glmnet_1.5 a poor default was set for the argument type which caused the program to be very slow or even crash when nvar (p) is very large. The argument type (now called type.gaussian) has two options, "covariance" or "naive", and is used for the default family="gaussion" model (squared erro

[R] Using image/contour with unevenly spaced data...

Is it possible to plot unevenly spaced data with image/contour function? Below is an example of the type uneven data that I'm trying to plot with image/contour functions. For example, I would like to have the x-distance on the x-axis and y-distance on the y-axis and then the temperature value

Re: [R] New Sampling question

I spent the whole afternoon on it, but there is still no progress. I wish I could take some courses... :( -- View this message in context: http://r.789695.n4.nabble.com/New-Sampling-question-tp3047885p3049796.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] About R and RKward

Hi Ajay, Tks for your URL > have you researched out all R commander E plugins before settling on R Kward Whether you meant the interfaces on; http://decisionstats.com/2010/10/05/interfaces-to-r/ No I have no urgent need installing RKward, just interested to know how it work. I tested Red-

[R] a rounding option for \Sexpr{}?

Hi all, Currently Sweave does not format the (especially numeric) value in \Sexpr{}, and I often have to round() the numbers explicitly, otherwise I will get more than 10 digits, which is not necessary for me in most cases. Is there a way to specify the number of digits to be kept for values in \S

Re: [R] conditional mean between two data frames with different levels

Hi Alberto, This should do it. 'x' is equivalent to your alfa vector. Cheers, Josh ## your example data, in a form R can read dat1 <- data.frame(score = c(1, 2, 0, 2, 1, 1, 3, 2, 1, 0), teams = c("a", "b", "c", "d", "e", "a", "b", "c", "d", "e")) dat2 <- data.frame(score = c(2, 3, 1, 0, 0,

[R] latex tables with quantreg

Hi all, I'm using the quantreg package to run some regressions. I'm trying to save the output as latex tables using the helpful "latex" command, which works fine. However, I've been unable to modify its behavior in some ways I'm insterested in. For instance, I'd like a table in which I had the ast

Re: [R] Bold plot axis

On 11/18/2010 4:27 PM, David Winsemius wrote: Well the ticks got wider on the quartz device (Mac) when I used axis(side=1,lwd=3), so only the box would seem to be faint. ?box box(lwd=3) Box was a new one for me. After removing a conflicting option, it seems to have worked. Thank you. J

[R] how do I build panel data/longitudinal data models with AR terms using the plm package or any other package

Hi All, I am doing econometric modeling of panel data (fixed effects). We currently use Eviews to do this, but I have discovered a bug in Eviews 7 and am exploring the use of R to build panel data models / longitudinal data models. I looked at the plm package but do not see how I can incorporat

Re: [R] RowSums Question

Thanks Jorge It works great. I solved it by using loop, but i like your way better. Thanks again Cameron -- View this message in context: http://r.789695.n4.nabble.com/RowSums-Question-tp3049261p3049682.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Problem in installing and starting Rattle

On Mon, Nov 15, 2010 at 10:57 AM, Graham Williams < graham.willi...@anu.edu.au> wrote: > On 16 November 2010 02:40, Feng Mai wrote: > > > > > I also have the problem trying to start rattle > > > > Windows 7 32-bit R 2.12.0 > > > > When I try library(rattle) I get an error message > > "The procedu

Re: [R] Bold plot axis

On Nov 18, 2010, at 5:10 PM, john polo wrote: Dear UseRs, I'd like to make the lines that the axes of my plot run on bold. I know "axis(side=1,lwd=3)" will make 0-max value on the axis bold, but I want the whole thing, from corner to corner (actually all the way around) and tick marks an

Re: [R] Playwith-problem with loading

On Tue, Nov 16, 2010 at 10:59 AM, Fencl, Martin wrote: > Helllo, > I am having trouble with running the library Playwith in the R-2.12.0. > running under 32bit Windows XP. After calling the library the error message > "The procedure entry point gdk_cairo_reset_clip could not be located in the > dy

Re: [R] Sample covariance matrix in R

perhaps ... library("fortunes") fortune("surgery") Cheers, Bert On Thu, Nov 18, 2010 at 1:57 PM, Ben Bolker wrote: > Alaios yahoo.com> writes: > >> >> Hello everyone. >> I would like to find the sample covariance matrix using R. >> >> So far I read on the wikipedia what a sample_covariance is

[R] Bold plot axis

Dear UseRs, I'd like to make the lines that the axes of my plot run on bold. I know "axis(side=1,lwd=3)" will make 0-max value on the axis bold, but I want the whole thing, from corner to corner (actually all the way around) and tick marks and the accompanying values for the ticks bold. In the

Re: [R] rasterImage and coordinate conversion

Hi, For curiosity's sake, below is another version with ggplot2 and Grid graphics, library(pixmap) logo <- read.pnm(system.file("pictures/logo.ppm", package="pixmap")[1]) library(ggExtra) # r-forge, requires gridExtra qplot(rnorm(100),rnorm(100)) + annotate("pixmap", x=-Inf, y=-Inf, picture=lo

Re: [R] Sample covariance matrix in R

Alaios yahoo.com> writes: > > Hello everyone. > I would like to find the sample covariance matrix using R. > > So far I read on the wikipedia what a sample_covariance is > http://en.wikipedia.org/wiki/Sample_covariance > > according to wikipedia one vector is > enough to calculate the sample

Re: [R] randomForest parameters for image classification

Thanks for your response. I agree with you with regards to memory. I am calling predict() that way as per the syntax in http://cran.r-project.org/web/packages/raster/raster.pdf, which has the raster object first, then the model... But in the end it seems to be much more efficient to loop through th

Re: [R] RowSums Question

Hi Cameron, May be this (untested)? rowSums(is.na(tsObj), na.rm = TRUE) HTH, Jorge On Thu, Nov 18, 2010 at 2:38 PM, cameron <> wrote: > > thanks Henrique > > I have another question. > > Lets say i have a timeSeries table >AB C > 1/1/90 NA 1 2 > 1/2/90

Re: [R] Looking up the directory a file is located in

The dirname() function looks very helpful; I hadn't heard of that one before. I'm still reading up to see how the parent.frame command does what it does; I didn't realize that the ofile variable might not be around in the future. Richard's suggestion -- source("c:/myfullpath/myfile.r", chdir=T

[R] Kalman filter

Hello, I have completed my kalman filter problem with more details. The transition- and the measurement equation is given by x[t]=A[t]*x[t-1]+B[t]*epsilon[t] y[t]=C[t]*x[t]+eta[t] A, y, B and C are Matrices. Y[t] is the data input vector with 800 elements (every t has one element) M

Re: [R] RowSums Question

thanks Henrique I have another question. Lets say i have a timeSeries table AB C 1/1/90 NA 1 2 1/2/90 NA 1 1 1/3/90 NA 1 -1 1/4/90 NA -1 1 1/5/901 1 1 1/6/901 51 1 I want to

Re: [R] rasterImage and coordinate conversion

Hi Rajarshi Guha wrote: Hi, I have a plot and I would like to overlay a PNG image over it. I'm using the rasterImage function to do this, but the problem I'm facing is working out the coordinates of the upper right corner of the final image in user coordinates. That is I can place the image so

Re: [R] dmultinomial

Alexander Shenkin ufl.edu> writes: > > Hello All, > > I'm trying to run a maximum likelihood analysis using dmultinomial (i'm > avoiding dmultinom as I'd like to run it with vectors for the ML stuff). > However, I'm having a hard time getting even the simplest example > running. Any help woul

Re: [R] sweep by levels of a factor

On Nov 18, 2010, at 1:23 PM, Henrique Dallazuanna wrote: Try this also: x$x - ave(x$x, x$f) Could also use scale for only its centering capacities: unlist( tapply(fx$x, fx$f, scale, scale=FALSE) # default for center == TRUE ) 11 12 21 22 23 1 -1 -2 2 0 ... but looking at the c

Re: [R] odfWeave/XML Windows issue

That seemed to do the trick, Thanks Ales!!! -- View this message in context: http://r.789695.n4.nabble.com/odfWeave-XML-Windows-issue-tp3036376p3049524.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] how exactly does 'identify' work?

Omg! yes, it is working now! tmp <- qqnorm( resid(test.lm) ) What a simple nice trick!!! Actually, i wasnt looking for the 'i'th label, I was looking for the 'row.names' as label, like I stated in the 1st post. > identify(tmp, , row.names(test) ) is the label i have been trying to get. T

Re: [R] how exactly does 'identify' work?

One additional point on your original post. You added row names to the test data frame, but did not specify the name of the data frame when you did the regression, rather you attached the data frame. When you did this lm found x and y, but did not find the rownames, so the diagnostic plot just

Re: [R] how exactly does 'identify' work?

Did you read the help page for qqnorm? The return value has the x and y coordinates used, you can just do something like: > tmp <- qqnorm( resid(test.lm) ) > identify(tmp, , names(resid(test.lm)) ) Or the plot.lm function has an argument id.n that automatically labels the n most extreme values

Re: [R] generalized linear regression - function glm - dismissed predictors - more information about simulated data

Bert Gunter gene.com> writes: [snip snip snip] > Cheers, > Bert > > P.S. Your results look fishy, even for a simulation. The extremely low > standard errors and rank deficiencies strongly suggest to me that your > methodology is faulty and that you are getting a bunch of baloney. I > would reco

Re: [R] generalized linear regression - function glm - dismissed predictors - more information about simulated data

... but do note the first line of the DESCRIPTION in the Help file: "Compute the rank of matrix, a well-defined functional in theory, somewhat ambigous in practice." So as the Help file makes very clear, what you get can depend on how you do it; and may be different than what you want; which will

[R] dmultinomial

Hello All, I'm trying to run a maximum likelihood analysis using dmultinomial (i'm avoiding dmultinom as I'd like to run it with vectors for the ML stuff). However, I'm having a hard time getting even the simplest example running. Any help would be greatly appreciated. > library(mc2d) > dmultin

Re: [R] generalized linear regression - function glm - dismissed predictors - more information about simulated data

On Nov 18, 2010, at 1:31 PM, Christine SINOQUET wrote: Hello, a) Thanks a lot for the answer relative to the dismissed coefficients. However, the result below has been obtained in the R console and there was no warning (please, see the full display below (***) if you wish). I guess my

[R] R Courses***Fundamentals & Advanced Programming ***December 2010 - in San Francisco, New York City and Washington DC by XLSolutions Corporation

Check out our Fundamentals and Advanced Programming R courses for December 2010 in Washington DC, San Francisco and New York City http://www.xlsolutions-corp.com/Rcourses *** R Fundamentals and Programming Techniques New York, Dec 13-14, 2010 Washington DC, Dec 16-17, 2010 San

Re: [R] how exactly does 'identify' work?

yes, i tried to modify the "2L" part in plot.lm ### if (show[2L]) { ylim <- range(rs, na.rm = TRUE) ylim[2L] <- ylim[2L] + diff(ylim) * 0.075 qq <- qqnorm(rs, main = main, ylab = ylab23, ylim = ylim, ...) if (qqline)

Re: [R] how exactly does 'identify' work?

On 18/11/2010 1:50 PM, casperyc wrote: Hi, I think the problem is 1 - when a linear model is fitted, ploting the qqnorm( test.lm$ res ) we dont 'know' what values are actually being used on the y-axis; and how do we refer to the ‘Index’ on the x-axis?? therefore, i dont know how to refer

Re: [R] odfWeave/XML Windows issue

Dave_F battelle.org> writes: > > > I am getting the following error when using odfWeave > > Error in xmlEventParse(infile, handlers = handlers, trim = FALSE, state = > state) : > File content_1.xml does not exist > > > Thanks, > Dave Dear Dave! I was getting the same message while tryi

[R] generalized linear regression - function glm - dismissed predictors - more information about simulated data

Hello, a) Thanks a lot for the answer relative to the dismissed coefficients. However, the result below has been obtained in the R console and there was no warning (please, see the full display below (***) if you wish). "yi ~ X$V1 + X$V2 + X$V3 + X$V4 + X$V5 + X$V6 + X$V7 + X$V8 + X$V9 + X$V

Re: [R] how exactly does 'identify' work?

Hi, I think the problem is 1 - when a linear model is fitted, ploting the qqnorm( test.lm$ res ) we dont 'know' what values are actually being used on the y-axis; and how do we refer to the ‘Index’ on the x-axis?? therefore, i dont know how to refer to the x and y coordinates in the identi

Re: [R] predict() an rpart() model: how to ignore missing levels in a factor

I don't think that, considering the mechanism behind recursive partitioning, that there is any way for you to ignore the crop factor if it is not in the original test set. What decision should be made if, for instance, the next split in a decision tree were on crops and output was 5 for apples,

Re: [R] Looking up the directory a file is located in

On 18/11/2010 11:56 AM, Cliff Clive wrote: Thanks, Gabor! So far I like this one best: https://stat.ethz.ch/pipermail/r-help/2005-November/082347.html So if my script is called "myRscript.r", I can do the following: this.file = parent.frame(2)$ofile this.dir = gsub("/myRscript.r", "", this.fi

Re: [R] Sample covariance matrix in R

The help and its examples are very comprehensive here. The usage you cite shows exactly what you need to do From: Alaios [mailto:ala...@yahoo.com] Sent: Thursday, November 18, 2010 1:30 PM To: Doran, Harold Subject: RE: [R] Sample covariance matrix in R Checked that Usage covr(x, y = NULL, na.r

Re: [R] RowSums Question

Try this: rowSums(tsObj, na.rm = TRUE) On Thu, Nov 18, 2010 at 3:58 PM, cameron wrote: > > > I have a question on RowSums. > > Lets say i have a timeSeries table >A B C > 1/1/90 NA 1 1 > 1/2/90 NA 1 1 > 1/3/90 NA 1 1 > 1/4/90 NA 1

Re: [R] Sample covariance matrix in R

Alex: ?cov > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of Alaios > Sent: Thursday, November 18, 2010 12:54 PM > To: Rhelp > Subject: [R] Sample covariance matrix in R > > Hello everyone. > I would like to find the sample cov

Re: [R] sweep by levels of a factor

Try this also: x$x - ave(x$x, x$f) On Thu, Nov 18, 2010 at 12:02 PM, Lancaster, Anthony < anthony_lancas...@brown.edu> wrote: > Hi, > I'd appreciate help with this. I have a data matrix with one column, called > f in the example below, a factor. I'd like to subtract the means from each > of >

[R] RowSums Question

I have a question on RowSums. Lets say i have a timeSeries table A B C 1/1/90 NA 1 1 1/2/90 NA 1 1 1/3/90 NA 1 1 1/4/90 NA 1 1 1/5/901 1 1 1/6/901 1 1 if i use RowSums, i will get 1/5/903 1/

Re: [R] Looking up the directory a file is located in

Thanks, Gabor! So far I like this one best: https://stat.ethz.ch/pipermail/r-help/2005-November/082347.html So if my script is called "myRscript.r", I can do the following: this.file = parent.frame(2)$ofile this.dir = gsub("/myRscript.r", "", this.file) setwd(this.dir) This will set the work

Re: [R] New Sampling question

Hi, Mike, thank you very much!!:) the following two functions are really helpful, which I didn't even know. Actually, I searched the forum for something like this, but failed. Now I am still trying to make up my own functions. :) sample(a[,2],2) #randomly draw two numbers from a column sample(

Re: [R] sweep by levels of a factor

Thank you very much, Tony L. On 18 November 2010 14:02, Lancaster, Anthony wrote: > Hi, > I'd appreciate help with this. I have a data matrix with one column, called > f in the example below, a factor. I'd like to subtract the means from each > of > other columns for each level of the factor.

Re: [R] conditional mean between two data frames with different levels

Thank you very much Josh, I guess you`re right. So this is an example: data frame 1 has 2 columns and 10 rows. The first column is "score" a variable indicating the number of goals scored by a football team score<-c(1,2,0,2,1,1,3,2,1,0), column 2 contains the football "teams " where teams<-c(a,b,c

Re: [R] Logistic regression with factorial effect

You would be better off posting to R-sig-mixed-models or R-sig-ecology -- Bert On Thu, Nov 18, 2010 at 9:32 AM, Billy.Requena wrote: > > Hello, > > I’d like to evaluate the temporal effect on the relationship between a > continuous variable (e.g. size) and the probability of mate success. > Init

Re: [R] Help with lmer, nested data and repeated measures

You would be better off posting to R-sig-mixed-models or R-sig-ecology -- Bert On Thu, Nov 18, 2010 at 9:58 AM, Simon Garnier wrote: > Dear all, > > I'm discovering the somehow confusing (for me) world of linear mixed models > after having been advised it could be the best option to analyze my

[R] Help with lmer, nested data and repeated measures

Dear all, I'm discovering the somehow confusing (for me) world of linear mixed models after having been advised it could be the best option to analyze my dataset. After several days of reading, I'm not sure that what I ended up with makes some sense and I'd greatly appreciate any help and expla

[R] Sample covariance matrix in R

Hello everyone. I would like to find the sample covariance matrix using R. So far I read on the wikipedia what a sample_covariance is http://en.wikipedia.org/wiki/Sample_covariance according to wikipedia one vector is enough to calculate the sample covariance matrix. In R I tried cov(myvector) a

[R] predict() an rpart() model: how to ignore missing levels in a factor

I am using an algorigm to split my data set into two random sections repeatedly and constuct a model using rpart() on one, test on the other and average out the results. One of my variables is a factor(crop) where each crop type has a code. Some crop types occur infrequently or singly. when the d

[R] Logistic regression with factorial effect

Hello, I’d like to evaluate the temporal effect on the relationship between a continuous variable (e.g. size) and the probability of mate success. Initially I was trying to do a logistic regression model incorporating the temporal effect, but I don’t know if that is the best option. I simulated s

[R] lme Random Effects and Covariates

1. I'm attempting to test for Random Effects. I've grouped the data on subject (grid) but want to use lme to build the model without subject as a RE then add it and do anova between the 2 models. This is the result I get and it appears it's adding Random Effects. tmp.dat4 <- groupedData(Trials

Re: [R] problems subsetting

On Nov 18, 2010, at 10:42 AM, David Winsemius wrote: On Nov 18, 2010, at 10:25 AM, Martin Tomko wrote: Hi Gerrit, indeed, that works. Excellent tip! For reference, I did this: subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) I am still not totally sure when one uses "&" am

[R] dmultinomial

Hello All, I'm trying to run a maximum likelihood analysis using dmultinomial (i'm avoiding dmultinom as I'd like to run it with vectors for the ML stuff). However, I'm having a hard time getting even the simplest example running. Any help would be greatly appreciated. > library(mc2d) > dmultin

Re: [R] kalman filter in sspir

To learn why sspir does not have a filter function, you need to ask the package maintainer, Claus Dethlefsen . My belief is that he, Soren Lundbye-Christensen and Anette Luther Christensen found other outlets for their time since they completed the package and the companion paper, "C. D

Re: [R] problems subsetting

Hi Martin, I think Steve meant dput(df), as I already told you Ivan Le 11/18/2010 17:40, Steve Lianoglou a écrit : Hi, On Thu, Nov 18, 2010 at 10:42 AM, Martin Tomko wrote: Thanks Steve, that explains it... Unfortunately, I did nto get that from my R docs... I am still searching for a sol

Re: [R] how exactly does 'identify' work?

I think that your problem comes from a misunderstanding. The general rule is that you give the plot command 2 vectors, x and y (though you can give it the vectors separately, or together in a list or matrix). If you give plot only a single vector then it will use this as the y vector and use t

Re: [R] stacking consecutive columns

On Wed, Nov 17, 2010 at 10:37 AM, Graves, Gregory wrote: > Follows is the exact solution to this: > > v <- NULL > > #note that decreasing is FALSE so preceding year preceeds > > for(i in 2:46) { >  kk <- melt(yearmonth[, c(1, i, i+1)], id.vars="month", variable_name="year") >  v[[i-1]] <- kk[order

Re: [R] conditional mean between two data frames with different levels

Hi Alberto, It would help if you could provide a small example. I might break the problem down into three parts: 1) create a vector that has the final subgroupings you want 2) find the conditional means by subgroup 3) replicate the means as needed. My first guess would be start with: "==" or

Re: [R] New Sampling question

You could try writing a loop a<-data.frame(c(1:10),c(21:30)) M<-10 #number of iterations- scale up to 1000 once you get your sampling function working res<-NULL #place to store your results for i in (1:M) { ares<-sample(a[,2],1) res<-c(res, ares) } res It's up to you how to

Re: [R] problems subsetting

Hi, On Thu, Nov 18, 2010 at 10:42 AM, Martin Tomko wrote: > Thanks Steve, that explains it... Unfortunately, I did nto get that from my > R docs... > > I am still searching for a solution for matching the entries in my matrix by > matching the rownames to the entires in a subset I got using Gerri

Re: [R] genralized linear regression - function glm - number of

On Nov 18, 2010, at 11:00 AM, Christine SINOQUET wrote: Hello, Performing a linear regression through the function glm ("yi ~ X$V1 + X$V2 + X$V3 + X$V4 + X$V5 + X$V6 + X$V7 + X$V8 + X$V9 + X$V10"), I then edit the information about the coefficients: print(coefficients(summary(fit))) I n

Re: [R] Updata Rdata File

Hi Jason, You can add as many objects as you want to save(). For instance: save(data1, data2, data3, data4, file = "whatever"). Once you have loaded the data, just add whatever data you want, and then save it (or edit your original object and resave it). Here is another example: > ls() charact

Re: [R] Updata Rdata File

?rbind "Jason Kwok" wrote: >Thanks for the response Duncan. > >I'm sorry I wasn't clear in my question. I need to add an additional >value >every day to the data and not update it. > >I know how to load and save but I dont' know how to add data. >load(file="C:\\datafile.rdata") >save(data,file=

Re: [R] how to find near neighbors?

Have a look at nncross() from the spatstat package. Best regards, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek team Biometrie & Kwaliteitszorg Gaverstraat 4 9500 Geraardsbergen Belgium Research

Re: [R] Odp: Sampling problem

I tried. yours works too. thanks a bounch!! Man -- View this message in context: http://r.789695.n4.nabble.com/Sampling-problem-tp3043804p3049013.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https:

Re: [R] New Sampling question

Also, I need some function at the end which would enable me to draw 1000 such random samples. thanks! :) -- View this message in context: http://r.789695.n4.nabble.com/New-Sampling-question-tp3047885p3048958.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] New Sampling question

Hi Wallace, Have you tried playing with sample()? Note that you can apply this function both to whole dataframes, as well as specific items within a vector. If you play with applying the function to different ways of indexing your sample data, you will likely arrive at your solution. for example:

[R] genralized linear regression - function glm - number of

Hello, Performing a linear regression through the function glm ("yi ~ X$V1 + X$V2 + X$V3 + X$V4 + X$V5 + X$V6 + X$V7 + X$V8 + X$V9 + X$V10"), I then edit the information about the coefficients: print(coefficients(summary(fit))) I note that the number of coefficients (7) is lower than the num

[R] conditional mean between two data frames with different levels

Hi guys, I have two data frames: one referred to 2008 and one to 2009. Their structure is identical except for the different data in them. I need to create a vector alfa of the same length of the dataframe 2009 and fill each element with the mean of 2008$var1 conditional to the subgroup indicated

Re: [R] aspect ratio 1 and blank space

Look at the squishplot function in the TeachingDemos package, it may do what you want, or if not quite then you could perhaps tweak the code to include the values you want and create the correct aspect ratio. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.

[R] how to find near neighbors?

I am looking for an efficient way to find near neighbors... More specifically... I have two sets of points: A & B and I want to find points in set B which are closer to set A than some cutoff (or n-closest) I will appreciate very much any pointers... Ryszard

Re: [R] Looking up the directory a file is located in

On Wed, Nov 17, 2010 at 4:08 PM, Cliff Clive wrote: > > Hello everyone, > > This should be an easy question, I think. > > I'd like to write a command in a program to set the working directory to > whatever directory the file is currently stored in.  Suppose I have a file > called "myRscript.r", an

Re: [R] Looking up the directory a file is located in

On Thu, Nov 18, 2010 at 3:26 PM, Cliff Clive wrote: > > So it sounds like the best we can do in R is to keep track of the script in a > sort of master file that runs the script, and set the working directory in > the master.  Is that accurate? Errr no. Maybe. What? Your script is somewhere. So

Re: [R] Looking up the directory a file is located in

On Thu, Nov 18, 2010 at 10:26 AM, Cliff Clive wrote: > In Python any time you run a script, there is a built-in "__file__" > variable > that can tell you the file name of the script itself. It would be nice to > have a feature like this in R. > R has that feature. source("c:/myfullpath/myfile.r

Re: [R] problems subsetting

On Nov 18, 2010, at 10:42 AM, Martin Tomko wrote: Thanks Steve, that explains it... Unfortunately, I did nto get that from my R docs... I am still searching for a solution for matching the entries in my matrix by matching the rownames to the entires in a subset I got using Gerrit's metho

Re: [R] Help me with R plotting

Thanks a lot :) nice one --- On Thu, 11/18/10, baptiste Auguié wrote: From: baptiste Auguié Subject: Re: [R] Help me with R plotting To: "Alaios" Cc: "Rhelp" Date: Thursday, November 18, 2010, 3:34 PM Hi, Try this, robots <- data.frame(id=letters[1:20],                      x=rnorm(20),

Re: [R] New Sampling question

Dear Ista Zahn-2, If you can give me some advice, I really appreciate it. I have been working on it for days. it seems hard for some novice of R like me to write flexible functions myself. This is for my dissertation. CSE and WSE are two scales of the same construct. The sampling strategy I wan

Re: [R] Looking up the directory a file is located in

So it sounds like the best we can do in R is to keep track of the script in a sort of master file that runs the script, and set the working directory in the master. Is that accurate? In Python any time you run a script, there is a built-in "__file__" variable that can tell you the file name of t

Re: [R] problems subsetting

On Nov 18, 2010, at 10:25 AM, Martin Tomko wrote: Hi Gerrit, indeed, that works. Excellent tip! For reference, I did this: subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) I am still not totally sure when one uses "&" amd when "&&" - I was under the impression that && stan

Re: [R] problems subsetting

On 2010-11-18 07:25, Martin Tomko wrote: Hi Gerrit, indeed, that works. Excellent tip! For reference, I did this: subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) I am still not totally sure when one uses "&" amd when"&&" - I was under the impression that&& stands for logical

Re: [R] problems subsetting

Thanks Steve, that explains it... Unfortunately, I did nto get that from my R docs... I am still searching for a solution for matching the entries in my matrix by matching the rownames to the entires in a subset I got using Gerrit's method. Any idea? Thanks Martin On 11/18/2010 4:35 PM, Ste

Re: [R] problems subsetting

Hi, On Thu, Nov 18, 2010 at 10:25 AM, Martin Tomko wrote: > Hi Gerrit, > indeed, that works. Excellent tip! > > For reference, I did this: > > subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) > > I am still not totally sure when one uses "&" amd when "&&"  - I was under > the impr

Re: [R] ncdf4 for Windows/R-2.12.0

Hi David, This link looks like it points to the resources you need: http://www.unidata.ucar.edu/software/netcdf/docs/faq.html#windows_netcdf4 Perry -- View this message in context: http://r.789695.n4.nabble.com/ncdf4-for-Windows-R-2-12-0-tp3047807p3048935.html Sent from the R help mailing li

Re: [R] Help me with R plotting

(Assuming you aren't going into ggplot2) You can have a look here, it might offer a nice solution: http://rgm2.lab.nig.ac.jp/RGM2/R_man-2.9.0/library/shape/man/colorlegend.html Also, you can use ?cut to change your coloring to specific chunks (although having a full gradient is probably nicer)

Re: [R] Exporting Dataframe to Excel spreadsheet

On Thu, Nov 18, 2010 at 10:26 AM, pankaj borah wrote: > I have a  data frame contains 30,000 rows and 105 columns.  How can i > extract each of the columns as separate  Excel spreadsheet of the same > excel file. See: http://rwiki.sciviews.org/doku.php?id=tips:data-io:ms_windows -- Statistics

Re: [R] Help me with R plotting

Hi, Try this, robots <- data.frame(id=letters[1:20], x=rnorm(20), y=rnorm(20), consumption=runif(20, 10, 100)) library(ggplot2) ggplot(robots) + geom_point(aes(x, y, colour=cut(consumption, c(0, 30, 50, 100 + geom_text(aes

Re: [R] Updata Rdata File

Thanks for the response Duncan. I'm sorry I wasn't clear in my question. I need to add an additional value every day to the data and not update it. I know how to load and save but I dont' know how to add data. load(file="C:\\datafile.rdata") save(data,file="C:\\datafile.rdata") Thanks. Jason

[R] Exporting Dataframe to Excel spreadsheet

Hi, I have a  data frame contains 30,000 rows and 105 columns.  How can i extract each of the columns as separate  Excel spreadsheet of the same excel file.  Regards, Pankaj Barah Department of Biology, Norwegian University of Science & Technology (NTNU) Realfagbygget, N-7491 Trondheim

Re: [R] problems subsetting

Hi Gerrit, indeed, that works. Excellent tip! For reference, I did this: subset1<-subset(summarystats,(Type==1)&(Class==1)&(Category==1)) I am still not totally sure when one uses "&" amd when "&&" - I was under the impression that && stands for logical AND Thanks a lot. Martin On 11

Re: [R] Help me with R plotting

Thanks a lot for your help I tried first this test<-(myvalues)-mean(myvalues) so to put the values close to the 0...30 scale... then I passed this as an argument plot(x,y,col=test) which plots the places where the robots are with a small color. I think this is ok for now. What I want next is a

[R] Plotting number of patients at risk below survival curve

Dear list, does anyone know of a R-package that has implemented the increasingly popular inclusion of the number of patients at risk below Kaplan-Meier curves like in http://bloodjournal.hematologylibrary.org/content/vol116/issue19/images/large/zh89991058760001.jpeg any hint (or negative answer)

  1   2   >