Re: [R] [FORGED] turning the output of cut into a waffle plot

2017-01-17 Thread Rolf Turner
On 18/01/17 12:21, Jim Lemon wrote: Hi all, A few days ago I offered a suggestion on how to display the initial values that were cut into a factor as a waffle plot. As Rolf Turner noted, a major problem for users would be constructing the matrix that was fed to the color2D.matplot function. Here

Re: [R] Adding regression line to each individual plot in a window with multiple plots

2017-01-17 Thread Jim Lemon
Hi Jake, My mistake, that should be: par(xpd=FALSE) abline(lm(...)) par(xpd=NA) I suspect that the odd looking regression lines are due to points off the plot as you have specified xlim and ylim for your plots. Jim On Wed, Jan 18, 2017 at 2:22 PM, Jake William Andrae

Re: [R] Adding regression line to each individual plot in a window with multiple plots

2017-01-17 Thread Jim Lemon
Hi Jake, In the second line of your script, you set xpd=NA. That means that abline will draw a line across the entire display region rather than restricting it to the plot region. If you must set this at the beginning, then add these lines: par(xpd=TRUE) abline(lm(...)) par(xpd=NA) for each

Re: [R] Adding regression line to each individual plot in a window with multiple plots

2017-01-17 Thread Jim Lemon
Hi Jake, As I don't have your data set, try this: attach(mtcars) plot(mpg~disp,xaxs="i",yaxs="i") abline(lm(mpg~disp)) Jim On Wed, Jan 18, 2017 at 12:04 PM, Jake William Andrae wrote: > Hi Everyone, > > > > I've constructed a script that adds multiple plots to the

[R] Adding regression line to each individual plot in a window with multiple plots

2017-01-17 Thread Jake William Andrae
Hi Everyone, I've constructed a script that adds multiple plots to the plot window, but I'm having a bit of trouble adding a regression line to each individual plot. Of course, the regression lines will vary depending on the variables plotted against one another. I've attached the script.

[R] the difference between R and other statistical software when using maxilla package

2017-01-17 Thread Alaa Sindi
Hi all, Is there any differences in estimation results when using maxlik package in R and other statistical software? Thanks Alaa __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help

[R] turning the output of cut into a waffle plot

2017-01-17 Thread Jim Lemon
Hi all, A few days ago I offered a suggestion on how to display the initial values that were cut into a factor as a waffle plot. As Rolf Turner noted, a major problem for users would be constructing the matrix that was fed to the color2D.matplot function. Here is a fairly general purpose function

Re: [R] Dimensionality reduction with ProDenICA

2017-01-17 Thread Rolf Turner
On 18/01/17 11:17, Neverstop . wrote: Hello, I have a dataset with many variables and I'd like to do dimensionality reduction with Independent Component Analysis. There are many statistical methods to estimate the latent variables of the ICA model. I'm trying the R package "ProDenICA" that

[R] Dimensionality reduction with proDenICA

2017-01-17 Thread Neverstop .
Hello, I have a dataset with many variables and I'd like to do dimensionality reduction with Independent Component Analysis. There are many statistical methods to estimate the latent variables of the ICA model. I'm trying the R package "proDenICA" that implements the penalized maximum

Re: [R] [FORGED] Re: output

2017-01-17 Thread Rolf Turner
On 18/01/17 10:28, Marc Schwartz wrote: Hi Val, Alternatively, there are other Excel file I/O packages available on CRAN if you would prefer to consider other options as well. They will typically require Java being installed, rather than Perl, as I and Greg do. The read_excel() function

Re: [R] Cronbach's Alpha

2017-01-17 Thread John Dougherty
On Tue, 17 Jan 2017 12:45:10 -0400 Jim Silverton wrote: > Dear all, > > I have 12 companies and I am developing a scale for innovation. I > want to check the reliability of my scale. I delivered it to > different (engineers, managers etc) people in each company and an >

Re: [R] output

2017-01-17 Thread Marc Schwartz
Hi Val, Presuming that the Excel file that you are trying to read has a second worksheet, which is what the read.xls() command you are using is trying to do, the problem may be that XLSX file support has not been installed for the gdata package, which is what I presume you are using. You do

Re: [R] output

2017-01-17 Thread Val
Hi Marc and all, Last time you suggest me to use WriteXLS function to write more than 65,000 row in excel. Creating the file worked fine. Now I wanted to read it using the WriteXLS function but have a problem,. The file has more than one sheets. Here is the script and the error message.

[R] Meta-regression

2017-01-17 Thread Carlos R. Moreira Maia
Dear all, I'm running a meta-regression but, as a newbie, I'm facing few problems to interpret the outputs from metareg. I wonder if you can inform some literature where I can find how to understand and interpret the meta-regression analysis from R. These are the books/papers I have, but they

Re: [R] Qvalue package: I am getting back 1, 000 q values when I only want 1 q value.

2017-01-17 Thread Jay Tanzman
What you're doing makes no sense. Given p-values p_i, i=1...n, resulting from hypothesis tests t_i, i=1...n, the q-value of p_i is the expected proportion of false positives among all n tests if the significance level of each test is α=p_i. Thus a q-value is only defined for an observed p-value.

[R] Current R terminal width

2017-01-17 Thread Jared Studyvin
Hello, I have found when using R in a unix environment that Sys.getenv("COLUMNS") will return the current R terminal width. This does not work on a Window OS. How can this be done? I'm looking for the Windows OS equivalent of options(width=Sys.getenv(" COLUMNS")) Thanks, *Jared Studyvin, PhD *

[R] Cronbach's Alpha

2017-01-17 Thread Jim Silverton
Dear all, I have 12 companies and I am developing a scale for innovation. I want to check the reliability of my scale. I delivered it to different (engineers, managers etc) people in each company and an unequal number of values. For example: Company 12 engineers 2 managers Company 24

Re: [R] Receiving NaN message

2017-01-17 Thread David L Carlson
And (re: "dat[,ID]"), the page on "Extract or Replace Parts of an Object": ? Extract - David L Carlson Department of Anthropology Texas A University College Station, TX 77840-4352 -Original Message- From: R-help [mailto:r-help-boun...@r-project.org]

Re: [R] Receiving NaN message

2017-01-17 Thread Rui Barradas
Hello, And in addition to what Jim and Petr said, take a look at ?%in%. Rui Barradas Em 17-01-2017 10:50, PIKAL Petr escreveu: Hi And above what Jim adviced you should look also at ?match help page. Cheers Petr -Original Message- From: R-help [mailto:r-help-boun...@r-project.org]

Re: [R] Receiving NaN message

2017-01-17 Thread PIKAL Petr
Hi And above what Jim adviced you should look also at ?match help page. Cheers Petr > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Robert > Piliero > Sent: Monday, January 16, 2017 9:10 PM > To: r-help@r-project.org > Subject: [R] Receiving NaN