[R] Two lines, two scales, one graph

o plot two variables against the same abscissa values. They have different scales. I've found how to make a second axis on the right for labeling, but not how to plot two lines at different scales. Thanks, - Rick r...@ece.pdx.edu __ R-help@r

Re: [R] Is there t.test with null hypothesis?

You need to include the argument "mu=1" (without parentheses). For example: > t.test(group1,group2, mu=1) for a two-sample independent groups t-test. If you type: > ?t.test you can see the help information for the t.test function. RIck On 09/08/2016 08:06 AM, Matti Vilj

Re: [R] Lavaan

Have you considered using the semPlot package? It works nicely with lavaan models (among other sem packages). There is also the DiagrammeR package. Rick On 06/23/2015 10:48 AM, DzR wrote: Dear Senior users of R/R Studio, I am very new to this environment hence am unable to plot the SEM

Re: [R] R vs. RStudio?

t within the first three hours of R development. Hope this helps, Rick Reeves Foxgrove Solutions, Ltd On 1/10/2015 6:47 PM, John Sorkin wrote: I urge you to try it. John John David Sorkin M.D., Ph.D. Professor of Medicine Chief, Biostatistics and Informatics University of Maryland School of M

Re: [R] HISTOGRAM

First, a histogram would not be appropriate (your data appear to be categorical - a histogram is for continuous numeric vales) - you would need a bar plot. You should make two vectors (one for the category names and the other for the frequencies) and use the barplot function. On Fri, Nov 9, 2018 a

Re: [R] chi-square test

, minimal, self-contained, reproducible code. You are using a Chi-squared test on a 124x2 matrix of values (not all integers) and many are zeros. The expected frequencies for many cells are very small (near zero, less than 1) hence the warning message. More importantly, does this application of the C

Re: [R] Logical operator in R

are different. If you want to have comparisons like this result in being considered equal, I think there is a way to use a fuzzy comparison but I don't remember the details. Rick -- Richard A. Bilonick, PhD Assistant Professor Dept. of Ophthalmology, School of Medicine Dept. of Biostatistic

[R] [R-pkgs] New Package: lrequire v0.1.3 - use modules to encapsulate and cache your R scripts

e work' hello <- function(person.name) { return (paste0('Hello, ', person.name, '!')) } module.exports = hello ____ File: main.R hello <- lrequire(welcome) hello('Rick') Note that the var

[R] nlsList {nlme} - control arguments problem

