[R] SKAT for survey data

2015-02-19 Thread Jeannette Simino
Hello all! Has anyone modified SKAT or another rare variant analysis packages to incorporate a subject's sampling weight? I am trying to analyze genetic data from a study that chose a) all people from a previous ancillary study; b) those with an atypical diagnosis; and c) a small number of

[R] looping multipanel plots to different figures

2015-02-19 Thread efisio solazzo
Dear, cannot find a way to direct multipanel plots to different figures (files) while within a loop. Say, the loop creates two plots each step: one plot should go to figure 1 and the other to figure 2. Same for the next steps of the loop: the plots should go to figure 1 and figure 2 in a

[R] Removing objects except user-defined functions

2015-02-19 Thread philippe massicotte
Dear R users. I would like to remove all object from my workspace except the function I have defined. However, is I use rm(list = ls()) everything is cleared. I was thinking to typeof to get information about objects, but I could not get it working right. Thank in advance, Phil

[R] find a previous command that starts with a certain letter

2015-02-19 Thread CasasPrat,Merce [Ontario]
Hi, I'm new to R (was using Matlab previously) Matlab has a very convenient utility that I haven't been able to find in R. When working interactively, you can find a previous command that starts with a certain letter by typing such a letter and then pressing the cursor up key. In R I've seen

Re: [R] API request from R

2015-02-19 Thread Robert Baer
On 2/19/2015 8:06 AM, Barry Rowlingson wrote: On Wed, Feb 18, 2015 at 11:44 AM, Mittal Ashra via R-help r-help@r-project.org wrote: Dear All, Apologies for mailing it to the whole crowd. This is Mittal, presently working in a Project where we have build a platform for displaying

Re: [R] API request from R

2015-02-19 Thread Barry Rowlingson
On 19 Feb 2015 20:11, Robert Baer rb...@atsu.edu wrote: On 2/19/2015 8:06 AM, Barry Rowlingson wrote: On Wed, Feb 18, 2015 at 11:44 AM, Mittal Ashra via R-help r-help@r-project.org wrote: Dear All, Apologies for mailing it to the whole crowd. This is Mittal, presently working in a Project

Re: [R] About Read in .csv Files with Chinese Characters

2015-02-19 Thread JS Huang
Hi Nicholas, I am not sure how much the following link can help but at least it is related to your question. http://r.789695.n4.nabble.com/Reading-Chinese-Language-GB2312-Input-td4647581.html Good luck! JS -- View this message in context:

Re: [R] Procrustes

2015-02-19 Thread Tara Dirilgen
Hi again, Just to clarify my question from previous email... Example script: A - data[, c(1,2,3,4)] B - data[, c(5,6,7,8)] library(vegan) vare.proc - procrustes(A,B, scale=FALSE, symmetric=FALSE) vare.proc summary(vare.proc) plot(vare.proc) #plot(vare.proc, kind=2) residuals(vare.proc)

Re: [R] Raster Help

2015-02-19 Thread Sven E. Templer
Without (example) code it is hard to follow... use ?dput to present some data (subset). But if it is data.frames you are dealing with (for sure with read.csv, but not so sure at all with raster maps), give this a try: ?merge On 19 February 2015 at 17:44, Simon Tarr simon.t...@adtrak.co.uk wrote:

[R] how to convert raster images to .kasc format

2015-02-19 Thread Angela Smith
Hi R user, I was trying to convert raster images into *kasc format using adehabitat package in R. but I could not convert it. I spent a lot of time but no luck. would you mind to give some hints to convert this example data?. I would really appreciate for your help. cheers, AS #

[R] Subsetting a list of lists using lapply

2015-02-19 Thread Aron Lindberg
Hi Everyone, I'm working on a thorny subsetting problem involving list of lists. I've put a dput of the data here: https://gist.githubusercontent.com/aronlindberg/b916dee897d051ac5be5/raw/a78cbf873a7e865c3173f943ff6309ea688c653b/dput I can get one intense of the element I want this

Re: [R] Averaging column scores when participants vary in number of observations

2015-02-19 Thread Chad Danyluck
I have a data set that includes the identity of a number of Video Coders who scored participants' behaviors in a video. Every participant was scored once, but some participants were randomly assigned to have their data scored twice so I could calculate inter-rater reliabilities. I have completed

