[R] gam problem

Hello, People. I ahve been trying to use mle for normal distribution data set but always reporting an erroe on gam object. is there a solution to this Victor __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] scalassoc package

Many new things at http://www.cuddyvalley.org/psychoR/ The scalassoc package, which fits exponential distance association models to indicator matrices, it now at version 1.0.0. It seems to be robust and can analyze large examples easily. It is a major improvement (in speed and

[R] Unfound objects in function

A couple of my functions that were working last week seem to have been changed over the weekend and no longer work, but I can't understand why not: it seems that objects defined at the start of the function are not located further on in the function, when this worked fine before. An example

[R] get the parameter estimates from lm model

R-help, Sorry for posting the same question (almost) once again. I'm using a lm model for some data. The code is below: fitData - glm(pm ~ age + mw , family = binomial(link = logit), maxit = 30, data = inputModel1) age is a factor and mw is a continous variable. The summary function does not

[R] get the parameter estimates from lm model - was glm

R-help, Sorry for the mistake in my previous message. I have written lm but in the function call is glm. Everything should be glm. Sorry for the inconvenience. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] get the parameter estimates from lm model

Do you have a look at names(fitData) fitData$coef Best, Matthias -Ursprüngliche Nachricht- Von: [EMAIL PROTECTED] [mailto:r-help- [EMAIL PROTECTED] Im Auftrag von Luis Ridao Cruz Gesendet: Dienstag, 18. April 2006 10:23 An: r-help@stat.math.ethz.ch Betreff: [R] get the parameter

Re: [R] get the parameter estimates from lm model

Look at: ?contrasts ?contr.treatment I hope it helps. 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] Unfound objects in function

It does not say `unfound' it says `missing value'. See ?NA. Note it is the value and not the object that is said to be missing. Looks like on your data that one or more of the values you are summing is NA. On Tue, 18 Apr 2006, Colin Beale wrote: A couple of my functions that were working

[R] Extra \input command in Sweave created figures.

Dear list, I am using Sweave to generate figures and tables for my thesis. Now, I have run into the problem that all of a sudden, Sweave produces not only an \includegraphics command, but a \include command as well. This breakes LaTeX processing. This is an example: The code:

[R] adding lines to scatterplot3d

Hi to all, I usually use R for simple 2D graphs so now I cannot solve the following problem (which I suppose very simple too). I'm using scatterplot3d. I want to add dotted lines parallel to the 3 axes, in order to create a sort of cage in the 3D box (so to represent subdomains in my

Re: [R] using betareg: problems with anova and predict

On Mon, 17 Apr 2006, Paul Johnson wrote: Dear R-helpers: We have had fun using betareg to fit models with proportions as dependent variables. However, in the analysis of these models we found some wrinkles and don't know where is the best place to start looking for a fix. The problems we

Re: [R] Subset dataframe based on condition

Thanx Steve and Tony for your help. Sachin Tony Plate [EMAIL PROTECTED] wrote: Works OK for me: x - data.frame(a=10^(-2:7), b=10^(10:1)) subset(x, a 1) a b 4 1e+01 1e+07 5 1e+02 1e+06 6 1e+03 1e+05 7 1e+04 1e+04 8 1e+05 1e+03 9 1e+06 1e+02 10 1e+07 1e+01 subset(x, a 1 b a) a b 8

Re: [R] Extra \input command in Sweave created figures.

Fredrik Karlsson dargosch at gmail.com writes: I am using Sweave to generate figures and tables for my thesis. Now, I have run into the problem that all of a sudden, Sweave produces not only an \includegraphics command, but a \include command as well. This breakes LaTeX processing. Check

Re: [R] adding lines to scatterplot3d

andrea valle valle at di.unito.it writes: I'm using scatterplot3d. I want to add dotted lines parallel to the 3 axes, in order to create a sort of cage in the 3D box (so to represent subdomains in my dataspace). I saw that points3d can be used to draw points or lines to the existing

