Re: [R] [FORGED] Multiple Integrals

2016-01-06 Thread Rolf Turner
On 07/01/16 13:03, Luísa Freitas wrote: Dear R-helpers, I'm looking for a quick way to calculate triple integrals. I have tried something like this example: f <- function(x,y,z) dnorm(x)*dnorm(y)*dnorm(z) llim <- -Inf ulim <- Inf integrate(function(z) { sapply(z,function(z) { integr

[R] exact trend test

2016-01-06 Thread li li
Hi all, Is there an R function that does exact randomization trend test? For example, consider the 2 by 5 contingency table below: dose0dose 0.15dose 0.5dose 1.5dose 5 row margin Yes 43 4 5 8

[R] Multiple Integrals

2016-01-06 Thread Luísa Freitas
Dear R-helpers, I'm looking for a quick way to calculate triple integrals. I have tried something like this example: f <- function(x,y,z) dnorm(x)*dnorm(y)*dnorm(z) llim <- -Inf ulim <- Inf integrate(function(z) { sapply(z,function(z) { integrate(function(y) { sapply(y, functio

Re: [R] par(las=1) wirkt nicht

2016-01-06 Thread David Winsemius
> On Jan 6, 2016, at 10:27 AM, Hans Haußmann wrote: > > Hallo R-Experten, > > ich habe seit Kurzem einige Grafiken entwickelt mit R. > Wenn ich Grafikparameter wie z.B. "las" oder "tcl" im plot-Befehl definiere, > funktioniert das. Wenn ich aber im Skript "par(las=1, tcl=0.5)" schreibe, > ble

Re: [R] Mixed Beta Disrubutions

2016-01-06 Thread mesude bayrakci
Thank you Bob. I already had it worked. Best, M. On Wed, Jan 6, 2016 at 1:19 PM, Bob Aronoff wrote: > I am working to understand the same issues with my datasets. > Adapting Dr Zeileis' posting I have written a humble function. > It creates a two cluster beta mix on a vector of data. >

Re: [R] Amelia: Inputing Integers

2016-01-06 Thread Jim Lemon
Hi Lorenzo, Perhaps rounding the values in the imputed dataset? Jim On Thu, Jan 7, 2016 at 8:42 AM, Lorenzo Isella wrote: > Dear All, > I can provide a numerical example if needed. > My problem is the following: I am using amelia to input some missing > data in a dataset. > The problem is that

Re: [R] HELP - as.numeric changing column data

2016-01-06 Thread William Dunlap via R-help
By the way, here is an example where the advice in FAQ 7.10 (change the factor columns to numeric) would give incorrect results. The incorrect header setting in the call to read.table causes an extra row of non-numeric data to appear at the start of the imported data. > txt <- "ColA ColB\n101 1

[R] Amelia: Inputing Integers

2016-01-06 Thread Lorenzo Isella
Dear All, I can provide a numerical example if needed. My problem is the following: I am using amelia to input some missing data in a dataset. The problem is that some columns consist of integer numbers only and amelia inputs some real numbers with decimals. Is there a way to tell amelia that cert

Re: [R] package broom

2016-01-06 Thread Uwe Ligges
Yes, it depends on psych whoich depends on mnormt and all perfectly well declared, henc install.packages() should have done the job for you and installed the dependencies. Probably you did not use it for installing the package. Best, Uwe Ligges On 06.01.2016 19:08, Sarah Goslee wrote: That

Re: [R] HELP - as.numeric changing column data

2016-01-06 Thread William Dunlap via R-help
You may have read in your data incorrectly - a column you expected to be numeric was not recognized as such so it was read in a character and then converted to a 'factor'. FAQ 7.10 tells how to work around the problem https://cran.r-project.org/doc/FAQ/R-FAQ.html#How-do-I-convert-factors-to-numer

[R] HELP - as.numeric changing column data

2016-01-06 Thread Andy Schneider
Hi - I'm trying to plot some data and having a lot of trouble! I have a simple dataset consisting of two columns - income_per_capita and mass_beauty_value. When I read the data in and plot it, I get the attached plot Mass Beauty Non-Numeric:

[R] par(las=1) wirkt nicht

2016-01-06 Thread Hans Haußmann
Hallo R-Experten, ich habe seit Kurzem einige Grafiken entwickelt mit R. Wenn ich Grafikparameter wie z.B. "las" oder "tcl" im plot-Befehl definiere, funktioniert das. Wenn ich aber im Skript "par(las=1, tcl=0.5)" schreibe, bleibt das ohne Wirkung. Woran kann das liegen? Mit freundlichen Grüß

Re: [R] HELP - Hausman Test + systemfit

