Re: [R] Update R 2.1.1-1 - 2.2.1 on UBUNTU

Maybe you need to spell cran.R-project.org in English, with a 'c'? On Wed, 1 Mar 2006, Florence Débarre wrote: I've installed Ubuntu 5.10 Breezy on my Mac. When using MacOSX, I used to work with R 2.2.1, but now, with Ubuntu, I can't download it with synaptic. I can only download R 2.1.1-1,

Re: [R] Help - lm, glm, aov results inconsistent with other stati stical package

BR == Ben Ridenhour [EMAIL PROTECTED] writes: BR Alright, I'll try to give some sample code. BR # create A with 2 levels - 2 and 4 O.k., let us add a set.seed(1) to make thinks reproducible. BR A-c(rep(2,times=30),rep(4,times=42)) BR # make A a factor BR

Re: [R] Update R 2.1.1-1 - 2.2.1 on UBUNTU

On Mer 1 mars 2006 8:58, Prof Brian Ripley a écrit : Maybe you need to spell cran.R-project.org in English, with a 'c'? Thanks, but... still doesn't work... W: Impossible de localiser la liste des paquets sources http://cran.R-project.org stable/ Packages

[R] Weighted networks and multigraphs

I would like to apply network measures (such as betweenness centrality, upper boundedness, etc.) to a weighted graph with non-integer weights, defined by a euclidean distance matrix. The package sna provides the measures that I want to use, but seems only to operate on binary graphs. I have read

Re: [R] Update R 2.1.1-1 - 2.2.1 on UBUNTU

In my opinion line in /etc/apt/sources.lst shoud be similar like this: deb http://cran.r-project.org/bin/linux/debian/stable ./ On Trešdiena, 1. Marts 2006 10:13, Florence Débarre wrote: On Mer 1 mars 2006 8:58, Prof Brian Ripley a écrit : Maybe you need to spell cran.R-project.org in

Re: [R] does svm have a CV to obtain the best cost parameter?

Thanks a lot Andy. I read that paper and followed the instructions, but met with a lot peculiarities: 1. In using tune function for svm, the best cost value turns out to be multi-peaks, and not with a single global peak. So I don't know which peak to follow in order to refine my search grid and

Re: [R] jpeg and pixels

On Tue, 28 Feb 2006, Erin Hodgess wrote: Dear R People: When using the jpeg function for plotting, is there a way to set the size in inches, please? There is an option for width and height in pixels, but not inches. Further to Marc Schwartz's answer, you can make use of the 'res' argument

Re: [R] Update R 2.1.1-1 - 2.2.1 on UBUNTU

On Mer 1 mars 2006 9:55, Andris Jankevics a écrit : In my opinion line in /etc/apt/sources.lst shoud be similar like this: deb http://cran.r-project.org/bin/linux/debian/stable ./ Doesn't work either... W: Impossible de localiser la liste des paquets sources http://cran.R-project.org ./

Re: [R] Help - lm, glm, aov results inconsistent with other stati stical package

Dear Ben Berwin is correct in his answer about different parameterizations. After changing the contrasts in R from treatment to sum options(contrasts = c(contr.sum, contr.poly )) test.lm-lm(y~A*x) summary(test.lm) I got similar results to JMP. Be careful in doing a correct interpretation of

Re: [R] heckit with a probit

Thank you both for your help. I've tried VGAM and it seems useful. Cheers David On 27/02/06, Robert Duval [EMAIL PROTECTED] wrote: I don't know if I understand your problem very well but the first reference that comes to mind is James Heckman, Dummy Endogenous Variables in a Simultaneous

Re: [R] Help - lm, glm, aov results inconsistent with other stati stical package

Perhaps it's worth mentioning that you can view and edit the contrasts for a particular factor by doing e.g. contrasts(A) contrasts(A) = contr.treatment(levels(A), base=2) # make the second category the base cheers David On 01/03/06, Christoph Buser [EMAIL PROTECTED] wrote: Dear Ben Berwin

[R] does ctest package exist?

R-help, I was searching the R archives for information on multiple comparisons and found that ctest package may contain something useful for my study. But there seems not to be such a package on CRAN. Anything I am missing? Thanks in advance __

Re: [R] Weighted networks and multigraphs

Dear Andrew, we've created a package called igraph, but the current version does not really support weighted graphs, so it's not much good for you. It can do weighted minimum spanning tree, but this is probably all. Also it can't handle multigraphs. I'm writing this email because if you end up