Re: [R] adding lines to scatterplot3d

Thanks a lot, I having it working. Is it possible to have dashed or dotted lines? I wasn't able to find a fitting type. Best -a- On 18 Apr 2006, at 14:42, Dieter Menne wrote: andrea valle valle at di.unito.it writes: I'm using scatterplot3d. I want to add dotted lines parallel to the 3

Re: [R] adding lines to scatterplot3d

andrea valle wrote: Thanks a lot, I having it working. Is it possible to have dashed or dotted lines? I wasn't able to find a fitting type. All as in standard plot functions: lty=dotted or lty=dashed should do the trick. Uwe Ligges Best -a- On 18 Apr 2006, at 14:42, Dieter

Re: [R] adding lines to scatterplot3d

Ah! It does. Thanks a lot. (I use R for batch processing, and I when I discover the right setup I don't change it for months... so I always forget basics...sorry) Best -a- On 18 Apr 2006, at 15:18, Uwe Ligges wrote: andrea valle wrote: Thanks a lot, I having it working. Is it possible

[R] how to change legend size in a figure

Dear All, I am producing a figure with many curves on it. How do I make the legends for all those curves smaller so that it can fit the figure itself? The commands I used for ploting are: plot(x1,y1,col=1,lty=1) lines(x2,y2,col=2,lty=2) ... legend(0.3,0.4,c(name1,name2,...),col=1:20,lty=1:20)

[R] lambda, uncertainty coefficient ( Somers D)

Dear colleagues in R, Has anybody implemented the 1) (Goodman Kruskal) lambda or the 2) (Thiel's) uncertainty coefficient statistics (in the asymmetric and symmetric forms), or is anyone aware that they might reside in some package? A search in the R archives does indicate that they are

Re: [R] how to change legend size in a figure

cex=0.8 (or whatever you like) as in legend(0.3,0.4,c(name1,name2,...),col=1:20,lty=1:20, cex=0.8) Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

[R] Nonlinear Regression model: Diagnostics

Hi, I am trying to run the following nonlinear regression model. nreg - nls(y ~ exp(-b*x), data = mydf, start = list(b = 0), alg = default, trace = TRUE) OUTPUT: 24619327 : 0 24593178 : 0.0001166910 24555219 : 0.0005019005 24521810 : 0.001341571 24500774 :

Re: [R] how to change legend size in a figure

Use a smaller cex within legend (e.g. cex=0.75) to adjust the size of the printed text. This controls the size of the legend. Also, you can play with xjust and yjust for finer control of how the legend box is justified at the specified position. -Christos -Original Message- From:

Re: [R] how to change legend size in a figure

probably the 'cex' and 'ncol' arguments of legend() will be useful here. 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

Re: [R] Unfound objects in function

Colin, Did you use correction = trans in Kinhom in the past? I think Ki$trans will contain a number of NAs. Try diff- sum (Ki$theo - Ki$trans, na.rm = TRUE) Peter Ehlers Colin Beale wrote: A couple of my functions that were working last week seem to have been changed over the weekend

Re: [R] MLE

vkatoma at cs.uct.ac.za writes: Hi , I want to compute the MLE for a simple sample of data, say 45,26,98,65,25,36,42,62,28,36,15,48,45, of which I obviously have the mean and the sd. Is there a way of calling the log normal and already diffrentiated formula other than entering the

Re: [R] lambda, uncertainty coefficient ( Somers D)

Antti Arppe wrote: Dear colleagues in R, Has anybody implemented the 1) (Goodman Kruskal) lambda or the 2) (Thiel's) uncertainty coefficient statistics (in the asymmetric and symmetric forms), or is anyone aware that they might reside in some package? A search in the R archives

[R] how funciton expression produces subscript

Dear R-users, does anyone know how funciton expression produces subscript? Best,Jing __ 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] how funciton expression produces subscript

See ?plotmath for the answer to your question and many more regarding math expression drawing. It might be worth printing this help page for future reference. -Christos -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jing Yang Sent: Tuesday, April 18,