2016-01-06 Thread Arne Henningsen
Dear Kateryna The hausman() method for objects of class "systemfit" tests multiple-equation models estimated by 3SLS against multiple-equation models estimated by 2SLS (see documentation, e.g. [1]). It seems to me that this is not the type of Hausman test that you want to conduct. [1] http://www.

Re: [R] HELP - Hausman Test + systemfit

2016-01-06 Thread David Winsemius
> On Jan 6, 2016, at 9:44 AM, Jeff Newmiller wrote: > > To post appropriately, use less apologizing and more reading of the Posting > Guide that is mentioned in the footer of all messages on this list (including > this one). Note that HTML does not work well on the list, so be sure to turn >

Re: [R] Mixed Beta Disrubutions

2016-01-06 Thread Bob Aronoff
I am working to understand the same issues with my datasets. Adapting Dr Zeileis' posting I have written a humble function. It creates a two cluster beta mix on a vector of data. It seems to be working well on my datasets. You are welcome to try on yours. regards, Bob _ bi.modal.beta<

[R] A question about corAR1 and grouping

2016-01-06 Thread Saptarshi Guha
Hello, I was under the impression that in a call to gls, if i specify gls(.., cor=AR1(form= ~ 1 | subject)) and suppose there are N subjects, then i would have obtained estimates for N Phi (s) (the AR1 coefficient) - but I only get 1 estimate. I then checked "Linear Mixed-Effects Models Using R

Re: [R] package broom

2016-01-06 Thread Sarah Goslee
That means you didn't install the dependencies when you installed broom, and now need to install the package mnormt. No big deal, just do what the error message tells you. Sarah On Wed, Jan 6, 2016 at 12:57 PM, James Henson wrote: > Dear R community > > My version is R version 3.2.3. > The pack

Re: [R] package broom

2016-01-06 Thread Nordlund, Dan (DSHS/RDA)
So, you need to download and install the mnormt package. Dan Daniel Nordlund, PhD Research and Data Analysis Division Services & Enterprise Support Administration Washington State Department of Social and Health Services > -Original Message- > From: R-help [mailto:r-help-boun...@r-proje

Re: [R] package broom

2016-01-06 Thread Kevin E. Thorpe
On 01/06/2016 12:57 PM, James Henson wrote: Dear R community My version is R version 3.2.3. The package "broom" appears to install, but it will not load. The error message is below. library(broom) Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) : ther

[R] package broom

2016-01-06 Thread James Henson
Dear R community My version is R version 3.2.3. The package "broom" appears to install, but it will not load. The error message is below. > library(broom) Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) : there is no package called ‘mnormt’ Error: package

Re: [R] HELP - Hausman Test + systemfit

2016-01-06 Thread Jeff Newmiller
To post appropriately, use less apologizing and more reading of the Posting Guide that is mentioned in the footer of all messages on this list (including this one). Note that HTML does not work well on the list, so be sure to turn it off in your email program (at least when sending emails here).

Re: [R] hello everyone and happy new year

2016-01-06 Thread Ulrik Stervbo
Melt your data.frame library(reshape2) s3.melted <- melt(s3) and plot using something along the lines of ggplot(s3.melted,aes(x=time2,y=value, colour=variable) + geom_line() (not tested) Hope this helps Ulrik On Wed, 6 Jan 2016 at 17:01 Matthias Worni wrote: > I got the following problem

[R] HELP - Hausman Test + systemfit

2016-01-06 Thread Kateryna Riabchenko
Hello, I am not advanced in R (that is why my questions can sound stupid). I apologize for that in advance. But that is how my brain works - I need to ask questions to understand. I was searching the answer everywhere - without result. So I am asking You. *Given:* I have 8 regression Models: Mod

[R] hello everyone and happy new year

2016-01-06 Thread Matthias Worni
I got the following problem in R studio. Im trying to make a plot of different time series using different colours, which so far worked finde with the ggplot tool. However I did not manage to accually write a name to the correct layers. Below you can see the code that I wrote and also woks. Remembe

Re: [R] Bezier to line segments

2016-01-06 Thread Adrian Dușa
I just found the package "bezier". Trying to find the needle, I missed the haystack... On Wed, Jan 6, 2016 at 9:56 AM, Adrian Dușa wrote: > Dear All, > > I am interested into transforming Bezier curves (or general splines) to a > series of line segments. > For simplicity, the Bezier curves are e

[R] The R Journal, Volume 7, Issue 2

2016-01-06 Thread Bettina Gruen
Dear all, The latest issue of The R Journal is now available at http://journal.r-project.org/archive/2015-2/ Many thanks to all contributors. Regards, Bettina -- --- Bettina Grün Department of Applied Statistics JOHANNES KEPLER U