[R] Issue writing correct code for Coursera assignment--performing matrix inversion

2015-02-19 Thread Samuel Colon
Hello All, I'm new to this mailing list, so please let me know if I've committed any posting faux-pas. I'm working on an assignment for my Coursera course; please see my code below in which I have tried to write two functions--to perform the task of matrix inversion and then caching that data.

Re: [R] Removing objects except user-defined functions

2015-02-19 Thread Stéphane Adamowicz
There is a function keep() in package gdata for this purpose Le 19 févr. 2015 à 10:25, philippe massicotte pmassico...@hotmail.com a écrit : Dear R users. I would like to remove all object from my workspace except the function I have defined. However, is I use rm(list = ls()) everything

Re: [R] find a previous command that starts with a certain letter

2015-02-19 Thread Duncan Murdoch
On 18/02/2015 5:06 PM, CasasPrat,Merce [Ontario] wrote: Hi, I'm new to R (was using Matlab previously) Matlab has a very convenient utility that I haven't been able to find in R. When working interactively, you can find a previous command that starts with a certain letter by typing such

Re: [R] looping multipanel plots to different figures

2015-02-19 Thread Jim Lemon
Hi efisio, I read this as wanting to start a new graphics device, then set some plot parameters, display two plots and then close the graphics device at each iteration of the loop. If so, plot_filenames-c(plot1.png,plot2.png,plot3.png) for(plotfn in plot_filenames) { png(plotfn)

Re: [R] looping multipanel plots to different figures

2015-02-19 Thread efisio solazzo
Thanks Jim, actually I need to keep open two devices at the same time, and within the loop access either of them in alternation. In MatLab there is the command Figure(#) which keeps track of the open devices and direct the output of the plot to whichever of them. For example:

Re: [R] looping multipanel plots to different figures

2015-02-19 Thread Jim Lemon
Hi efisio, Okay, you can switch devices using the dev.* functions in the grDevices package. If you only have two devices open at one time, this is not too difficult: #open both devices png(...) par(mfrow=c(1,2)) png(...) par(mfrow=c(1,2)) hist(...) dev.set(dev.next()) hist(...)

Re: [R] Metro_Hastings

2015-02-19 Thread hms Dreams
Any suggestions :( ??From: cute_loo...@hotmail.com To: r-help@r-project.org Subject: Metro_Hastings I wrote my code again Date: Sun, 15 Feb 2015 21:47:25 +0300 Hi again :) I wrote my code here: library(MHadaptive) baysianlog=function (param,data) { alpha=param[1]

[R] How to analyse nonlinear response to categorical and quantitative explanatory variables?

2015-02-19 Thread Jan-Ulrich Kreft
Dear list I have data from a collaborator who has used DesignExpert to design the experiment and analyse the data but no longer has access to this software and does not know exactly what the software did and why. So I’m now trying to analyse the data in R but can't quite decide what to do.

Re: [R] API request from R

2015-02-19 Thread Barry Rowlingson
On Wed, Feb 18, 2015 at 11:44 AM, Mittal Ashra via R-help r-help@r-project.org wrote: Dear All, Apologies for mailing it to the whole crowd. This is Mittal, presently working in a Project where we have build a platform for displaying recommendations and the results are based on the

Re: [R] About Read in .csv Files with Chinese Characters

2015-02-19 Thread JS Huang
Hi, I tried your file in my Windows7 pc. It worked fine except the display of the variable name. R version 3.1.2 (2014-10-31) -- Pumpkin Helmet Copyright (C) 2014 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) Hello - read.csv(Hello.csv);Hello

[R] Noob: How to upgrade from 3.0.2 on Ubuntu 14.04 LTS?

2015-02-19 Thread brian piercy
One of my favorite blogs (AnalyzeCore) uses the dplyr library which fails in my environment (R 3.0.2.) I've tried upgrading my R package to get around this issue, to no avail. I know I'm making a simple mistake but haven't solved it. A simple $ sudo apt-get install r-base --upgrade didn't get the

Re: [R] Noob: How to upgrade from 3.0.2 on Ubuntu 14.04 LTS?

2015-02-19 Thread John Kane
This worked for me the last time I tried it and I'm running Ubuntu 14.04. (Note to self-time to upgrade to 14.10) Installing or updating R sudo add-apt-repository ppa:marutter/rrutter sudo apt-get update sudo apt-get install r-base r-base-dev from

[R] 'java' had status 1

2015-02-19 Thread Erica Cseko Nolasco
Dear all, I´m trying to perform a modeling with Maxent, but I keep getting the Warning message: running command 'java' had status 1. I reinstalled Java jdk, I raised the memory to 1024, but it still gives me the same Warning message. I´m testing with the data from the

Re: [R] Change error bar length in barplot2

2015-02-19 Thread John Kane
https://github.com/hadley/devtools/wiki/Reproducibility http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example John Kane Kingston ON Canada -Original Message- From: jmadi...@yahoo.fr Sent: Tue, 17 Feb 2015 16:46:34 + (UTC) To:

Re: [R] Procrustes

2015-02-19 Thread Pierrick Bruneau
Hi Tara, Providing a simple example script that reproduces your case and using it to support your question would increase your chances to obtain an answer. Best, Pierrick On Thu, Feb 19, 2015 at 5:43 PM, Tara Dirilgen tara.diril...@ucdconnect.ie wrote: I have been using R to calculate the

[R] Raster Help

2015-02-19 Thread Simon Tarr
Hello everyone, I need a little help with some R syntax to complete what (I think) is a fairly straightforward task- hopefully someone can assist! I have a raster map of the UK which is split into postcode areas (e.g. DE, NG, NR etc. 127 postcodes in total). I have installed the package

Re: [R] ggplot2 shifting bars to only overlap in groups

2015-02-19 Thread John Kane
Ah yes, with the full data set that looks rather good, indeed. I tried it with the simulated data but there was not enough variation in the data to give any useful indication of what was happening. RE the quartile information, I don't know, it might be of use but it also might add to much

[R] Procrustes

2015-02-19 Thread Tara Dirilgen
I have been using R to calculate the significance of Procrustes correlations. With one series of data, where there are five cases, the value returned for the correlation coefficient is one although there are differences as shown by the procrustes error graph. Is there a statistical reason for

Re: [R] Averaging column scores when participants vary in number of observations

2015-02-19 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Chad Danyluck Sent: Thursday, February 19, 2015 12:33 PM To: r-help@r-project.org Subject: Re: [R] Averaging column scores when participants vary in number of observations I have a data set that

Re: [R] Subsetting a list of lists using lapply

2015-02-19 Thread Rolf Turner
On 20/02/15 08:45, Aron Lindberg wrote: Hi Everyone, I'm working on a thorny subsetting problem involving list of lists. If you think this is thorny you ain't seen nothin' yet! But note that you've got a list of lists of lists ... i.e. the nesting is at least 3 deep. I've put a dput of

Re: [R] Issue writing correct code for Coursera assignment--performing matrix inversion

2015-02-19 Thread Duncan Murdoch
On 19/02/2015 3:13 PM, Samuel Colon wrote: Hello All, I'm new to this mailing list, so please let me know if I've committed any posting faux-pas. Yes, you should probably be using the Coursera support resources instead of this list. Duncan Murdoch I'm working on an assignment for my

Re: [R] API request from R

2015-02-19 Thread Mittal Ashra via R-help
Dear All Thanks for the reply. RegardsMittal On Friday, 20 February 2015 1:56 AM, Barry Rowlingson b.rowling...@lancaster.ac.uk wrote: On 19 Feb 2015 20:11, Robert Baer rb...@atsu.edu wrote: On 2/19/2015 8:06 AM, Barry Rowlingson wrote: On Wed, Feb 18, 2015 at 11:44 AM, Mittal

Re: [R] Averaging column scores when participants vary in number of observations

2015-02-19 Thread JS Huang
Hi, Another implication: data1 Observation Participant.ID Video.Coder Score 1 A 1 Donald 4 2 B 1 Tracy 5 3 C 2 Donald 6 4 D 3 Sam 2 5 E

[R] Custom function causing an error with returning a ggplot object and Error in eval(expr, envir, enclos) : object '...' not found

2015-02-19 Thread jcrosbie
I'm trying to create a custom function to return a chart object. This function seems to be having an error with calculating min/max/etc in the ggplot object. If I run the code for the ggplot not inside a custom function it works. To reproduce this error after I need to clear the memory with

Re: [R] Subsetting a list of lists using lapply

2015-02-19 Thread Charles Berry
Aron Lindberg aron.lindberg at case.edu writes: Hi Everyone, I'm working on a thorny subsetting problem involving list of lists. I've put a dput of the data here: https://gist.githubusercontent.com/aronlindberg/b916dee897d051ac5be5/

[R] Fwd: How to analyse nonlinear response to categorical and quantitative explanatory variables?

2015-02-19 Thread Bert Gunter
Statistical methodology questions are generally off topic here. This list is about R programming. I suggest you post to a statistics site like stats.stackexchange.com instead. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 Data is not information. Information is not

[R] Density plots for Chi-square and F-distribution on my data

2015-02-19 Thread varin sacha
Dear R-helpers, I want to compute the density plot (probability plot) of the Chi-square distribution. My 2 categorical variables are gender (male, female) and colors of the eyes (blue, green and brown). The sample size n = 100. The proportions are the following : male and blue eyes : 10% male

Re: [R] Density plots for Chi-square and F-distribution on my data

2015-02-19 Thread varin sacha
Arg, just been caught ! No, just joking, it is a long time I don't have homework anymore. I keep on giving myself homeworks life long yes this is true. Anyways thanks for your response. Best, - Mail original - De : Bert Gunter gunter.ber...@gene.com À : varin sacha

Re: [R] find a previous command that starts with a certain letter

2015-02-19 Thread Loris Bennett
CasasPrat,Merce [Ontario] merce.casasp...@ec.gc.ca writes: Hi, I'm new to R (was using Matlab previously) Matlab has a very convenient utility that I haven't been able to find in R. When working interactively, you can find a previous command that starts with a certain letter by typing such

Re: [R] Density plots for Chi-square and F-distribution on my data

2015-02-19 Thread Bert Gunter
I think you might do well to consult a local statistical expert, as you appear to be out of your statistical depth here. Also, is this homework? If so, this is not a homework site (although posters sometimes get help anyway). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650)

[R] package 'corrplot'

2015-02-19 Thread Antonello Preti
I'm using the package 'corrplot' for a figure to be used in an academic article. The figure must be black-and-white. I was able to produce a figure with some shades of grey (less than fifty, admitedly). However, the legend is still in red color. I need the legend be in black color. How can I get

Re: [R] Removing objects except user-defined functions

2015-02-19 Thread Henrik Bengtsson
On Thu, Feb 19, 2015 at 1:25 AM, philippe massicotte pmassico...@hotmail.com wrote: Dear R users. I would like to remove all object from my workspace except the function I have defined. However, is I use rm(list = ls()) everything is cleared. I was thinking to typeof to get information

[R] non-terminal token lacking children from utils::getParseData

2015-02-19 Thread B Tyner
Hi, I have run across a source file for which the return value of getParseData() includes a record having FALSE for $terminal, yet it is not the parent of any other tokens. Before I spend time constructing a reproducible example, I wanted to verify that this is in fact unexpected behavior (under

Re: [R] package 'corrplot'

2015-02-19 Thread Duncan Murdoch
On 19/02/2015 6:13 PM, Antonello Preti wrote: I'm using the package 'corrplot' for a figure to be used in an academic article. The figure must be black-and-white. I was able to produce a figure with some shades of grey (less than fifty, admitedly). However, the legend is still in red color.

Re: [R] non-terminal token lacking children from utils::getParseData

2015-02-19 Thread Benjamin Tyner
I tried to reduce the offending portion as best I could to a more-or-less minimal example (1136 bytes), which can be downloaded via: wget https://www.dropbox.com/s/74rgxr5x2aalr99/badstring.R then once in R, b - parse(file = ~/badstring.R, keep.source = TRUE) d - getParseData(b,

Re: [R] package 'corrplot'

2015-02-19 Thread David Winsemius
On Feb 19, 2015, at 3:13 PM, Antonello Preti wrote: I'm using the package 'corrplot' for a figure to be used in an academic article. The figure must be black-and-white. I was able to produce a figure with some shades of grey (less than fifty, admitedly). However, the legend is still in red

Re: [R] non-terminal token lacking children from utils::getParseData

2015-02-19 Thread Duncan Murdoch
On 19/02/2015 6:31 PM, B Tyner wrote: Hi, I have run across a source file for which the return value of getParseData() includes a record having FALSE for $terminal, yet it is not the parent of any other tokens. Before I spend time constructing a reproducible example, I wanted to verify that