[R] what happen?

Hi! I want to choose part of the data. But it has a question. attach(PIKO) PIKO[1:5,] code species dbh x y tag status branch 92 10142PIKO 38.9 6.05 12.81 10165 alive 0 109 10213PIKO 41.0 6.71 26.21 10202 alive 0 135 10222PIKO 48.3 18.98 21.28 10214

[R] Question about PLS regression

Hello useRs, I am new user to R and also statistics. Why predicted results in this example are different? Is the order of variables in X matrix important? library (pls) set.seed (1) Y1 - c(1,2,3,4,5,6,7,8,9,10) Y2 - c(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0) X1 - rnorm(10,sd=0.2) X2 -

Re: [R] what happen?

Look at: ?Comparison; you need: PIKO[status == snag, ] 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] Question about PLS regression

I am sorry. It was my fault. My example is wrong. I need also rearrange a validation data set too. But I have a sligthy different results with my real data. Where can the problem be? Andris Jankevics On Otrdiena, 18. Aprīlis 2006 17:55, Andris Jankevics wrote: Hello useRs, I am new user to

[R] 3D pie

Hi all, Is there a way to draw 3D pie with R (like excel does)? I know how to do it in 2D, just by using pie(something)... I know it isn't the best way to represent data, but people are sometimes more interested by the look and feel than by the accuracy of the results... If there

Re: [R] lambda, uncertainty coefficient ( Somers D)