Re: [R] does ctest package exist?

On Wed, 2006-03-01 at 10:08 +, Luis Ridao Cruz wrote: R-help, I was searching the R archives for information on multiple comparisons and found that ctest package may contain something useful for my study. But there seems not to be such a package on CRAN. Anything I am missing?

Re: [R] does ctest package exist?

Luis Ridao Cruz [EMAIL PROTECTED] writes: R-help, I was searching the R archives for information on multiple comparisons and found that ctest package may contain something useful for my study. But there seems not to be such a package on CRAN. Anything I am missing? It was fused with the

Re: [R] R-help Digest, Vol 37, Issue 1

Mi nueva dirección de correo es: [EMAIL PROTECTED] New e-mail address: [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Bounding Box in pdf files

Mikael, I have encountered the same problem just recently and I was also stumbled upon it. I knew all that has been said by previous replies and therefore I was very confused. Now I have come up with the solution. You have to use option pdftex in either \usepackage[pdftex]{graphicx} or

[R] Maximum Likelihood Estimation

hi all, Can anyone tell me how to do Maximum Likelihood Estimation in R? Thanks in advance Arun -- [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] any more direct-search optimization method in R

Given that information, I think a genetic algorithm should probably do well with your problem. You may want to try the rgenoud package (R-GENetic Optimization Using Derivatives) which is on CRAN. For more information see: http://sekhon.berkeley.edu/rgenoud/ It works well for these kinds of

[R] interrupted time series analysis using ARIMA models

Hi R-users, I am using arima to fit a time series. Now I would like to include an intervention component It (0 before intervention, 1 after) using different types of impacts, that is, not only trying the simple abrupt permanent impact (yt = w It ) with the xreg option but also trying with a

Re: [R] Update R 2.1.1-1 - 2.2.1 on UBUNTU

Here is a copy of the relevant line in my /etc/apt/sources.list deb http://cran.R-project.org/bin/linux/debian sarge/ Hope this helps, Clint Florence Débarre wrote: I've installed Ubuntu 5.10 Breezy on my Mac. When using MacOSX, I used to work with R 2.2.1, but now, with Ubuntu, I can't

[R] Problems to get a ctree plot in a file via jpeg/png

Hello All, I am using library party and I have found a curious/strange behaviour when trying to save the output of a ctree in a file via jpeg/png command. If you use: library(party) airq - subset(airquality, !is.na(Ozone)) airct - ctree(Ozone ~ ., data = airq) plot(airct,

[R] mle.

hi all, suppose I have: r[i] ~ N(0, h[i]) h[i] = a + b*r[i-1] + c*h[i-1] for all i=2..n I want to get estimates of a, b, c by mle. Can you tell me how to do that? thanks in advance, Arun [[alternative HTML version deleted]] __

[R] C-Code in R (running Windows)

Hello, I need to include C (or C++)-Code in R. I'm searching for a tutorial, how to do this. It is only a very small algorith (one function). Using google, I read somethig about including dll-Files and so on (on windows). Is there no possibility to include a C-File directly?!? -- View this

[R] Class Method, S3 / S4

Hi I'm trying to create a class. However, I've some problems... I am used to program php. There I can create a new object and call a specific function of this object. In R, I hoped to create a class similar and then call the function like: Creating the Class:

[R] linear lists, creation, insertion, deletion, traversal *someone?*

Hi, In a second try I will ask this list to give me some useful pointers. Linear lists, as described e.g. by N.Wirth in Algorithms and Data Structures, seem not to be implemented in S/R, although in lisp we have cons, car, cdr. Nevertheless I want to implement an algorithm using such linear

Re: [R] C-Code in R (running Windows)

The simplest way to add C (C++) code in R, that I know is by creating a package ( even if you are not planning on distributing it). See 'Writing R Extensions' manual (http://cran.r-project.org/manuals.html) for details. Download source codes of packages for more examples, many of them will have .c

Re: [R] mice library / survival analysis

Dear Brian, The mice() function returns an object of class mids that contains information for constructing the multiply-imputed data sets; you can extract these data sets using complete(). The general idea is then to perform your analysis in parallel on the completed data sets. Take a look at

Re: [R] Help - lm, glm, aov results inconsistent with other stati stical package

Dear Prof. Ripley Please appologize. I was a little too short in my first answer. My guess is that the line equations have been produced by Ben (since JMP does not give me that line equation as output), using the coefficients from JMP with the wrong assumption that JMP uses the constraint that

[R] matrix matching NA

Dear R list members, I am seeking for an elegant way for matching NA elements of a matrix. Everthing I tried, the result was a vector. Here ist an example with vectors that works: a and b are vectors of same lenght representing pairs of data (a[i] b[i]) with some NA elements: ## R Code a -

Re: [R] does svm have a CV to obtain the best cost parameter?

Do you know that there are (pseudo-)randomness involved in CV? Even if you fix the parameters and run multiple times, you're going to get different answers, let alone changing parameters each time. Hoping to narrow the optimal parameters down to that fine a resolution is generally not realistic.

[R] Drop1 and weights

Hi, If I used drop1 in a weighted lm fit, it seems to ignore the weights in the AIC calculation of the dropped terms, see the example below. Can this be right? Yan library(car) unweighted.model - lm(trSex ~ (river+length +depth)^2- length:depth, dno2)

[R] SOLVED Re: Update R 2.1.1-1 - 2.2.1 on UBUNTU

deb http://cran.R-project.org/bin/linux/debian sarge/ This line worked, after having typed apt-get update in a terminal. worked means that synaptic didn't crash. But I still wasn't able to download the newest version of R. I've done it in a coarse way : I've downloaded the .tar file. I had to

Re: [R] C-Code in R (running Windows)

The simplest interface is to use .C() in R to call the C function. The details (with examples) are in the Writing R Extensions manual that comes with every version of R. To compile and load the .dll/.so into R, do: - Type R CMD SHLIB file1.c file2.c (at the command prompt, not the R prompt)

Re: [R] linear lists, creation, insertion, deletion, traversal *someone?*

On 3/1/2006 8:09 AM, Christian Hoffmann wrote: Hi, In a second try I will ask this list to give me some useful pointers. Linear lists, as described e.g. by N.Wirth in Algorithms and Data Structures, seem not to be implemented in S/R, although in lisp we have cons, car, cdr. They are in

Re: [R] matrix matching NA

probably you'll find ?complete.cases() useful Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015 Web:

Re: [R] matrix matching NA

Dear list members, Berwin gave me the right hint, it works! Look at this: ## R Code a - (1:5) a[2] - NA b - (6:10) b[3] - NA a 1 NA 3 4 5 b 6 7 NA 9 10 matrix - array(1:5*2, dim=c(5,2)) matrix[,1] - a matrix[,2] - b matrix[ !apply(is.na(matrix), 1, any),] [,1] [,2] [1,]

Re: [R] matrix matching NA

Is this what you are looking for? m - matrix(1:10, 5, 2) m[3, 1] - m[5, 1] - m[2, 2] - NA m [,1] [,2] [1,]16 [2,]2 NA [3,] NA8 [4,]49 [5,] NA 10 m[complete.cases(m),] [,1] [,2] [1,]16 [2,]49 Andy From: Fabian Lienert Dear R list

Re: [R] C-Code in R (running Windows)

I got it!!! As Andy said, I generated an dll-File (using Dev-C++, it was quite easy) and imported some R-header. This could be opened in R with dyn.load(dll-file) und executed with .C()!! Thanks for your help!! PS: I've taken the Files of my (Dev C++ -)project behind: dllmain.c: #include

Re: [R] linear lists, creation, insertion, deletion, traversal *someone?*

Christian Hoffmann wrote: Hi, In a second try I will ask this list to give me some useful pointers. Linear lists, as described e.g. by N.Wirth in Algorithms and Data Structures, seem not to be implemented in S/R, although in lisp we have cons, car, cdr. This has been discussed on the

Re: [R] linear lists, creation, insertion, deletion , traversal *someone?*

Christian Hoffmann christian.hoffmann at wsl.ch writes: Hi, In a second try I will ask this list to give me some useful pointers. [There is a general problem with the list, because of its volunteer nature: easy questions and questions that catch peoples' attention get answered. Others

[R] glm binomial with zero proportions

Hello. I must fit a logistic regression to data in the form of proportions, but in which some of the proportions are zero. I therefore cannot use the glm function with a binomial link since the link function is not defined for p=0 or 1. What other solutions are available? Any references to

Re: [R] Maximum Likelihood Estimation

Arun Kumar Saha arun.kumar.saha at gmail.com writes: hi all, Can anyone tell me how to do Maximum Likelihood Estimation in R? Unfortunately this question is ***way*** too vague for us to answer adequately. A short answer is that R provides general-purpose minimization functions

Re: [R] mle.

Arun Kumar Saha arun.kumar.saha at gmail.com writes: hi all, suppose I have: r[i] ~ N(0, h[i]) h[i] = a + b*r[i-1] + c*h[i-1] for all i=2..n I want to get estimates of a, b, c by mle. Can you tell me how to do that? (1) can you convince us this isn't a homework

[R] Problems to get a ctree plot (library party) in a file via jpeg/png

Hello All, I am using library party and I have found a curious/strange behaviour when trying to save the output of a ctree in a file via jpeg/png command. If you use: library(party) airq - subset(airquality, !is.na(Ozone)) airct - ctree(Ozone ~ ., data = airq) plot(airct,

[R] lattice-Internal

Hi, The functions prepanel.default.bwplot() and lpretty() are not running in Deepayan's barley example concerning vertical bars with the lattice function barchart(). Why, is there a restricted use for the package lattice-Internal? Urs Simmen mailto:[EMAIL PROTECTED]

Re: [R] glm binomial with zero proportions

On Wed, 1 Mar 2006, Bill Shipley wrote: Hello. I must fit a logistic regression to data in the form of proportions, but in which some of the proportions are zero. I therefore cannot use the glm function with a binomial link since the link function is not defined for p=0 or 1. It is a

Re: [R] ex-Gaussian survival distribution

On Tue, 28 Feb 2006, Jonathan Williams wrote: Dear R-Helpers, I am hoping to perform survival analyses using the ex-Gaussian distribution. I understand that the ex-Gaussian is a convolution of exponential and Gaussian distributions for survival data. I checked the survreg.distributions

Re: [R] linear lists, creation, insertion, deletion, traversal *someone?*

Also check out ?append and the after= argument, in particular. On 3/1/06, Ben Bolker [EMAIL PROTECTED] wrote: Christian Hoffmann christian.hoffmann at wsl.ch writes: Hi, In a second try I will ask this list to give me some useful pointers. [There is a general problem with the list,

[R] How to do a proc summary in R?

Hi, I'm a SAS user trying to convert myself to R but I still have problems with some pretty simple commands. First I wanted to add up a number of red-tailed hawks seen per day (julian day) per year. So I tried: RTyrday - tapply(RThr,list(year,julian),sum) And then I tried the following

[R] inconsistency between anova() and summary() of glmmPQL

Dear All, Could anyone explain me how it is possible that one factor in a glmmPQL model is non-significant according to the anova() function, whereas it turns out to be significant (or at least some of its levels differ significantly from some other levels) according to the summary()

Re: [R] How to do a proc summary in R?

Emilie Berthiaume wrote: Hi, I'm a SAS user trying to convert myself to R but I still have problems with some pretty simple commands. First I wanted to add up a number of red-tailed hawks seen per day ( julian day) per year. So I tried: RTyrday - tapply(RThr,list(year,julian),sum)

[R] Using read.table() with Excel data

Help I apologize if this message posts twice, but I do not think my last attempt went through. I am trying to use the read.table() function to read in a massive alphanumeric data set. I thought I might scan in each column individually and do a cbind() to reconstitute the array in R, so

Re: [R] inconsistency between anova() and summary() of glmmPQL

To quote one of my professors, it usually doesn't make sense to ask questions like `Is variable X significant?' (Or, sort of more formally, testing H0: beta_j = 0 vs. H1: beta_j != 0.) If `X' is the _only_ variable you will ever consider, then the question can make sense. Otherwise, you need

[R] Width of bars in barplot2

I'm using barplot2 to plot some data. Is there any way to determine the width of the bars in the generated plot? I know that barplot2 returns a list of the coordinates of the center of each bar, but since there is some white space between each bar, I don't know how to get the width of each bar.

Re: [R] Maximally independent variables

In case others are interested I did get a reply offlist regarding the escouf function in the pastecs package. See: library(pastecs) ?escouf Also see pages 47-52 of system.file(doc/pastecs.pdf, package = pastecs) (in French). On 3/1/06, Gabor Grothendieck [EMAIL PROTECTED] wrote:

Re: [R] How to do a proc summary in R?

Emilie Berthiaume wrote: Hi, I'm a SAS user trying to convert myself to R but I still have problems with some pretty simple commands. First I wanted to add up a number of red-tailed hawks seen per day (julian day) per year. So I tried: RTyrday - tapply(RThr,list(year,julian),sum)

[R] Selecting points

Dear all, I have point data from an irregular survey. Starting with a random point, I would like to select as many other points as possible, so long as the minimum inter-point distance is never below a given threshold. I would appreciate help in finding a nice solution. Thankyou in advance.

Re: [R] Width of bars in barplot2

On Wed, 2006-03-01 at 12:01 -0500, Jamieson Cobleigh wrote: I'm using barplot2 to plot some data. Is there any way to determine the width of the bars in the generated plot? I know that barplot2 returns a list of the coordinates of the center of each bar, but since there is some white space

[R] a strange problem with integrate()

Dear all, I am stuck on the following problem with integrate(). I have been out of luck using RSiteSearch().. My function is g2-function(b,theta,xi,yi,sigma2){ xi-cbind(1,xi) eta-drop(xi%*%theta) num-exp((eta + rep(b,length(eta)))*yi) den- 1 + exp(eta +

Re: [R] Selecting points

I just asked a similar question yesterday except I was using correlation-based distances rather than Euclidean. Check out that thread: https://www.stat.math.ethz.ch/pipermail/r-help/2006-March/087782.html On 3/1/06, Matt Sakals [EMAIL PROTECTED] wrote: Dear all, I have point data from an

Re: [R] lattice-Internal

Hi Urs, Are you referring to the ?barley help topic example that calls dotplot() ?. I just tried it via example(barley) and cut-and-paste it and it ran with now complaints on my system: R 2.2.1, Windows 2000 SP4, lattice version 0.12-11 Both the functions you mention seem to be there:

Re: [R] Using read.table() with Excel data

On 3/1/06 11:31 AM, Nealy, David [EMAIL PROTECTED] wrote: Help I apologize if this message posts twice, but I do not think my last attempt went through. I am trying to use the read.table() function to read in a massive alphanumeric data set. I thought I might scan in each

Re: [R] a strange problem with integrate()

vito muggeo wrote: Dear all, I am stuck on the following problem with integrate(). I have been out of luck using RSiteSearch().. My function is g2-function(b,theta,xi,yi,sigma2){ xi-cbind(1,xi) eta-drop(xi%*%theta) num-exp((eta + rep(b,length(eta)))*yi)

Re: [R] Width of bars in barplot2

Thanks. That was all I needed. Jamie On 3/1/06, Marc Schwartz (via MN) [EMAIL PROTECTED] wrote: Unless you have modified the 'width' argument, the default width is 1. Thus the sides of the bars are the centers +/- 0.5. If you modified the width argument, then the widths are set to the

Re: [R] Class Method, S3 / S4

The object paradigm in R is different; I like to think of methods as orthogonal to classes, instead of nested within them. Probably what you want to try is # a generic function, to operate on different classes setGeneric(test, function( obj ) standardGeneric( test )) # a class, containing data

Re: [R] question about lmer--different answers from different versions of R?

Yes, I believe that lmer has been under active development during the past several months, and an intentional change in the code might have produced change you reported. In particular lmer was changed to use nlminb instead of optim last year, and this could easily generate the

Re: [R] inconsistency between anova() and summary() of glmmPQL

My apologies: I got the descriptions of tests in summary() and anova() backward. Cheers, Andy From: I.Szentirmai Dear Andy, Thanks a lot for clarifying this for me. However, to me it seems that anova is the one that provides the results of a sequential test since whether a factor

Re: [R] lattice-Internal

On 3/1/06, Urs Simmen [EMAIL PROTECTED] wrote: Hi, The functions prepanel.default.bwplot() and lpretty() are not running in Deepayan's barley example concerning vertical bars with the lattice function barchart(). What example? Please be more specific. Why, is there a restricted use for

Re: [R] does svm have a CV to obtain the best cost parameter?

Can svm/tune in R do an automatic grid search with increasingly finer grid and smaller focus for me? On 3/1/06, Liaw, Andy [EMAIL PROTECTED] wrote: Do you know that there are (pseudo-)randomness involved in CV? Even if you fix the parameters and run multiple times, you're going to get

[R] Wavelet ordenates

Hello, does someone knows how can I get the ordenates of a father wavelet using the package wavethresh? And using other package? Any hints? For example, I'd like knowing what's the value of \phi{j,k}(t), that is, the value of a father wavelet (from some family) at scale j, shift k, point t.

[R] help regarding permieter/boundary detection in spatial data set

I have a large data set that consists of x,y coordinates and number of plants counted at each point. The sampling grid consists of contiguous 1 m^2 quadrats placed over a patch of plants such that the patch is competely within the grid and the perimeter quadrats straddle the patch edge. Each 1

[R] fixed variance of lmer object

Dear expRts, could anybody nudge me toward some documentation about how to extract the variance-covariance matrix of the fixed parameters from an lmer object? Best Johannes -- Johannes HüsingThere is something fascinating about science. [EMAIL PROTECTED] One gets such wholesale

[R] library file for R's nmath routines

Hi, I am wondering where the library file for R's nmath routines are? Doing a search on libR gave me the following: /usr/lib/libRKC16.so.1.2.0 /usr/lib/libRKC.so.1.2.0 /usr/lib/R/lib/libRlapack.so /usr/lib/R/lib/libR.so /usr/lib/libRKC.so.1 /usr/lib/libRKC16.so.1 None of these have the

[R] Optimization problem: selecting independent rows to maximize the mean

Dear R community, I have a dataframe with 500,000 rows and 102 columns. The rows represent spatial polygons, some of which overlap others (i.e., not all rows are independent of each other). Given a particular row, the first column contains a unique RowID. The second column contains the Variable

Re: [R] Optimization problem: selecting independent rows to maximizethe mean

This sounds either easy via a greedy algorithm or NP-hard. Moreover, it is not clear to me that 1) A subset of 500 indpendent rows exists, where I presume independent means pairwise nonoverlapping; 2) That the mean and sd can be simultaneously optimized as you describe-- what if the subset with

[R] Book: Multilevel Modeling in R ETA?

Hi R folks (Dr. Bates in particular), In August 2005, Dr. Bates mentioned that the documentation for lme4 will be in the form of a book with the working title 'Multilevel Modeling in R' and I'm just wondering if there is an estimated date of publication or if it's still a long way off. The

Re: [R] Optimization problem: selecting independent rows to maximize the mean

Package lpSolve might help. On 3/1/06, Mark [EMAIL PROTECTED] wrote: Dear R community, I have a dataframe with 500,000 rows and 102 columns. The rows represent spatial polygons, some of which overlap others (i.e., not all rows are independent of each other). Given a particular row, the

Re: [R] Maximally independent variables

Gabor Grothendieck ggrothendieck at gmail.com writes: That's basically what I already do but what I was wondering was if there were any other approaches such as connections with clustering, PCA, that have already been developed in R that might be applicable. Have you considered finding

[R] postscript bounding box in trellis/lattice plot is wrong ?

Hi, a problem involving postscript bounding boxes: I'm composing three scatterplots into a single figure, postsript for publication. The individual scatterplots should be square, so the overall figure should have a roughly 1:3 sort of aspect ratio. By default however, the overall figure comes

[R] parameterization of gnls model

I am working with estimates of vegetation height derived from radar data. We have a nonlinear model to correct these estimates for errors associated with viewing geometry. I am trying to estimate a single parameter in this model while accounting for spatial (spherical structure) autocorrelation.

Re: [R] postscript bounding box in trellis/lattice plot is wrong ?

On 3/1/06, context grey [EMAIL PROTECTED] wrote: Hi, a problem involving postscript bounding boxes: I'm composing three scatterplots into a single figure, postsript for publication. The individual scatterplots should be square, so the overall figure should have a roughly 1:3 sort of

[R] installing R on OS 10.4.5

Dear all, I'm trying to install R on Mac OS 10.4.5, on a new MacBook Pro (Intel processor). Everything is standard at the moment in my setup, except I switched to tcsh and customized my environment variables, aliases, etc in the tcsh directory. I installed the developer tools directly off the

Re: [R] Weighted networks and multigraphs

Hi, Some implementations are available in the packages graph and RBGL from Bioconductor (www.bioconductor.org). I think you will find that they do what you would like, but let the maintainers know if not. best wishes Robert Wilson, Andrew wrote: I would like to apply network

Re: [R] SOLVED Re: Update R 2.1.1-1 - 2.2.1 on UBUNTU

Florence, On 1 March 2006 at 14:56, Florence Débarre wrote: | deb http://cran.R-project.org/bin/linux/debian sarge/ | This line worked, after having typed | apt-get update | in a terminal. worked means that synaptic didn't crash. But I still | wasn't able to download the newest version of R. I

Re: [R] obtaining means/SD after fitting a mixed model

Have you consulted Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer)? The MASS function glmmPQL translates the problem into approximate gaussian terms and passes the translation into lme. Therefore, I would expect much of what appears in Pinheiro and Bates

[R] [R-pkgs] New Quantian releases with almost all of CRAN / BioC

[ I hope this is not considered off-topic for the list: A new Quantian release was just announced, and it contains even more R. As before, we have R (now 2.2.1), ESS, Ggobi, Rpy, several BioInfo tools -- but now we also have RSPerl, JGR, rJava and more. A full sync to CRAN and BioC as

[R] How can I use r-cran-lmtest?

Hi, I installed the package r-cran-lmtest in a Debian Sid but I can't use it. I typed lmtest but nothing occur. Any help? Thanks. -- Sávio Martins Ramos - Arquiteto Rio de Janeiro ICQ 174972645 Pirataria não! Seja livre: Linux http://www.debian.org

Re: [R] How can I use r-cran-lmtest?

At 12:28 AM 3/2/2006, Savio Ramos wrote: Hi, I installed the package r-cran-lmtest in a Debian Sid but I can't use it. I typed lmtest but nothing occur. Any help? Thanks. Hi Savio! The Lmtest is a package for Testing Linear Regression Models. This is a organized group of rotines, data

Re: [R] How can I use r-cran-lmtest?

On 2 March 2006 at 00:28, Savio Ramos wrote: | I installed the package r-cran-lmtest in a Debian Sid but I can't use it. I | typed lmtest but nothing occur. You need to issue a library(foo) command for every package foo you want to use. So try library(lmtest) before you do

[R] Autoregressive Model with Independent Variable

Hey, all, I may just be missing something, but I'm trying to construct a temporal autoregression with an independant variable other than just what is happened at a previous point in time. So, the model structure would be something like y(t)=b0+b1*y(t-1)+b2*y(t-2)...+a*x(t) I'm even

Re: [R] Autoregressive Model with Independent Variable

On 1 March 2006 at 20:06, Jarrett Byrnes wrote: | Hey, all, I may just be missing something, but I'm trying to construct | a temporal autoregression with an independant variable other than just | what is happened at a previous point in time. So, the model structure | would be something like |

Re: [R] Autoregressive Model with Independent Variable

It looks like just simple ARMA model; there is no visible I-part in Jarrett's specification. Unless it's hidden in the '...' part :) Janusz. On Wed, 1 Mar 2006, Dirk Eddelbuettel wrote: On 1 March 2006 at 20:06, Jarrett Byrnes wrote: | Hey, all, I may just be missing something, but I'm trying

[R] Overfitting in snr (ASSIST)

Dear R-Help group, I'm trying to fit some smoothing spline models using the snr function in the ASSIST library. My model is the following: y~log(exp(a0-a1*x)+exp(b0-b1(x))) where b1(x) is a nonparametric function. It's basically a biexponential decay, with wiggliness in 'phase 2' decay. I

Re: [R] Autoregressive Model with Independent Variable

On Wed, 1 Mar 2006, Janusz Kawczak wrote: It looks like just simple ARMA model; there is no visible I-part in Jarrett's specification. Unless it's hidden in the '...' part :) If it were an ARMA model, x would be (unknown) white noise but we were told its [integer] values. So in no useful

Re: [R] library file for R's nmath routines

On Wed, 1 Mar 2006, someone claiming to be `Globe Trotter' wrote: Hi, I am wondering where the library file for R's nmath routines are? Doing a search on libR gave me the following: /usr/lib/libRKC16.so.1.2.0 /usr/lib/libRKC.so.1.2.0 /usr/lib/R/lib/libRlapack.so /usr/lib/R/lib/libR.so

[R] Doubly Non-Central F-Distribution

Dear Professor I have read your questions in the website on Doubly non-central F-distribution. I am looking for source code for evaluating this function now. I have matlab code but it's only accurate up to the 4th decimal point. Dataplot is more accurate, but it is not user friendly as you