HS, data) {   xy     <- sortedXyData(mCall[["x"]],LHS,data)   min.s  <- min(y)   dif.s  <- max(y)-min(y)   dplt.s <- 0.5   p.s    <- -.20   value  <- c(min.s, dplt.s, dif.s, p.s)   names(value) <- mCall[c("min","dplt","dif","p&

[R] Tck/tk help

I am trying to create a simple window that displays a dataframe. I was able to put together a function from a couple of examples on the web. For the most part my function works. The only problem is with resizing the window. I found a post where someone suggested using with(env,tkpack("configur

[R] -network- package in R: removing nodes

Dear all, If you have experience with the -network- package, can you show me how to remove nodes based on their degrees? For instance, removing all nodes with an out-degree or in-degree of 1? Thank you very much for your help! Rick __ R-help@r

[R] -network- package in R: removing nodes

Dear all, If you have experience with the -network- package, can you show me how to remove nodes based on their degrees? For instance, removing all nodes with an out-degree or in-degree of 1? Thank you very much for your help! Rick __ R-help@r

[R] Installing Rmpi on Centos 5.2

Rmpi: R CMD INSTALL /home/rick/Rmpi_0.5-7.tar.gz --with-mpi=/usr/lib64/lam/include But I get the following error message: /usr/bin/ld: skipping incompatible /usr/lib/lam/lib/libmpi.so when searching for -lmpi I'm not sure what else to install/unins

Re: [R] CFA with lavaan or with SEM

Not sure if you are aware of the OpenMx SEM package (http://openmx.psyc.virginia.edu/). It's a very full-featured structural equation modeling package. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

[R] SEM and count data

ndle SEM models. I'm thinking I may have to do a two-step process and run SEM package to get factor scores, and then put in as a predictor in count outcomes package. Thanks! Rick [[alternative HTML version deleted]] __ R-help@r-proj

[R] Problem with coordinates when trying to draw lines into a raster (image) file

used these extended graphics capabilities much at all, so I really am getting frustrated Regards and thanks in advance, Rick -- The code segment in que

[R] R - Time Series Regression with a p-value check for each additional added point.

ate. Thanks for your help, Rick [[alternative HTML version deleted]] __ 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 commen

[R] plotCI overlay

I'm using the plotCI function and I'd like to overlay additional means with CIs onto an existing plotCI-created plot in a different color. Is this possible? Thanks. Rick __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

[R] Generate multivariate normal data with a random correlation matrix

Hi All. I'd like to generate a sample of n observations from a k dimensional multivariate normal distribution with a random correlation matrix. My solution: The lower (or upper) triangle of the correlation matrix has n.tri=(d/2)(d+1)-d entries. Take a uniform sample of n.tri possible correlations

[R] Generate multivariate normal data with a random correlation matrix

Hi All. I'd like to generate a sample of n observations from a k dimensional multivariate normal distribution with a random correlation matrix. My solution: 1) The lower (or upper) triangle of the correlation matrix has n.tri=(d/2)(d+1)-d entries. 2) Take a uniform sample of n.tri possible correl

Re: [R] Generate multivariate normal data with a random correlation matrix

correlation matrices of size k equally likely using your approach? It would seem so based on your description. If so, it is a way cool solution. Rick On Thu, Feb 10, 2011 at 12:18 PM, wrote: > If you want a random correlation matrix, why not just generate random data > and acce

[R] Missing times chron

70 (the default date by chron) and making a subset and only replacing values that are not equal to zero doesn't work either: d <- x ## above coding for chron with x replaced with d ## mat <- match(x$ID, d$ID) x$dateofdeath <- d$dateofdeath[mat]

Re: [R] FIML with missing data in sem package

On 12/01/2011 07:18 AM, John Fox wrote: > To:r-help@r-project.org > Subject: [R] FIML with missing data in sem package > You should check out the OpenMx R package. Just search for OpenMx and SEM. You can download from the web site. It does FIML and is an excellent SEM package.

[R] args() function does not list function prototype for locally-produced/installed R package

quot;function" not found. Question: If it is possible, what changes do I make to my package and install process to get args() to display the function prototype(s) for the functions in my local package? Thanks, Rick Reeves __ R-help@r-proje

Re: [R] args() function does not list function prototype for locally-produced/installed R package

Duncan: Indeed, adding the function names to the NAMESPACE file (in the form: export(functionName)) solved the problem, and the args() function displays the function argument lists. Thanks, Rick Reeves On 12/5/2011 2:49 PM, Duncan Murdoch wrote: On 11-12-05 3:04 PM, Rick Reeves wrote

[R] PMML for Cox Regression

this cox regression model. Can someone point us to the proper documentation for scoring this model? Thanks in advance, Rick Pechter http://www.dmg.org/PMML-3_2"; xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"; xsi:schemaLocation="http://www.dmg.org/PMML-3_2 http://

Re: [R] How to use IML with R and SAS

The SAS/IML documentation on calling R is at http://support.sas.com/documentation/cdl/en/imlug/64248/HTML/default/viewer.htm#r_toc.htm Other ways to get help include calling SAS Technical Support (it's free) and posting your question to the SAS/IML Discussion forum at http://communities.sas.com/com

[R] Getting R and x11 to work

nd do computations. When I try to do any graphics using x11, I get the message: unable to load shared library '/usr/lib/R/modules//R_X11.so': /usr/lib/R/modules//R_X11.so: undefined symbol: cairo_image_surface_get_data Does anyone have an idea on how to fi

[R] Installing R-2.6.2 on Windows multiuser

l/reinstall and had the same problem. Thanks for any suggestions. -- Rick Voland [EMAIL PROTECTED] __ 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

Re: [R] Installing R-2.6.2 on Windows multiuser

Files\R\R-2.6.2 (used for previous versions of R) and still had the same permissions issue. I also tried setting R_USER environment variable as I previously indicated. -- Rick Voland [EMAIL PROTECTED] __ R-help@r-project.org mailing list https

[R] NLS plinear question

582 residual sum-of-squares: 555915 Number of iterations to convergence: 11 Any idea why having a zero for the first value of X causes this problem? Thanks in advance, Rick DeShon [[alternative HTML version deleted]] __ R-help@r-projec

[R] Efficient computation of average covariance matrix over a list

cient way to compute the average covariance matrix over the list members in "lcov"? Thanks in advance, Rick [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLE

[R] Efficient computation of average covariance matrix over a list

cient way to compute the average covariance matrix over the list members in "lcov"? Thanks in advance, Rick DeShon __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-

[R] Ordinal Package Errors

e: 'rmutil' The following object(s) are masked from package:base : units Error in library.dynam("ordinal", pkg, lib) : shared library 'ordinal' not found Error in library(ordinal) : .First.lib failed for 'ordinal' Does anyone know ho

Re: [R] Ordinal Package Errors

I corrected it as described above. It would have been a whole lot easier if the package was in the contributed package repository. Rick B. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] Power analysis for MANOVA?

retrospective" power analysis leads to logical contradictions with respect to the confidence intervals and p-values from the analysis of the data. In other words, DON'T DO IT! All the information is contained in the confidence intervals which are based on the observed data - an after the f

[R] replace zeros in a block diagonal matrix with small random values?

Hi All. Imagine you have a large block diagonal matrix. I'd like to replace the zeros in this matrix with small random (runif) numbers. Any ideas for a simple and efficient way to do this? Best regards, Rick DeShon __ R-help@r-project.org ma

[R] nlsList (nlme) error

;"8"<"30"<..: 41 11 33 22 4 27 5 2 37 19 ... $ trial: int 1 1 1 1 1 1 1 1 1 1 ... $ ACC : int 1 0 1 1 1 0 1 1 1 1 ... $ RT : int NA 1309 544 654 NA 441 882 1097 898 ... $ block: int 1 1 1 1 1 1 1 1 1 1 ... - attr(*, "formula")=Class 'formul

[R] Errors running gam examples

I was able to get step.gam to work on its example but it doesn't work on my model. Has anyone else run into this? Rick B. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide h

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

incorporate a Spatial(Points or Polygons)DataFrame into the spplot display list? It seems as though many scientists would like to create such plots without resorting to GRASS or another GIS. Thanks for any advice, Rick Reeves Rick Reeves Scientific Programmer / Analyst National Center for

[R] This site provided an excellent answer to my raster / vector plot question

I can recommend it highly, esp for those of us in the geo spatial realm: http://casoilresource.lawr.ucdavis.edu/drupal/node/442 Thanks do Dylan for answering... RR -- Rick Reeves Scientific Programmer / Analyst National Center for Ecological Analysis and Synthesis UC Santa Barbara

[R] Cannot Install rimage

ges are in /tmp/RtmpGDx12u/downloaded_packages Warning message: installation of package 'rimage' had non-zero exit status in: install.packages("rimage") I cannot find "config.log" that is mentioned. Can someone plea

[R] Executing a Function in a Loop With a Changing Value for an Argument

) } where "vars" is a vector of names for items in a data.frame. If I use "as.symbol" to create names, the evaluation never works correctly. Is there a way to do this? Any suggestions would be appreciated. Rick B. __ R-help@r-proj

Re: [R] Executing a Function in a Loop With a Changing Value foran Argument

t will allow my function to work normally (just as it would if I wrote out the function calls with the names manually inserted)? (I would assume that any function that generates a character variable from the variable name is not going to produce the desired result.) Rick B. __

Re: [R] Executing a Function in a Loop With a ChangingValue foran Argument

n vars) with(dat, eval(subst(ex, x = as.symbol(i),y=2))) > myfunc <- function(y,x) { + nam <- deparse(substitute(x)) + val <- mean(x) + cat("mean(", nam, ") =", val,"y=",y, "\n") + invisible(val) + } > ex <- quote(myfunc(y,

[R] Confidence Intervals for Random Effect BLUP's

d. Is the approach shown below correct? Rick B. # Orthodont is from nlme (can't have both nlme and lme4 loaded at same time!) # OrthoFem<-Orthodont[Orthodont$Sex=="Female",] # http://tolstoy.newcastle.edu.au/R/help/06/03/23758.html library(lme4) fm1OrthF. <- lmer(distan

Re: [R] Confidence Intervals for Random Effect BLUP's

ess the information from lmer to compute intervals for the BLUPS that reflect the uncertainty in these estimates. The code I showed was taken from someone who had the same question and the various replies. Rick B. __ R-help@r-project.org mailing list ht

Re: [R] Confidence Intervals for Random Effect BLUP's

.hokudai.ac.jp/~kubo/Rdoc/library/lme4/html/lmer-class.html bVar: A list of the diagonal inner blocks (upper triangles only) of the positive-definite matrices on the diagonal of the inverse of ZtZ+Omega. With the appropriate scale factor (and conversion to a sy

[R] Using lme (nlme) to find the conditional variance of the random effects

appropriate scale factor (and conversion to a symmetric matrix) these are the conditional variance-covariance matrices of the random effects. Is there anything similar in the nlme package using the lme function? Rick B. __ R-help@r-project.org mailing list https

Re: [R] Using lme (nlme) to find the conditional variance of therandom effects

ilto:[EMAIL PROTECTED] On Behalf Of Rick Bilonick > > Sent: Monday, November 12, 2007 4:40 PM > > To: R Help > > Subject: [R] Using lme (nlme) to find the conditional > > variance of therandom effects > > > > Using lmer in the lme4 package, you can compute t

Re: [R] Using lme (nlme) to find the conditional variance of therandom effects

On Tue, 2007-11-13 at 01:03 -0500, Rick Bilonick wrote: > > Is there some way to get ranef with postVar=TRUE to show what the > variances are, or what the lower and upper bounds are? qqmath makes nice > plots but I need to obtain the numerical values. > > Rick B. > I fou

[R] How to remove index from list after split?

In the following example, how can I drop the group index from the list after I perform a split? n <- 3 nn <- 10 g <- factor(round(n * runif(n * nn))) x <- rnorm(n * nn) + sqrt(as.numeric(g)) df<- data.frame(g,x) df.s <- split(df,g) Thanks! Rick DeShon

[R] Rserve - Request assistance with installation, confirmation thereof and starting (R 3.5.1 on 64-bit Win7 Pro)

Hello, I am attempting to install two R packages with specific version numbers. These are Rserve_1.8-0.zip and MASS_7.3-45.zip. As the file extension suggests, I am installing packages from local zipfiles in the GUI Packages menu. R Console displays positive feedback when MASS is loaded in this