Re: [R] Rolling Sample VAR

2012-05-26 Thread bantex
But after looking through rollapply I still don't seem to be able to implement it to my problem. Could you make it more explicit for me to understand? B -- View this message in context: http://r.789695.n4.nabble.com/Rolling-Sample-VAR-tp4631328p4631430.html Sent from the R help mailing list

Re: [R] How to draw y axis in horizontal bar plot?

2012-05-26 Thread Manish Gupta
continued to previous mail ... http://r.789695.n4.nabble.com/file/n4631429/Image1.png After plotting y axis it will look like -- View this message in context: http://r.789695.n4.nabble.com/How-to-draw-y-axis-in-horizontal-bar-plot-tp4631428p4631429.html Sent from the R help

Re: [R] How to change width of bar when there are very few bars?

2012-05-26 Thread Manish Gupta
If i assign a vector with same small width for each bar (here) still both are same broad. I think it takes ratio for whole area. Is there any other parameter to control it.? Any working example will help me. -- View this message in context:

[R] How to draw y axis in horizontal bar plot?

2012-05-26 Thread Manish Gupta
http://r.789695.n4.nabble.com/file/n4631428/Image2.png How can i draw y axis in bar plot at x = 0? After plotting axis it will look like. -- View this message in context: http://r.789695.n4.nabble.com/How-to-draw-y-axis-in-horizontal-bar-plot-tp4631428.html Sent from the R help mailing list

Re: [R] How to change width of bar when there are very few bars?

2012-05-26 Thread Jim Lemon
On 05/26/2012 12:26 PM, Manish Gupta wrote: If i assign a vector with same small width for each bar (here) still both are same broad. I think it takes ratio for whole area. Is there any other parameter to control it.? Any working example will help me. Hi Manish, If you are using the barplot

[R] Assessing interaction effects in GLMMs

2012-05-26 Thread Luke Duncan
Dear R gurus I am running a GLMM that looks at whether chimpanzees spend time in shade more than sun (response variable 'y': used cbind() on counts in the sun and shade) based on the time of day (Time) and the availability of shade (Tertile). I've included some random factors too which are the

[R] How to transform OLS covariance matrix to White standard errors?

2012-05-26 Thread Dunken
Hi! I am working with a regression of a log-log model that suffers from heteroskedasticity. I have calculated the White standard errors. I would like to use these White standard errors in a RESET test instead of the originally OLS standard errors calculated by the regression. How can I transform

[R] Values in scatterplot??

2012-05-26 Thread Jhope
Hi R-listers, I am having trouble plotting the values for the R2 line, lowess smooth, smothered conditional spread, outlier identification. I have tried the function text and legend but was unsuccessful. Please see below. Your help would be appreciated. Thanks in advance. Jean library(car)

Re: [R] Java problem - XLConnect/xlsx package

2012-05-26 Thread rhaidinger
Hello Martin, thank you for your reply. Yes, it is maybe an architecture clash, we run win 7 64 bit and java 32 bit. But R is 32-bit and the path variable is set to the 32-bit java folder. So maybe the problem don't lies there. Best regards Rainer -- View this message in context:

[R] how to check given number seq. is time series or not?

2012-05-26 Thread sagarnikam123
i have following numbers 0.889046409368551 1.22726162946495 1.22726162946495 1.35785109728356 1.35785109728356 1.10704609982913 1.4424189950435 1.2277843378837 1.35785109728356 0.970883941918588

Re: [R] File format for single channel analysis of Agilent microarray data in Limma?

2012-05-26 Thread Uwe Ligges
Please ask BioC related questions on the BioC mailing list. Best, Uwe Ligges On 23.05.2012 11:44, Parisa wrote: Hi, I am following the protocol outlined here for analysis of single channel Agilent microarray data:

Re: [R] Multiple cbind according to filename

2012-05-26 Thread Uwe Ligges
On 25.05.2012 16:42, Ken wrote: Matthew Ouellettemouellette89at gmail.com writes: Hi all, I'm just a beginner with R but I have not been able to search for any relevant answer to my problem. I apologize if it has in fact been asked before. Recently I've realized that I need to combine

Re: [R] how a latent state matrix is updated using package R2WinBUGS

2012-05-26 Thread Uwe Ligges
On 23.05.2012 17:20, Jean V Adams wrote: I'm trying to understand how a latent state matrix is updated by the MCMC iterations in a WinBUGS model, using the package R2WinBUGS and an example from Kery and Schaub's (2012) book, Bayesian Population Analysis Using WinBUGS. The example I'm using is

Re: [R] Rolling Sample VAR

2012-05-26 Thread Uwe Ligges
On 26.05.2012 07:34, bantex wrote: But after looking through rollapply I still don't seem to be able to implement it to my problem. Could you make it more explicit for me to understand? PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented,