On Tue, 2006-04-18 at 16:40 +0300, Antti Arppe wrote: Dear colleagues in R, Has anybody implemented the 1) (Goodman Kruskal) lambda or the 2) (Thiel's) uncertainty coefficient statistics (in the asymmetric and symmetric forms), or is anyone aware that they might reside in some

[R] R and ViM

Dear all, I'm starting to learn R, but I'm already programing for a while, using ViM as editor. Therefore I'd like to be able to use R together with ViM. My question now is, whether there are already people out there knowing how to do this in a similar easy way as with Emacs, and if those

Re: [R] R and ViM

Michael Graber wrote: Dear all, I'm starting to learn R, but I'm already programing for a while, using ViM as editor. Therefore I'd like to be able to use R together with ViM. My question now is, whether there are already people out there knowing how to do this in a similar easy way

[R] store levels in a string

Dear R-users, i need to store in a variable a string made from levels of a factor e.g. a-(a,a,b,b) af-factor(a) mylevels- ...a string with all the levels(af) Thanks -- DM __ R-help@stat.math.ethz.ch mailing list

Re: [R] lambda, uncertainty coefficient ( Somers D)

On Tue, 2006-04-18 at 10:30 -0500, Marc Schwartz (via MN) wrote: On Tue, 2006-04-18 at 16:40 +0300, Antti Arppe wrote: Dear colleagues in R, Has anybody implemented the 1) (Goodman Kruskal) lambda or the 2) (Thiel's) uncertainty coefficient statistics (in the

Re: [R] store levels in a string

On Tue, 2006-04-18 at 17:52 +0200, Daniele Medri wrote: Dear R-users, i need to store in a variable a string made from levels of a factor e.g. a-(a,a,b,b) The above should be: a - c(a,a,b,b) ^ af-factor(a) mylevels- ...a string with all the levels(af) Thanks -- DM

Re: [R] store levels in a string

Il giorno mar, 18/04/2006 alle 11.00 -0500, Marc Schwartz (via MN): mylevels- ...a string with all the levels(af) mylevels - levels(af) mylevels [1] a b I don't need to store these two levels, but a string with the values ab. Thanks Cheers -- DM

Re: [R] R and ViM

On 18-Apr-06 Sundar Dorai-Raj wrote: Michael Graber wrote: I'm starting to learn R, but I'm already programing for a while, using ViM as editor. Therefore I'd like to be able to use R together with ViM. [...] Try RSiteSearch('vim editor') at the R prompt. I got 28 hits. HTH, --sundar

Re: [R] store levels in a string

On Tue, 2006-04-18 at 18:07 +0200, Daniele Medri wrote: Il giorno mar, 18/04/2006 alle 11.00 -0500, Marc Schwartz (via MN): mylevels- ...a string with all the levels(af) mylevels - levels(af) mylevels [1] a b I don't need to store these two levels, but a string with the values

Re: [R] second try; writing user-defined GLM link function

Is this something that is fairly widely useful? If so, it is something that we could add to the stats package. (We ought in any case to make it easier to extend glm families, and this could be an example. If so, I'd like a reference for its use.) There is a good reason why logit is

Re: [R] R and ViM

Sundar == Sundar Dorai-Raj [EMAIL PROTECTED] on Tue, 18 Apr 2006 10:43:30 -0500 writes: Sundar Michael Graber wrote: Dear all, I'm starting to learn R, but I'm already programing for a while, using ViM as editor. Therefore I'd like to be able to use R together

[R] Any R/S course coming up this April in the bay area?

Please email me info on dates, etc __ 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] R and ViM

This was just posted this week: https://www.stat.math.ethz.ch/pipermail/r-help/2006-April/092457.html and I think there are some relevant macros on the vim web page I use vim but do not use anything special with it. It gives syntax highlighting out of the box. On 4/18/06, Michael Graber [EMAIL

Re: [R] R and ViM

yOn Tue, 18 Apr 2006, [EMAIL PROTECTED] wrote: On 18-Apr-06 Sundar Dorai-Raj wrote: Michael Graber wrote: I'm starting to learn R, but I'm already programing for a while, using ViM as editor. Therefore I'd like to be able to use R together with ViM. [...] Try RSiteSearch('vim editor') at

Re: [R] R and ViM

On 18-Apr-06 Prof Brian Ripley wrote: [...] RSiteSearch('vim editor') A search query has been submitted to http://search.r-project.org The results page should open in your browser shortly which it duly did. However -- and this is what's interesting -- I entered the command in R running

Re: [R] Specifying an appropriate error term in a hierarchical regression

On 17 Apr 2006 23:55:14 +0200, Peter Dalgaard [EMAIL PROTECTED] wrote: Confusing nesting and nesting... The / operator is designed to handle cases like this a b 1 1 1 2 2 1 2 2 in which the numbering of b only makes sense within a - no connection between b=1 when a=1 and when a=2

Re: [R] R debugging options

On Monday April 17 2006 21:08, Francisco J. Zagmutt wrote: RSiteSearch(debug) or RSiteSearch(debugging) will give you a lot or relevant information. I personally use library(debug) extensivelly and it should do all the taks you asked about. There is a nice article describing the debug

Re: [R] April Course***R/S-plus in San Francisco

Check out our upcoming 2-day R/S-plus Fundamentals and Programming Techniques in San Francisco: www.xlsolutions-corp.com/Rfund.htm San Franciso,April 27-28 Payment due AFTER the class Course Description: This two-day beginner to intermediate R/S-plus

[R] Compare two Power law or Exponential distributions

Dear All, I'd like to compare exponential or power-law distributions. To do so, people are often referred to the ks.test. However, I imagine ks.test wouldn't be as powerful as a test specifically designed for a distribution type. So my question is, is there a more specific test for each of

[R] Jeffries-Matusita distance

I've calculated as follows: matusita-function(f1,f2,lw,up){ fmatu-function(x) ( sqrt(f1(x)) - sqrt(f2(x)) )^2 sqrt (integrate( fmatu, lower=lw, upper=up, subdivisions=100)$value) } Josué Almansa Unitat de Recerca en Serveis Sanitaris Institut Municipal

Re: [R] 3D pie

You can see my opinion of 3D piecharts at http://www.burns-stat.com/pages/Tutor/spreadsheet_addiction.html Patrick Burns [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User) COMTE Guillaume wrote: Hi all, Is there a way to

[R] predict.nls confidence intervals

Hello- It has been several years since anyone has asked, so i am asking again- has anyone created a routine to estimate confidence intervals for predictions from nls models (ala Bates and Watts 1988)? Thanks - Alice Shelly [[alternative HTML version deleted]]

[R] installation of package randomForest failed

Hello I'd like to try out some functions in the package randomForest. Therefore, I did install this package. However, it is not possible to load the library, although I have R-Version 2.1.1 (i.e. later than 2.0.0). The commands I used and the Answers/Error from R is as follows:

Re: [R] what happen?

zhang jian wrote: * n=PIKO[status=snag,] Error in [.data.frame(PIKO, status = snag, ) : unused argument(s) (status ...) = is the assignment operator, you should use == as in: n=PIKO[status==snag,] __ LLama Gratis a

Re: [R] installation of package randomForest failed

Ruth, Does your R know about the library C://Programme//R//rw2011//library//randomForest ? I suspect that you want lib = C://Programme//R//rw2011//library in your install.packages() call. Or you could just leave the 'lib=' argument at its default. Peter Ehlers Ruth Meili wrote: Hello

Re: [R] store levels in a string

Il giorno mar, 18/04/2006 alle 11.29 -0500, Marc Schwartz (via MN): paste(levels(af), collapse = ) [1] ab See ?paste. :) more simple: toString(levels(x)) Thanks, Cheers -- DM __ R-help@stat.math.ethz.ch mailing list

[R] nls - nonlinear regression

Hi, I am trying to run the following nonlinear regression model. nreg - nls(y ~ exp(-b*x), data = mydf, start = list(b = 0), alg = default, trace = TRUE) OUTPUT: 24619327 : 0 24593178 : 0.0001166910 24555219 : 0.0005019005 24521810 : 0.001341571

Re: [R] store levels in a string

On Tue, 2006-04-18 at 20:53 +0200, Daniele Medri wrote: Il giorno mar, 18/04/2006 alle 11.29 -0500, Marc Schwartz (via MN): paste(levels(af), collapse = ) [1] ab See ?paste. :) more simple: toString(levels(x)) Thanks, Cheers -- DM Actually, you don't quite get the same

