Re: [R] HISTOGRAM

2018-11-09 Thread Rick Bilonick
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

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

2016-09-08 Thread Rick Bilonick
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 Viljamaa wrote: I’m

Re: [R] Logical operator in R

2016-01-22 Thread Rick Bilonick
On 01/22/2016 10:46 AM, li li wrote: Hi all, I encountered the following strange phenomenon. For some reason, the obs_p[1] and res1$st_p[89] have the same value but when I run "==", it returns FALSE. Can anyone help give some explanation on this? Thanks very much! Hanna obs_p[1]

Re: [R] Lavaan

2015-06-24 Thread Rick Bilonick
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] chi-square test

2014-09-15 Thread Rick Bilonick
On 09/15/2014 10:57 AM, eliza botto wrote: Dear useRs of R, I have two datasets (TT and SS) and i wanted to to see if my data is uniformly distributed or not?I tested it through chi-square test and results are given at the end of it.Now apparently P-value has a significant importance but I

Re: [R] CFA with lavaan or with SEM

2013-02-01 Thread Rick Bilonick
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

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

2011-12-01 Thread Rick Bilonick
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. Rick B .

[R] Installing Rmpi on Centos 5.2

2009-03-26 Thread Rick Bilonick
I am trying to install the R package Rmpi which needs libmpi. I've installed openmpi and lam in Centos 5.2: [r...@rab45-1 /]# rpm -qv openmpi openmpi-1.2.5-5.el5 openmpi-1.2.5-5.el5 [r...@rab45-1 /]# rpm -qv lam lam-7.1.2-14.el5 lam-7.1.2-14.el5 I'm using this to install Rmpi: R CMD INSTALL

Re: [R] Power analysis for MANOVA?

2009-02-01 Thread Rick Bilonick
On Wed, 2009-01-28 at 21:21 +0100, Stephan Kolassa wrote: Hi Adam, first: I really don't know much about MANOVA, so I sadly can't help you without learning about it an Pillai's V... which I would be glad to do, but I really don't have the time right now. Sorry! Second: you seem to be

[R] Ordinal Package Errors

2009-01-14 Thread Rick Bilonick
I'm trying to install the ordinal package (http://popgen.unimaas.nl/~plindsey/rlibs.html). I downloaded ordinal03.tgz and untarred it. rmutil was previously installed (and appears to work ok.) Then I installed ordinal: [r...@localhost ~]# R CMD INSTALL /home/chippy/Download/ordinal * Installing

Re: [R] Ordinal Package Errors

2009-01-14 Thread Rick Bilonick
On Wed, 2009-01-14 at 16:38 +, Prof Brian Ripley wrote: You need to ask the author (as the posting guide asked you to). I'm tempted to not help further given the (almost complete) lack of cooperation of that author with R's recommendations, but note 'ordinal..so' in your log and look

[R] Errors running gam examples

2008-11-09 Thread Rick Bilonick
I'm running R 2.8.0 under Fedora 8 (32-bit). I installed the gam package. I can fit gam models, but I get error messages when I try to use step.gam and plot.gam, even for examples: library(gam) ?plot.gam data(gam.data) gam.object - gam(y ~ s(x,6) + z,data=gam.data) plot(gam.object,se=TRUE)

[R] Getting R and x11 to work

2008-06-05 Thread Rick Bilonick
I'm using Suse Linux Enterprise Desktop 10.2 (SP2) on an HP 2133 (x86) mini-notebook. (There apparently are a LOT of bugs in 10.1!) I downloaded R-base from the openSuse 10.2 repository and was (finally) able to install it (after installing blas and gcc-fortran). I can start an R session and do

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

2007-11-12 Thread Rick Bilonick
Using lmer in the lme4 package, you can compute the conditional variance-covariance matrix of the random effects using the bVar slot: 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

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

2007-11-12 Thread Rick Bilonick
On Mon, 2007-11-12 at 16:45 -0500, Doran, Harold wrote: No, don't reach into the bVar slot. Use the proper extractor function ranef() with postVar=T. There is no similar function for lme() -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Rick

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

2007-11-12 Thread Rick Bilonick
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 found a way: attr(ranef(lmer.13

[R] Confidence Intervals for Random Effect BLUP's

2007-11-09 Thread Rick Bilonick
I want to compute confidence intervals for the random effect estimates for each subject. From checking on postings, this is what I cobbled together using Orthodont data.frame as an example. There was some discussion of how to properly access lmer slots and bVar, but I'm not sure I understood. Is

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

2007-11-09 Thread Rick Bilonick
On Fri, 2007-11-09 at 10:01 -0800, Bert Gunter wrote: Ummm... Define: Confidence interval for BLUP . I know what a confidence interval for a parameter or function of parameters (which is what a predicted value is) is; but a BLUP is neither, so I don't get what a confidence interval for it

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

2007-11-09 Thread Rick Bilonick
On Fri, 2007-11-09 at 18:55 +, Prof Brian Ripley wrote: I think Bert's point is important: I picked up a student on it in a case study presentation on this week because I could think of three interpretations, none strictly confidence intervals. I think 'tolerance interval' is fairly

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

2007-10-24 Thread Rick Bilonick
I want to run a function in a loop and replace one of the arguments from a large list each time through the loop. If I was writing it out manually: myfunc(x=var1) myfunc(x=var2) etc. But I want to do this in a loop where x is replaced by a new name. Something like: for(i in vars) {

[R] Cannot Install rimage

2007-10-02 Thread Rick Bilonick
I'm trying to install rimage in R version 2.5.1 running on Fedora 6 (kernel 2.6.22.7-57.fc6 with the headers and gcc installed, along with fftw2 and libjpeg and headers): install.packages(rimage) Warning in install.packages(rimage) : argument 'lib' is missing: using '/usr/lib/R/library' trying