Re: [R] quotes in Rscript -e through system

2012-05-26 Thread Uwe Ligges
On 24.05.2012 12:58, abhagwat wrote: I figured out how to use quotes and parentheses when using Rscript -e (on a bash shell): Rscript -e write\(1,\a.txt\\) -- Question 1: why do the parentheses need to be escaped in the shell? (More a shell than an R question) Because the shell shopudl not

Re: [R] How to transform OLS covariance matrix to White standard errors?

2012-05-26 Thread David Winsemius
On May 26, 2012, at 3:09 AM, Dunken wrote: Hi! I am working with a regression of a log-log model that suffers from heteroskedasticity. I have calculated the White standard errors. I would like to use these White standard errors in a RESET test instead of the originally OLS standard

Re: [R] how to check given number seq. is time series or not?

2012-05-26 Thread R. Michael Weylandt
Thank you for your list of numbers. Next time I run out I'll know where to find some. Now what exactly is your question? Any series of data can be a time series -- it's just a matter of interpretation. Michael On Sat, May 26, 2012 at 5:13 AM, sagarnikam123 sagarnikam...@gmail.com wrote: i have

Re: [R] How to draw y axis in horizontal bar plot?

2012-05-26 Thread R. Michael Weylandt
I'm not sure there's a direct way to do it, but you could most likely use abline(v = 0) to add one to a plot. E.g., barplot(sample(-5:5), horiz = TRUE); abline(v = 0, lwd = 3, col = 2) Michael On Fri, May 25, 2012 at 10:55 PM, Manish Gupta mandecent.gu...@gmail.com wrote: continued to

Re: [R] shift + arrow in getGraphicsEvent

2012-05-26 Thread Uwe Ligges
On 26.05.2012 00:33, cgenolin wrote: Hi the list, Do you know if getGraphicsEvent can detect the use of keys combinaison like shift+arrow or Crlt+arrow ? Some key combinations, but not those mentioned above. Just write a simple test: getGraphicsEvent(onKeybd = function(key) cat(key, \n))

Re: [R] Values in scatterplot??

2012-05-26 Thread Uwe Ligges
On 26.05.2012 08:17, Jhope wrote: Hi R-listers, I am having trouble plotting the values for the R2 line, lowess smooth, smothered conditional spread, outlier identification. I have tried the function text and legend but was unsuccessful. Please see below. Your help would be appreciated.

Re: [R] Selecting with mouse the lines drawn by matplot()