Re: [R] predict.nls confidence intervals

It is _just two days_ since someone asked: https://stat.ethz.ch/pipermail/r-help/2006-April/092437.html and that question was answered several times in the recent past, so I at least declined to answer it again. A quick search via RSiteSearch() shows several explanations

[R] questioin about cluster in R

Hi,All.Sorry for the group mail. I recently met a question and I have struggled on that for a while but failed to found the solution. I have a distance matrix as below. --- 03533936 350103251 331003049 93230035 365149

Re: [R] installation of package randomForest failed

Ruth Meili wrote: Hello I'd like to try out some functions in the package randomForest. Therefore, I did install this package. However, it is not possible to load the library, although I have R-Version 2.1.1 (i.e. later than 2.0.0). The commands I used and the Answers/Error from R is as

Re: [R] second try; writing user-defined GLM link function

A couple more comments: The null deviance is wrong here, as the code assumes that the link function maps constant vectors to constant vectors, which it does not here. You can circumvent that by setting an offset. Even setting dispersion = 1 I get slightly different se's. Here's a more robust

[R] R and tcltk

Hi people, I'm trying to run simple commands using tcltk in R but I got strange errors. Look at this: library(tcltk) Loading Tcl/Tk interface ... done tclRequire(Iwidgets) Tcl 4.0.1 tt - tktoplevel() tn - tkwidget(tt, iwidgets::tabnotebook) Erro em structure(.External(dotTclObjv, objv, PACKAGE

[R] R Reference Card (especially useful for Newbies)

Newbies (and others!) may find useful the R Reference Card made available by Tom Short and Rpad at http://www.rpad.org/Rpad/Rpad-refcard.pdf or through the Contributed link on CRAN (where some other reference cards are also linked). It categorizes and organizes a bunch of R's basic, most used

