[R] multidimensional integration with adapt

2007-11-04 Thread Helene
Hello, I apologize for eventual double-posting. I am trying to integrate a 2-dimensional function that already calls the function adapt. More precisely, I am calling adapt(2,lower=c(-100,-100),upper=c(100,100),functn=function(s){1-exp(-50*Unc enteredGauss(c(-10,10,-10,10),60,s)}) w

[R] multidensional integration with adapt

2007-11-04 Thread Helene
Hello, I am trying to integrate a 2-dimensional function that already calls the function adapt. More precisely, I am calling adapt(2,lower=c(-100,-100),upper=c(100,100),functn=function(s){1-exp(-50*Unc enteredGauss(c(-10,10,-10,10),60,s)}) where UncenteredGauss is given by the followin

Re: [R] Sending E-Mail from R

2007-11-04 Thread Benilton Carvalho
On Nov 4, 2007, at 9:24 PM, (Ted Harding) wrote: > On 05-Nov-07 01:45:50, Josh Kalish wrote: >> Hi, >> Has anyone used R to send out an email via an SMTP server? >> I run R from a script and would like to be able to send out >> an email from there also. I could shell out to an email program, >>

Re: [R] perl module for R

2007-11-04 Thread Lux Zhang
On 04/11/2007, Andrew Perrin <[EMAIL PROTECTED]> wrote: > > http://www.omegahat.org/RSPerl/ > > -- > Andrew J Perrin - andrew_perrin (at) unc.edu - > http://perrin.socsci.unc.edu > Associate Professor of Sociology; Book Review Edit

Re: [R] Sending E-Mail from R

2007-11-04 Thread Ted Harding
On 05-Nov-07 01:45:50, Josh Kalish wrote: > Hi, > Has anyone used R to send out an email via an SMTP server? > I run R from a script and would like to be able to send out > an email from there also. I could shell out to an email program, > but was hoping maybe someone had a package for R or a smoo

Re: [R] Sending E-Mail from R

2007-11-04 Thread Gabor Grothendieck
Try SiteSearch("blat") from within R for examples using blat. Please read the posting guide and, in particular, give your platform next time. On Nov 4, 2007 8:45 PM, Josh Kalish <[EMAIL PROTECTED]> wrote: > Hi, > > > > Has anyone used R to send out an email via an SMTP server? I run R from

Re: [R] structure vs. matrix

2007-11-04 Thread Gabor Grothendieck
In general, you want to use as high level a construct as possible since then the details implementing the construct can be modified and improved without your code changing. This is really just an aspect of modularity. Thus matrix would be preferred. Actually since R allows access to the lower le

[R] Sending E-Mail from R

2007-11-04 Thread Josh Kalish
Hi, Has anyone used R to send out an email via an SMTP server? I run R from a script and would like to be able to send out an email from there also. I could shell out to an email program, but was hoping maybe someone had a package for R or a smooth way to do it. Thanks, Josh

Re: [R] structure vs. matrix

2007-11-04 Thread Edna Bell
I was reading "S Programming" and a set of the code used structure to set up a matrix. I just wondered which was better and why. On 11/4/07, Mark Difford <[EMAIL PROTECTED]> wrote: > > Hi Edna, > > >> When creating a matrix, is it better to use the structure function or the > >> matrix function

[R] R package "segmented"

2007-11-04 Thread Maura E Monville
Has anyone successfully used R package "segmented" ? I tried on my data and could not get the syntax right. In addition, I wonder whether it only operates with continuous variables as I deal with discrete variables. I was contacted by its author and, upon his request, I forwarded one of my files wh

Re: [R] install problem 2.6 mac os x 10.5

2007-11-04 Thread Chung-hong Chan
Refer to R Mac SIG list, the mpkg installer in 2.6 binary is broken in Leopard. Until an updated version is release, you can install R by install each package in "Packages" folder one by one. Hope it helps. Regards, CH On 11/4/07, Emilio Gagliardi <[EMAIL PROTECTED]> wrote: > Hi everyone, > > I

Re: [R] Help in error of mixed models

2007-11-04 Thread Andrew Robinson
Hi Bernardo, it may sound simple-minded, but it looks to me as though you have a problem with the model. Unfortunately it's difficult for us to diagnose the problem because you didn't send an example that we can reproduce. We don't have the "base" data. Based on my understanding of your code, y

Re: [R] install problem 2.6 mac os x 10.5

2007-11-04 Thread Matthew Keller
Hi Oscara and Emelio, You guys might want to do a search on the R mailing list - this issue has been discussed and I think there's a pretty simple work-around. Best, Matt On 11/4/07, Oscar Moreno <[EMAIL PROTECTED]> wrote: > I am having the very SAME problem. I removed Version 2.5.1 before > tr

Re: [R] unable to install package ff

2007-11-04 Thread Matthew Keller
Hi Adrian et al, Thanks for the advice. I haven't yet tried R.huge but will do so soon. I have used filehash and it is working well now. A note: for some reason, R objects were not being 'flushed' from memory on R 2.5.1 - even when removed and after an explicit call to gc(). Originally I thought

Re: [R] Problems with garch() function tseries package R 2.6.0

2007-11-04 Thread José Augusto Morais de Andrade J únior - JAMAJ
Hi, Yes, you are right. With Rterm it does work. With RGui, the newer version of R and tseries package, the function Garch() does not work anymore. Thanks. -- Cabeçalho original --- De: "Prof Brian Ripley" [EMAIL PROTECTED] Para: "José Augusto Morais de Andrade J únior - JAMAJ"

Re: [R] Problems with garch() function tseries package R 2.6.0

2007-11-04 Thread José Augusto Morais de Andrade J únior - JAMAJ
This is the sessionInfo() output immediatly before the hanging command garch(): sessionInfo() R version 2.6.0 (2007-10-03) i386-pc-mingw32 locale: LC_COLLATE=Portuguese_Brazil.1252;LC_CTYPE=Portuguese_Brazil.1252;LC_MONETARY=Portuguese_Brazil.1252;LC_NUMERIC=C;LC_TIME=Portuguese_Brazil.1252 at

[R] hierarchical mixed model

2007-11-04 Thread Irene Mantzouni
I would like to fit a 2-level mixed model: yit=a+a[i]+a[it] +(b+b[i]+b[it])*xit+eps[it] However, the variance of the second level components should depend on the group, i.e. sigma for a[it] and b[it] should be [i] specific. I do not know whether this is conceptually right in the mixed model

Re: [R] R CMD SHLIB gives error: bad value (generic) for -mtune= switch

2007-11-04 Thread McKay Curtis
Thanks Prof. Ripley. I removed that snippet from my path, and R CMD SHLIB works like a dream... Thanks again, McKay On Nov 4, 2007 2:32 AM, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > This is clearly a compiler configuration error, not an R error. I note > that you have c:\Rtools\MinGW\libex

[R] install problem 2.6 mac os x 10.5

2007-11-04 Thread Oscar Moreno
I am having the very SAME problem. I removed Version 2.5.1 before trying to install 2.6.0. Please help. TIA, Oscar __ 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.

Re: [R] Problems with garch() function tseries package R 2.6.0

2007-11-04 Thread José Augusto Morais de Andrade J únior - JAMAJ
Hi, I have just downloaded R 2.6.0 in other machine. The problem is the same. So, does anyone know how to make this function work again? For now, i will use my old version of R. Thanks in advance. José Augusto Jr. -- Cabeçalho original --- De: "Prof Brian Ripley" [EMAIL PROT

Re: [R] [Spam] Re: Problems with garch() function tseries package R 2.6.0

2007-11-04 Thread José Augusto Morais de Andrade J únior - JAMAJ
Hi, What other information should i post? I now have the two versions on the same machine. With R 2.3.1 and tseries 0.10-7 it works fine. With R 2.6.0 and tseries 0.10-11 the function garch doesnt work. The test was made with the same data, in the same machine. The maintainer did not say any

[R] Problem implementing adapt()

2007-11-04 Thread Sergey Goriatchev
that's the bvnpdf() code: bvnpdf <- function(xy, mu, Siginv, detSig){ f<-numeric() x <- xy[[1]] if(xy[[1]]>xy[[2]] & (1==2)) { f <- 0} else { v <- rbind(xy[1], xy[2]) e <- t(v-mu) %*% Siginv %*% (v-mu) f <- as.numeric(

[R] Problem implementing adapt()

2007-11-04 Thread Sergey Goriatchev
Hello, I am trying to compute rectangle probability of bivariate normal distribution with the following function: bvnrectangle <- function(mu, Sig, xmin, xmax, ymin, ymax){ library(adapt) Siginv <- solve(Sig) detSig <- det(Sig) areal <- adapt(ndim=2, lower=c(xmin,ymin), upper=c(xmax,ymax), minpts=

Re: [R] Why can repeated measures anova with within & between subjects design not be done if group sizes are unbalanced?

2007-11-04 Thread Charles C. Berry
'nowhere documented' ?? As the posting guide suggests, you could perform a search using > RSiteSearch("repeated measures", restric="functions") say. That generates an inventory of functions that pertain to repeated measures designs. And once you have browsed through that you might focus o

Re: [R] markov regime switching models

2007-11-04 Thread Liviu Andronic
Hello, On 10/30/07, valentina bonetti <[EMAIL PROTECTED]> wrote: > Hi, > > I am looking for a package to estimate regime switching models (states > following a markov chain). > I found packages for Hidden Markov Models but I am looking for something a > little different: In the HMM the conditional

[R] Why can repeated measures anova with within & between subjects design not be done if group sizes are unbalanced?

2007-11-04 Thread Gilbert G
Dear R people: I wish to switch from SPSS to R, but there is one particular type of ANOVA design that cannot be done in R. Or more likely, it can be done, but it is nowhere documented. The problem is typical for psychologists: You have a repeated measures design with different groups of subjects

Re: [R] structure vs. matrix

2007-11-04 Thread Prof Brian Ripley
On Sun, 4 Nov 2007, Duncan Murdoch wrote: > On 04/11/2007 1:28 AM, Edna Bell wrote: >> Hi R Gurus! >> >> When creating a matrix, is it better to use the structure function or >> the matrix function, please? > > I'd use the matrix() function. You can do it with structure(), and your > code might b

Re: [R] Problems with garch() function tseries package R 2.6.0

2007-11-04 Thread Prof Brian Ripley
Hmm, did you forget to tell us the minimal information required by the posting guide? And what did the package maintainer say? A guess: you are using Windows and doing this in Rgui. I believe 'tseries' does Fortran I/O, and this is a symptom we have seen before. The workaround is to use Rter

Re: [R] help on error message!!!

2007-11-04 Thread Bernardo Rangel Tura
On Sun, 2007-11-04 at 15:06 +0100, david csongor wrote: > Could anyone please help me with this one!!! > I am looping some simple stuff but when I enter an if - function I get this > error message: > > > Error in if (y < (2/(1 * pi) * sqrt(1 - ((x - 2)/1)^2))) c[i] = x else c[i] > = 0 : >

Re: [R] help on error message!!!

2007-11-04 Thread Duncan Murdoch
On 04/11/2007 9:06 AM, david csongor wrote: > Could anyone please help me with this one!!! > I am looping some simple stuff but when I enter an if - function I get this > error message: > > Error in if (y < (2/(1 * pi) * sqrt(1 - ((x - 2)/1)^2))) c[i] = x else c[i] > = 0 : > missing value

Re: [R] structure vs. matrix

2007-11-04 Thread Duncan Murdoch
On 04/11/2007 7:43 AM, Mark Difford wrote: > Hi Edna, > >>> When creating a matrix, is it better to use the structure function or the >>> matrix function...? > > I hope you have a huge (empty) jar in the kitchen, and that your pantry is > empty. > > R isn't too difficult, except if you're trying

Re: [R] structure vs. matrix

2007-11-04 Thread Duncan Murdoch
On 04/11/2007 1:28 AM, Edna Bell wrote: > Hi R Gurus! > > When creating a matrix, is it better to use the structure function or > the matrix function, please? I'd use the matrix() function. You can do it with structure(), and your code might be a little faster, but it will be harder to read, wi

[R] help on error message!!!

2007-11-04 Thread david csongor
Could anyone please help me with this one!!! I am looping some simple stuff but when I enter an if - function I get this error message: Error in if (y < (2/(1 * pi) * sqrt(1 - ((x - 2)/1)^2))) c[i] = x else c[i] = 0 : missing value where TRUE/FALSE needed WHAT does it mean, wha

Re: [R] Where to ask general stats questions?

2007-11-04 Thread Patrick Drechsler
francogrex <[EMAIL PROTECTED]> writes: > But does anyone know of a list/forum where we can ask help/questions > regarding general statistics problems (and where solutions are > offered possibly in the R language) . You could try the newsgroup sci.stat.math or sci.stat.consult. HTH Patrick -- W

[R] Problems with garch() function tseries package R 2.6.0

2007-11-04 Thread José Augusto Morais de Andrade J únior - JAMAJ
Hi all, I recently updated my to R 2.6.0 and tseries package ‘tseries’ version: 0.10-11. When i was using R Version 2.3.1 (2006-06-01) with tseries 'tseries' version: 0.10-7, the code > garch(dflnRCLC1) * ESTIMATION WITH ANALYTICAL GRADIENT * Call: garch(x = dflnRCLC1) Coefficient

Re: [R] structure vs. matrix

2007-11-04 Thread Mark Difford
Hi Edna, >> When creating a matrix, is it better to use the structure function or the >> matrix function...? I hope you have a huge (empty) jar in the kitchen, and that your pantry is empty. R isn't too difficult, except if you're trying to do stats (and don't know what you are doing --- though

[R] Where to ask general stats questions?

2007-11-04 Thread francogrex
I see that some people ask general statistics questions here sometimes. As I understand this list is only for "technical" help with R-related issues. But does anyone know of a list/forum where we can ask help/questions regarding general statistics problems (and where solutions are offered possibly

[R] Help in error of mixed models

2007-11-04 Thread Bernardo Rangel Tura
Hi R-masters I read the article: Bivariate analysis of sensitivity and specificity produces informative summary measures in diagnostic reviews. In this paper i proposed a bivariate mixed model and use SAS proc mixed to adjust the estimates. I thinks use R to make the same and try with this cod

Re: [R] mantel tests

2007-11-04 Thread Mark Difford
Hi Simon, >> This is a general statistics question so I'm sorry if its outside the >> field of r help. >> Could anyone please provide a starting point? This doesn't directly answer your question, but one thing I would do is get ade4 and look at Chessel's co-inertia analysis method. It's a very