Re: [R] Modality Test

Here is a function that tests for equality of however many distributions as you like: equaldist - function(...){ ## ... numeric sample vectors from the possibly different distributions to be tested ## returns TRUE only if the distributions are the same FALSE } ;-) -- Bert Gunter Genentech

[R] figure margins too large issue

Hi everyone, I updated R from 2.6.2 to 2.7.2 recently but keep getting the error figure margins too large when plot pictures which never happened to me when using 2.6.2. After googling and searching the mailing list, I still have no idea how to optimize it. One of my packages fixed bugs in 2.7.2

[R] test for a single variance

Dear R Gurus: Is there a test for a single variance available in R, please? Thanks, Edna Bell __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] New member with question on multiple comparisons in mixed effects models

Dear fellow R.users/.lovers, I am very new to both R and this list, so I hope you will be patient with me in the beginning if my enquiries are inappropriate/unclear. I am trying to perform some rather complex statistical modelling using mixed-effects models. I have, after a rather

[R] Package installation in Windows

As a sys admin, how do I install packages so that there is one common library for all users of the MS Windows computer, instead of the default individual location for each user? I've done this for Linux. Thanks, in advance, for your help. Mike __

[R] fBasics package: dnig

Hi, I am trying to calculate probability density of normal inverse gaussian distribution. I am using dnig function of fBasics package. However, I am getting following result. The density at x = 0.003042866 is: dnig(x= 0.003042866, alpha=5.184868, beta= 0.11841, delta= 0.06038513, mu=

[R] Genmod in SAS vs. glm in R

Hello, I have different results from these two softwares for a simple binomial GLM problem. From Genmod in SAS: LogLikelihood=-4.75, coeff(intercept)=-3.59, coeff(x)=0.95 From glm in R: LogLikelihood=-0.94, coeff(intercept)=-3.99, coeff(x)=1.36 Is there anyone tell me what I did wrong? Here

[R] probably easy methods question

Dear R Gurus: I want to look at the code for the t.test function. I did the following: t.test function (x, ...) UseMethod(t.test) environment: namespace:stats getAnywhere(t.test) A single object matching 't.test' was found It was found in the following places package:stats registered S3

Re: [R] probably easy methods question

On Tue, 9 Sep 2008, Edna Bell wrote: Dear R Gurus: I want to look at the code for the t.test function. I did the following: t.test function (x, ...) UseMethod(t.test) environment: namespace:stats getAnywhere(t.test) A single object matching 't.test' was found It was found in the following

Re: [R] Genmod in SAS vs. glm in R

For one thing your call to glm() is wrong --- didn't you notice the warning messages about ``non-integer #successes in a binomial glm!''? You need to do either: glm(r/k ~ x, family=binomial(link='cloglog'), data=bin_data, offset=log(y), weights=k) or: glm(cbind(r,k-r) ~ x,

[R] bootstrapping - number of items to replace is not a multiple of replacement length

Hello, I'm new to boostrapping and I'd need some help to understand the error message that pops up when I run my script. I have a data.frame with 73 lines and 21 column. I am running a stepwise regression to find the best model using the R function step. I apply bootstrapping to obtain model

Re: [R] NMDS and varimax rotation

At 9:54 PM +0200 9/9/08, Bernd Panassiti wrote: hello, subsequently to a NMDS analysis (performed with metaMDS or isoMDS) is it possible to rotate the axis through a varimax-rotation? Thanks in advance. Bernd Panassiti Bernd, Yes. The output of isoMDS is an object with points and

Re: [R] test for a single variance

Edna Bell edna.bell01 at gmail.com writes: Dear R Gurus: Is there a test for a single variance available in R, please? Thanks, Edna Bell Do you mean a test for homogeneity of variance? If so, try RSiteSearch(variance homogeneity) and see what you get ... Ben Bolker

Re: [R] Problem with starting and using R

The patched version (r46512) solves the problem!! Thanks! On Tue, Sep 9, 2008 at 12:20 AM, Thomas Lo [EMAIL PROTECTED] wrote: Hi Brian and Duncan, Many thanks for your responses. Setting the 'Current format' in 'Regional and language options' under Control Panel to English (Singapore)

Re: [R] densities with overlapping area of 0.35

Hi, I just findout that the KL Divengence depends on the variances of the distrubutions too and I'm trying to figureout how to do this. please help me. let's say X ~n(0,1) and Y~n(2*a, sigma^2). Is there any way I can calculate the sigma and a corresponding to an overlapping area of 0.35.

[R] Substitute Values in a Matrix?

Hi, I'm searching for a function to subistute Values in a Matrix to new Values. For example: old value new value 1.1 6 1.2 7 .. .. .. 1.9 14 2.0 15 and 2.1 15.5 2.2 16 .

[R] help: error handling in try

First time to post and searched archive for this problem with no clue. My version is 2.5.1. Below is a function to check if a given date is a valid date to a given date function object. It uses try (also tried tryCatch but with same problem). When given an invalid date, I am hoping try will

Re: [R] figure margins too large issue

on 09/09/2008 05:19 PM Hui-Yi Chu wrote: Hi everyone, I updated R from 2.6.2 to 2.7.2 recently but keep getting the error figure margins too large when plot pictures which never happened to me when using 2.6.2. After googling and searching the mailing list, I still have no idea how to

Re: [R] figure margins too large issue

It is not the default setting for par(mar) that changed. There are other device specific changes that occurred in 2.7.0, including, for some devices (eg. pdf and bitmap), the default device dimensions. See the bullets in the NEWS file under SIGNIFICANT USER-VISIBLE CHANGES and GRAPHICS CHANGES

Re: [R] bootstrapping - number of items to replace is not a multiple of replacement length

Hello, Your theta() function is returning different sets of coefficients depending on the results of step(). You'll need to add code to theta() to figure out which variables were selected, and store them into the right positions of a vector of length 20 (the apparent number of covariates you

Re: [R] Substitute Values in a Matrix?

on 09/09/2008 06:12 PM Chris82 wrote: Hi, I'm searching for a function to subistute Values in a Matrix to new Values. For example: old value new value 1.1 6 1.2 7 .. .. .. 1.9 14 2.0 15

Re: [R] extract variance components

Hi, Thanks for all. I now know how to extract the \sigma's. For the unbalanced model y_{ijk}=x\beta+\alpha_i+\beta_{ij}+e_{ijk} i=1,2,\dots,a, j=1,2,\dots,b_i, k=1,2,\dots,n_{ij} How can I extract the variance matrix $V$? The variance for the ith group is also of help. Suppose the ith group has

Re: [R] yahoo finance into R

thanks a lot. i thought it was more difficult. however i am dissapointed that i get so little data. thanks again thomas Peter Dalgaard wrote: thomastos wrote: Hi R, I am familiar with the basics of R. To learn more I would like how to get data from Yahoo!finance directly into R.

Re: [R] figure margins too large issue

Thank you Marc, I have checked two versions following your suggestions and the default setting of mar in both versions are same to each other as c(5,4,4,2)+0.1. Besides, I dont change any default settings about the graphics devices. For more info, following is my two sessions of R. *1.* R version

Re: [R] building a package that contains S4 classes and methods

Marie Pierre Sylvestre wrote: Hello R users, I am trying to make a my first package and I get an error that I can understand. The package is build out of three files (one for functions, 1 for s4 classes and 1 for s4 methods). Once I source them I run package.skeleton( name=TDC ) within a R

Re: [R] Genmod in SAS vs. glm in R

Rolf Turner wrote: For one thing your call to glm() is wrong --- didn't you notice the warning messages about ``non-integer #successes in a binomial glm!''? You need to do either: glm(r/k ~ x, family=binomial(link='cloglog'), data=bin_data, offset=log(y), weights=k) or: glm(cbind(r,k-r) ~

Re: [R] Problem with starting and using R

On Wed, 10 Sep 2008, Thomas Lo wrote: The patched version (r46512) solves the problem!!  Thanks! Thank you for confirming this. In the meantime I tracked down the exact problem. The standards for a locale name are 'language_country.encoding', as in 'en_GB.utf8' or 'English_United

<    1   2