Re: [R] April Course***R/S-plus in San Francisco

Check out our upcoming 2-day R/S-plus Fundamentals and Programming Techniques in San Francisco: www.xlsolutions-corp.com/Rfund.htm San Franciso,April 27-28 Payment due AFTER the class Course Description: This two-day beginner to intermediate R/S-plus

Re: [R] R debugging options

Dear Larry, I'm not aware of an existing function that lists functions, but here's a simple solution: listFunctions - function(all.names=FALSE, envir=.GlobalEnv){ # all.names=TRUE: include names beginning with . # envir: environment to search Objects - objects(envir,

Re: [R] R and ViM

(Ted Harding) [EMAIL PROTECTED] writes: On 18-Apr-06 Prof Brian Ripley wrote: [...] RSiteSearch('vim editor') A search query has been submitted to http://search.r-project.org The results page should open in your browser shortly which it duly did. However -- and this is what's

[R] Different results under Windows and Linux

Hi everyone, disturbingly, a PCA yields different results when I run the analysis under Windows and Linux. I open the same workspace (which is on an independent partition) and perform the exact same function (by using the history file). Interestingly, PC1 is the same, but PC2 comes out totally

Re: [R] R and ViM

Hmm, how timely. I posted yesterday my solution to integrate R and vim. The message is in my sent box but I don't think it showed up in the list... Here it is again: Hi All, If you use vim to edit R code, you may be interested in this. I have put together a personalized syntax file, some code

Re: [R] R and ViM

Jose Quesada [EMAIL PROTECTED] writes: Hmm, how timely. I posted yesterday my solution to integrate R and vim. The message is in my sent box but I don't think it showed up in the list... Here it is again: Hi All, If you use vim to edit R code, you may be interested in this. I have

Re: [R] Different results under Windows and Linux

Thanks Bert, I did not describe the problem accurately enough. Actually, by now, I found out that floating-point imprecision is the reason for my confusion. The problem does not arise with the prcomp() function but rather with my data preprocessing: means=apply(dev,2,mean,na.rm=T) means

[R] How to efficiently extract or construct the residual covariance matrix from lme()?

Dear R-community, I'm trying to get the estimated residual covariance matrices from an lme object. If we write the model as: Y = X \beta + Z b + \epsilon and assume that b ~ N(0, P) and \epsilon ~ N(0, \Sigma), where P is non-diagonal and \Sigma might have correlation and weights components,

[R] about intersection set and union set

Hello sir: group1:1,2,3,4 group2:1,3,4,5 group3:2,4,8,9 .. ... group1000:9,3,8,2 I wanna get the intersection set and union set. I've tried command setdiffand union,but only two groups is permited.How can I deal with multi groups to find the intersection set and union set? Thanks a

Re: [R] about intersection set and union set

This was just discussed recently. See the following thread: http://tolstoy.newcastle.edu.au/R/help/06/04/24923.html On 4/18/06, XinMeng [EMAIL PROTECTED] wrote: Hello sir: group1:1,2,3,4 group2:1,3,4,5 group3:2,4,8,9 .. ... group1000:9,3,8,2 I wanna get the intersection set and union

[R] isoMDS and 0 distances

Hi, I'm trying to do a non-metric multidimensional scaling using isoMDS. However, I have some '0' distances in my data, and I'm not sure how to deal with them. I'd rather not drop rows from the original data, as I am comparing several datasets (morphology and molecular data) for the same

Re: [R] second try; writing user-defined GLM link function

It's probable that Roger's pdf that he provided is much more appropriate (or perhaps may have a good reference). However, if interested, the specific example used for my code, is from the following paper (and is now widely being used in many avian nesting studies, but I leave to others to

Re: [R] R and ViM

[Michael Graber] [...] I'd like to be able to use R together with ViM. [...] My question now is, whether there are already people out there knowing how to do this in a similar easy way as with Emacs [...] I've been an Emacs user for a very long time, and then, switched to Vim. See