2012-05-26 Thread Uwe Ligges
On 25.05.2012 22:52, servet cizmeli wrote: Dear all, I have a dataset of spectral measurements of sunlight. One row for each different observation, one column for each spectral channel. I would like to select the spectral curves drawn with matplot() (I only know matplot() as a way of easily

Re: [R] R base package grid does not output raster image

2012-05-26 Thread Uwe Ligges
On 24.05.2012 23:04, Patrick Nicholson wrote: I am running 64-bit R 2.15.0 on a Windows Server 2008 R2 Amazon EC2 instance. grid does not produce output. For example, the following code should print the R logo to the window() device: library(grid) library(png) img.path- system.file(img,

Re: [R] change colors/ symbols of lda plots

2012-05-26 Thread Uwe Ligges
On 25.05.2012 21:50, dga...@huskers.unl.edu wrote: Greetings R experts, I am running a simple lda on some simulation data of mine to show an illustration of my multivariate response data, since it is a simulation I have a very large amount of data and the default on plot seems to plot the

Re: [R] how to check given number seq. is time series or not?

2012-05-26 Thread Rui Barradas
Hello, You can't. A collection of observations does not make a time series, you need an index variable, such as time. Rui Barradas Em 26-05-2012 10:13, sagarnikam123 escreveu: i have following numbers 0.889046409368551 1.22726162946495 1.22726162946495

Re: [R] Kolmogorov-Smirnov test and the plot of max distance between two ecdf curves

2012-05-26 Thread Rui Barradas
Hello, Try the following. (i've changed the color of the first ecdf.) loga - log10(a+1) # do this logb - log10(b+1) # only once f.a - ecdf(loga) f.b - ecdf(logb) # (2) max distance D x - seq(min(loga, logb), max(loga, logb), length.out=length(loga)) x0 - x[which( abs(f.a(x) - f.b(x)) ==

[R] Kolmogorov-Smirnov test and the plot of max distance between two ecdf curves

2012-05-26 Thread maxbre
Hi all, given this example #start a-c(0,70,50,100,70,650,1300,6900,1780,4930,1120,700,190,940, 760,100,300,36270,5610,249680,1760,4040,164890,17230,75140,1870,22380,5890,2430) length(a) b-c(0,0,10,30,50,440,1000,140,70,90,60,60,20,90,180,30,90,

[R] avoid error within for loop, try, trycatch, while, move to next iteration, unlist

2012-05-26 Thread Λεωνίδας Μπαντής
Hi there, I would like to ask something about how to avoid a possible error message within a for loop. I am running a simulation and in some repetitions there may be an error that will cause a crash and stop the whole procedure, what I want is to simply move on to the next iteration

[R] Problem with readHTMLTable

2012-05-26 Thread Sven D
Hello All, i was trying to simply run the readHTMLTable on the example published in the package. And on a page I was working on. So running: u = http://en.wikipedia.org/wiki/List_of_countries_by_population; tables = readHTMLTable(u) returns the following error: Error in tb[[thead]] : subscript

Re: [R] how to check given number seq. is time series or not?

2012-05-26 Thread David Winsemius
On May 26, 2012, at 7:50 AM, R. Michael Weylandt wrote: Thank you for your list of numbers. Next time I run out I'll know where to find some. Now what exactly is your question? Any series of data can be a time series -- it's just a matter of interpretation. After I stopped chuckling at

[R] Plotting interactions from lme with ggplot

2012-05-26 Thread dadrivr
I'm fitting a lme growth curve model with two predictors and their interaction as predictors. The multilevel model is nested so that level 1 is time within the individual, and level 2 is the individual. I would like to plot the mean group-level trajectories at plus and minus 1 SD from the mean of

[R] How to assign height value on bar plot?

2012-05-26 Thread Manish Gupta
Hi, How to assign height value on bar plot for each bar? http://r.789695.n4.nabble.com/file/n4631457/Barplot.png -- View this message in context: http://r.789695.n4.nabble.com/How-to-assign-height-value-on-bar-plot-tp4631457.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] avoid error within for loop, try, trycatch, while, move to next iteration, unlist

2012-05-26 Thread William Dunlap
You could put the try() into a while loop, inside for for or lapply loop, as in: f - function(x){ stopifnot(x0.5) ; x } o - numeric(100); for(i in seq_along(o)) { while(inherits(try(tmp - f(runif(1)), silent=TRUE), try-error)) {} ; o[[i]] - tmp } range(o) [1] 0.5000588 0.9991261 Bill Dunlap

Re: [R] How to assign height value on bar plot?

2012-05-26 Thread Uwe Ligges
On 26.05.2012 18:07, Manish Gupta wrote: Hi, How to assign height value on bar plot for each bar? See ?barplot. Uwe Ligges http://r.789695.n4.nabble.com/file/n4631457/Barplot.png -- View this message in context:

Re: [R] avoid error within for loop, try, trycatch, while, move to next iteration, unlist

2012-05-26 Thread William Dunlap
You could encapsulate this idiom in a function: untilOK - function (expr, silent = TRUE, maxIter = 1000, quotedExpr = substitute(expr), envir = parent.frame()) { while ((maxIter - maxIter - 1) = 0 inherits(tmp - try(eval(quotedExpr, envir = envir), silent = silent),

Re: [R] How to measure level of similarity of two data frames

2012-05-26 Thread Özgür Asar
Hi, You can try rdist function under the package fields, or dist package under stats package. Best Ozgur - Ozgur ASAR Research Assistant Middle East Technical University Department of Statistics 06531, Ankara Turkey Ph: 90-312-2105309

Re: [R] R base package grid does not output raster image

2012-05-26 Thread Patrick Nicholson
Thanks Uwe, I had tried this before but hearing it again made me realize what I was doing incorrectly. I still can't get it to work with grid.draw(), but here is the solution with ggplot(): library(ggplot2) library(grid) library(maps) library(mapproj) library(png) library(RgoogleMaps)

Re: [R] Problem with readHTMLTable

2012-05-26 Thread R. Michael Weylandt
Odd... it doesn't throw an error for me. Two things possibly come to mind, but I don't know enough of the XML internals to know if they're reasonable. First is that something seems funny about your version number -- mine is simply XML_3.9-4 (no added -1); do you have any idea why yours would be

Re: [R] Assessing interaction effects in GLMMs

2012-05-26 Thread Ben Bolker
Luke Duncan luke.mangaliso.duncan at gmail.com writes: Dear R gurus I am running a GLMM that looks at whether chimpanzees spend time in shade more than sun (response variable 'y': used cbind() on counts in the sun and shade) based on the time of day (Time) and the availability of shade

[R] How do I modify sytem shell environment from within R?

2012-05-26 Thread Michael
Hi all, Could you please tell me how do I modify the environment variable (as below), from within R session? LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/mypath export LD_LIBRARY_PATH Thanks a lot! [[alternative HTML version deleted]] __

[R] path analysis and diagram (structural model)

2012-05-26 Thread Kristi Glover
Hi R- USER, I have been struggling to compute path analysis in R. I installed package sem and I tried to run the analysis but could not get a diagram. I have very big data set, here I just copied a sample of my data and I wanted to see how I can do the path analysis and get the patch diagram

Re: [R] How do I modify sytem shell environment from within R?

2012-05-26 Thread Jeff Newmiller
You don't, at least not for building libraries. Do it at the command line before you invoke the R commands. If you are using Cygwin, then be sure to use the Bash shell export command to set the variables. --- Jeff Newmiller

[R] barplot problem

2012-05-26 Thread jack hietpas
Hello- I am having trouble making a barplot... The y-values are: data=c(-0.0002129061,0.00,-0.0002699561,0.0163883061,0.04). The x-values are distance=c(0, 71, 172, 206, 292). The desired x-range is 0 and ~300+ (kilometers), range=seq(0,300,1) I would like to make a bar plot

Re: [R] barplot problem

2012-05-26 Thread R. Michael Weylandt
It sounds like you want something of the form plot(distance, data, type = h) But the lines might be too thin: try adding lwd = NUM to increase the thickness. Hope this helps, Michael On Sat, May 26, 2012 at 11:26 PM, jack hietpas mikros...@gmail.com wrote: Hello- I am having trouble making a

[R] a simple mixed model

2012-05-26 Thread array chip
Hi, I was reviewing a manuscript where a linear mixed model was used. The data is simple: a response variable y was measured for each subject over 3 time points (visit 1, 2 and 3) that were about a week apart between 2 visits. The study is a non-drug study and one of the objectives was to

[R] How to measure level of similarity of two data frames

2012-05-26 Thread Lamke
Hi group, I've been thinking of calculating euclidean distance between each column of a data frames that each consists of standardized numerical columns. However, I don't know if there's a way of summarizing the overall distance by some kind of metrics. If anyone know a proper way of doing so

Re: [R] subset columns from list with variable substitution

2012-05-26 Thread Jon Ween, MD
Thanks Don Please forgive my poor mail-liost etiquette. I had a couple of errors: 1) the counter logic in the loop was i in name list, my typo in the post. 2) There were typos in the variable list that caused the loop to crash when they were encountered. Thanks for your help, really

Re: [R] How to assign height value on bar plot?

2012-05-26 Thread Manish Gupta
My Question is how to write height length on each bar. -- View this message in context: http://r.789695.n4.nabble.com/How-to-assign-height-value-on-bar-plot-tp4631457p4631468.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Assessing interaction effects in GLMMs

2012-05-26 Thread arun
HI Luke, It would be better to ask this question on R mixed models (r-sig-mixed-mod...@r-project.org).  Just for curiosity (as I am doing a similar kind of light/dark response in insects), I am interested in the response variable (y) using cbind().  If I understand it correctly, you are using

Re: [R] how to check given number seq. is time series or not?

2012-05-26 Thread sagarnikam123
Yes,sir index of above/below numbers is time (both dataset are same) x 1 0.8890464 2 1.2272616 3 1.2272616 4 1.3578511 5 1.3578511 6 1.1070461 7 1.4424190 8 1.2277843 9 1.3578511 10 0.9708839 11 0.8221709 12 1.3578511 13 0.3588158 14 0.7742342 15 0.8221709 16

Re: [R] Problem with readHTMLTable

2012-05-26 Thread arun
Hi Sven, The function is working perfectly fine for me.  Could it be due to package masking the function?  Sometimes, these problems will get solved by restarting R console.  It happened to me a couple of times. A.K. - Original Message - From: Sven D sd...@hotmail.com To:

Re: [R] How to measure level of similarity of two data frames

2012-05-26 Thread Lamke
Thanks Ozgur. I will definitely try them out. Cheers, Kel -- View this message in context: http://r.789695.n4.nabble.com/How-to-measure-level-of-similarity-of-two-data-frames-tp4631466p4631478.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] how to check given number seq. is time series or not?

2012-05-26 Thread sagarnikam123
is this is time series ? x 1 2 2 2 3 2 4 2 5 2 6 2 7 2 8 2 9 2 10 2 11 2 12 2 -- View this message in context: http://r.789695.n4.nabble.com/how-to-check-given-number-seq-is-time-series-or-not-tp4631434p4631483.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] how to check given number seq. is time series or not?

2012-05-26 Thread Jeff Newmiller
could be. how would you know? --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...

Re: [R] How to assign height value on bar plot?

2012-05-26 Thread Jim Lemon
On 05/27/2012 04:30 AM, Manish Gupta wrote: My Question is how to write height length on each bar. Hi Manish, If don't have stacked bars, it is fairly simple: barpos-barplot(height,...) text(barpos,height,labels=height) usually with a bit of fooling around with the cex argument. If you want