[R] orca binary?

2004-05-28 Thread Laura Holt
Dear R People: Is there a binary version for rorca, please? Thanks, Laura mailto: [EMAIL PROTECTED] R Windows 1.9.0 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] gauss.hermite?

2004-05-28 Thread Dimitris Rizopoulos
Dear Spencer, In Dr. Gray's course notes for Advanced Statistical Computing (Appendix of Chapter 6) there are some functions for computing abscissas and corresponding weights for several Gaussian integration rules: http://icommons.harvard.edu/~hsph-bio248-01/Lecture_Notes/ I hope this helps.

[R] distance in the function kmeans

2004-05-28 Thread n\.bouget
Hi, I want to know which distance is using in the function kmeans and if we can change this distance. Indeed, in the function pam, we can put a distance matrix in parameter (by the line pam-pam(dist(matrixdata),k=7) ) but we can't do it in the function kmeans, we have to put the matrix of data

[R] Building packages and Rcmd check

2004-05-28 Thread Laura Holt
Hello yet again! I'm trying to build a package. The Rcmd build works fine. But the Rcmd check produces an error because I don't have latex on my laptop. Any suggestions, please? R Windows 1.9.0 thanks again, Laura _ Learn to

Re: [R] privileged slots

2004-05-28 Thread Matthias Kohl
Martin Maechler schrieb: Matthias == Matthias Kohl [EMAIL PROTECTED] on Thu, 27 May 2004 14:01:51 +0100 writes: Matthias Hi all, in the help for RClassUtils I found the Matthias expression privileged slots in function Matthias checkSlotAssignment with the explanation:

Re: [R] Building packages and Rcmd check

2004-05-28 Thread Uwe Ligges
Laura Holt wrote: Hello yet again! I'm trying to build a package. The Rcmd build works fine. But the Rcmd check produces an error because I don't have latex on my laptop. Any suggestions, please? Rcmd check --help tells you that you can use the option: --no-latexdo not run LaTeX on

Re: [R] distance in the function kmeans

2004-05-28 Thread Uwe Ligges
n.bouget wrote: Hi, I want to know which distance is using in the function kmeans and if we can change this distance. Indeed, in the function pam, we can put a distance matrix in parameter (by the line pam-pam(dist(matrixdata),k=7) ) but we can't do it in the function kmeans, we have to put the

Re: [R] orca binary?

2004-05-28 Thread Uwe Ligges
Laura Holt wrote: Dear R People: Is there a binary version for rorca, please? Given I recall correctly, I tried to build rorca for Windows one or two years ago without success. But maybe Tony (In this context I recall that I'm still owing him a pizza) knows of a binary? Uwe Ligges Thanks,

Re: [R] distance in the function kmeans

2004-05-28 Thread Christian Hennig
On Fri, 28 May 2004, Uwe Ligges wrote: n.bouget wrote: Hi, I want to know which distance is using in the function kmeans and if we can change this distance. Indeed, in the function pam, we can put a distance matrix in parameter (by the line pam-pam(dist(matrixdata),k=7) ) but we

Re: [R] gauss.hermite?

2004-05-28 Thread Tamas Papp
On Fri, May 28, 2004 at 09:04:53AM +0200, Dimitris Rizopoulos wrote: Dear Spencer, In Dr. Gray's course notes for Advanced Statistical Computing (Appendix of Chapter 6) there are some functions for computing abscissas and corresponding weights for several Gaussian integration rules:

Re: [R] distance in the function kmeans

2004-05-28 Thread [EMAIL PROTECTED]
n.bouget wrote: Hi, I want to know which distance is using in the function kmeans and if we can change this distance. Indeed, in the function pam, we can put a distance matrix in parameter (by the line pam-pam(dist(matrixdata),k=7) ) but we can't do it in the function kmeans, we

Re: [R] Sorting Data?

2004-05-28 Thread Martin Klaffenboeck
Am 27.05.2004 22:11:04 schrieb(en) Berton Gunter: Your question does not make a lot of sense to me. If you only want to compute distances between mothers and THEIR daughters, you cannot have a matrix, because different mothers have different numbers of daughters. Thats true, but we took

[R] R CMD check, latex

2004-05-28 Thread TEMPL Matthias
Hello, I'm trying to build a package and the R CMD build works fine. The R CMD check produces an error, because no latex is found. But I have MikTeX installed and it works fine (only?) in combination with WinEdt. After reinstallation of MikTeX (no errors) I have the same problem. I´m not sure

Re: [R] Sorting Data?

2004-05-28 Thread Martin Klaffenboeck
Am 27.05.2004 22:17:04 schrieb(en) Patrick Burns: First comment: F is a bad name to use for an object as it masks the commonly used meaning of FALSE. Oh, thanks, I just took F because 'Fragebogen' is such a long variable. But you are right, I will change that. You'll probably get some better

RE: [R] R CMD check, latex

2004-05-28 Thread TEMPL Matthias
It was an (for me not understandable) problem on miktex I´m solving my problem, by copying ..\texmf\miktex\bin\etex.exe to latex.exe I´m wondering that now R CMD check works now. Sorry, for my question to R help. Matthias -Ursprüngliche Nachricht- Von: TEMPL Matthias Gesendet:

[R] Pr(|z|) in lme4

2004-05-28 Thread Dieter Menne
Dear List, I am struggling understanding S4 classes. For example, when GLMM summary(glmmML( whatever)) outputs the following line: Estimate Std. Error DF z value Pr(|z|) (Intercept)0.856 0.319 452.68 0.0073 How do I access the Pr column? Dieter

Re: [R] Statistics Avec R (Vincente Zoonekynd): portuguese translation

2004-05-28 Thread Ernesto Jardim
On Fri, 2004-05-28 at 00:36, José Cláudio Faria wrote: Dears friends, I have been translating and extended, with educational porposes only, the available material in the www (http://zoonek2.free.fr/UNIX/48_R/all.html) Statistcs avec R for the Portuguese of Brazil. I consider the

Re: [R] distance in the function kmeans

2004-05-28 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: n.bouget wrote: Hi, I want to know which distance is using in the function kmeans and if we can change this distance. Indeed, in the function pam, we can put a distance matrix in parameter (by the line pam-pam(dist(matrixdata),k=7) ) but we can't do it in the function

Re: [R] Pr(|z|) in lme4

2004-05-28 Thread Martin Maechler
Dieter == Dieter Menne [EMAIL PROTECTED] on Fri, 28 May 2004 11:04:13 +0200 writes: Dieter Dear List, I am struggling understanding S4 Dieter classes. For example, when GLMM Dieter summary(glmmML( whatever)) Dieter outputs the following line: Dieter

Re: [R] Pr(|z|) in lme4

2004-05-28 Thread Uwe Ligges
Martin Maechler wrote: Dieter == Dieter Menne [EMAIL PROTECTED] on Fri, 28 May 2004 11:04:13 +0200 writes: Dieter Dear List, I am struggling understanding S4 Dieter classes. For example, when GLMM Dieter summary(glmmML( whatever)) Dieter outputs the following line: Dieter

RE: [R] Pr(|z|) in lme4

2004-05-28 Thread Dieter Menne
Martin, s - summary(glmmML( whatever)) s ## to give the above str(s) ## to show the internal structure str() is my favorite, it's essentially how I learned R. It was no problem to find my way in nlme, but still lost in lme4. Uwe Liggges wrote: Well, I wonder what Dieter Menne

Re: [R] distance in the function kmeans

2004-05-28 Thread Martin Maechler
n\ == n\ bouget n on Fri, 28 May 2004 09:37:35 +0200 writes: n\ Hi, I want to know which distance is using in the n\ function kmeans and if we can change this distance. n\ Indeed, in the function pam, we can put a distance n\ matrix in parameter (by the line n\

Re: [R] distance in the function kmeans

2004-05-28 Thread Christian Hennig
On Fri, 28 May 2004, Martin Maechler wrote: n\ == n\ bouget n on Fri, 28 May 2004 09:37:35 +0200 writes: n\ Hi, I want to know which distance is using in the n\ function kmeans and if we can change this distance. n\ Indeed, in the function pam, we can put a distance

Re: [R] Sorting Data?

2004-05-28 Thread Martin Klaffenboeck
Am 27.05.2004 22:43:02 schrieb(en) Jonathan Baron: On 05/27/04 21:58, Martin Klaffenboeck wrote: Hello, Im reading through some manuals, but I cannot find my answer. I have a file containing many data: Vpn CodeFamily Age F1 F2 ... F17 1 1 M 46 1 2

Re: [R] distance in the function kmeans

2004-05-28 Thread Thomas Petzoldt
n.bouget wrote: Hi, I want to know which distance is using in the function kmeans and if we can change this distance. Indeed, in the function pam, we can put a distance matrix in parameter (by the line pam-pam(dist(matrixdata),k=7) ) but we can't do it in the function kmeans, we have to put the

Re: [R] distance in the function kmeans

2004-05-28 Thread [EMAIL PROTECTED]
I don't exactly understand what you do, could you show me the program that you execute to do that? n.bouget wrote: Hi, I want to know which distance is using in the function kmeans and if we can change this distance. Indeed, in the function pam, we can put a distance matrix in

[R] useR! 2004 keynote lecture slides

2004-05-28 Thread Achim Zeileis
Dear useRs, there has been some interest on this list into material about the presentations at useR! 2004. The slides from the keynote lectures are now available from the conference web page at http://www.ci.tuwien.ac.at/Conferences/useR-2004/ thanks to the R-core team members who were

Re: [R] distance in the function kmeans

2004-05-28 Thread Gavin Simpson
Thomas Petzoldt wrote: n.bouget wrote: Hi, I want to know which distance is using in the function kmeans and if we can change this distance. Indeed, in the function pam, we can put a distance matrix in parameter (by the line pam-pam(dist(matrixdata),k=7) ) but we can't do it in the function

Re: [R] distance in the function kmeans

2004-05-28 Thread Gavin Simpson
Gavin Simpson wrote: ... An R example is given in the help file of decostand() in Jari Oksanen's vegan library for two of the transformations mentioned above. ^^^ Pre-empting the usual response about proper terminology, I of course meant package not library. Gav --

Re: [R] - making a Windows library from Unix source code

2004-05-28 Thread Wolski
Hi! Normally you are building R packages using cd c:\Program Files\R\rw1081\src R CMD build fatigue Sincerely. *** REPLY SEPARATOR *** On 5/11/2004 at 6:32 PM Livingstone, Paul wrote: Hi All, I'm using R1.8.1 on Windows XP. I'm having trouble producing an R library from

Re: [R] privileged slots,

2004-05-28 Thread Duncan Murdoch
On Thu, 27 May 2004 22:55:20 +0200, Torsten Steuernagel [EMAIL PROTECTED] wrote : 2. Is it true that one can gain execution speed by accessing slots using the function slots(object,name,check=FALSE)-value would speed up execution? Try this: get(@-) function (object, name, value) {

[R] Adding key to simple plot() function

2004-05-28 Thread Laura Quinn
I am generating a simple x,y plot to show geographical positions where point data has been taken. Rather than plot a graph and lay points on top I've opted to generate the map from the point data, ie plot(-x,-y). This works fine, though I am wanting to label each point 1:20 and have a key

[R] Regression model type II

2004-05-28 Thread Benjamin PLANQUE
I am trying to fit regression models type II with R, but it seems to me that most (all) of the linear model functions are for type I regressions. Does anyone knows whether type II regressions functions exist in R. Benjamin -- Benjamin Planque IFREMER Laboratoire d'Ecologie Halieutique BP 21105

Re: [R] Adding key to simple plot() function

2004-05-28 Thread Sean Davis
Laura, For your axis issues, look at plotting without axes [plot(,axes=F) and then see ?axis for making them on your own. As for labeling your data, look at ?points and ?text. Sean On 5/28/04 8:27 AM, Laura Quinn [EMAIL PROTECTED] wrote: I am generating a simple x,y plot to show

[R] optim(method=SANN)

2004-05-28 Thread Robin Hankin
Hello List I'm working on a combinatoric problem in which the object is to minimize the badness() of a vector. I think this class of problem is only soluble by optim() using method=SANN. The badness() of anything is = 0, and when I've found a solution with zero badness, I want optim() to stop

Re: [R] distance in the function kmeans

2004-05-28 Thread Thomas Petzoldt
[EMAIL PROTECTED] wrote: I don't exactly understand what you do, could you show me the program that you execute to do that? I did such things sometimes ago, so the following is (as usual) without warranty. There are several methods, e.g. using Choleski factorization, singular value decomposition

[R] Simple list manipulation question

2004-05-28 Thread Sean Davis
I have a list of vectors $A AB BC CD $B GF HG FH FJ and I want to convert it into a dataframe of form A AB A BC A CD B GF B HG B FH B FJ Just can't quite come up with a nice R solution for it. Thanks, Sean __ [EMAIL PROTECTED] mailing list

Re: [R] Simple list manipulation question

2004-05-28 Thread Wolski
Hallo! x-unlist(yourlist) names(x) # maybee you have to postprocess the names a little dataframe(names(x),x) Sincerely Eryk *** REPLY SEPARATOR *** On 5/28/2004 at 8:49 AM Sean Davis wrote: I have a list of vectors $A AB BC CD $B GF HG FH FJ and I want to convert it into a

Re: [R] Simple list manipulation question

2004-05-28 Thread Gabor Grothendieck
R L - list(A=c(AB,BC,CD),B=c(GF,HG,FH,FJ)) R stack(L) values ind 1 AB A 2 BC A 3 CD A 4 GF B 5 HG B 6 FH B 7 FJ B Sean Davis sdavis2 at mail.nih.gov writes: : : I have a list of vectors : : $A : AB BC CD : : $B : GF HG FH FJ : : and I want to

[R] any simple examples for SJAVA please

2004-05-28 Thread Laura Holt
Hi! Are there any simple examples for SJAVA please? Thanks, Laura R Version 1.9.0 Windows _ Is your PC infected? Get a FREE online computer virus scan from McAfee® Security. http://clinic.mcafee.com/clinic/ibuy/campaign.asp?cid=3963

Re: [R] any simple examples for SJAVA please

2004-05-28 Thread Wolski
Hallo Laura! What I have heard at UseR there are no stable running SJava at the Windows platform. ( I wish that I am wrong!) I knew that Simon Urbanek is using and developing a java interface to R. google: simon urbanek Sincerely Eryk *** REPLY SEPARATOR *** On 5/28/2004 at

Re: [R] privileged slots

2004-05-28 Thread Wolski
Hi! The simplest way to find the document is ?checkSlotAssignmen A function which is described there and which Matthias Kohl mentioned in his e-mail You will find it also using the html search form of the html documentation searching for RClassUtils. Its true that the document is a little

Re: [R] any simple examples for SJAVA please

2004-05-28 Thread Vincent MUTEAUD
Hi, like you I'm trying to use SJava and I meet with many difficulties. I try to run examples from Calling R from Java written by Duncan Temple Lang available from http://www.omegahat.org/SJava (there are many examples on this web site). Now I have an error that my R version is R 1.6.1 while I

Re: [R] Sorting Data?

2004-05-28 Thread Jonathan Baron
On 05/28/04 12:21, Martin Klaffenboeck wrote: Am 27.05.2004 22:43:02 schrieb(en) Jonathan Baron: On 05/27/04 21:58, Martin Klaffenboeck wrote: Hello, Im reading through some manuals, but I cannot find my answer. I have a file containing many data: Vpn CodeFamily Age F1 F2

[R] Merging nlme output

2004-05-28 Thread Harold Doran
Dear list: I am trying to merge two files together from output I get based on the coef() command. Here is what I am running into. I have two simple linear mixed models mod1.lme-lme(math~year, data=sample, random=~year|childid/schoolid) mod2.lme-lme(math~year, data=sample,

Re: [R] Crossed random effects in lme

2004-05-28 Thread Douglas Bates
Søren Højsgaard [EMAIL PROTECTED] writes: Dear all, In the SASmixed package there is an example of an analysis of a split-plot experiment. The model is fm1Semi - lme( resistance ~ ET * position, data = Semiconductor, random = ~ 1 | Grp) where Grp in the Semiconductor dataset is

Re: [R] orca binary?

2004-05-28 Thread A.J. Rossini
Uwe Ligges [EMAIL PROTECTED] writes: Laura Holt wrote: Dear R People: Is there a binary version for rorca, please? Given I recall correctly, I tried to build rorca for Windows one or two years ago without success. But maybe Tony (In this context I recall that I'm still owing him a pizza)

RE: RE: [R] ANOVA and contrasts

2004-05-28 Thread Jos Cludio Faria
Dear FZ In one way I have not problem... The problem arise to establish the contrasts for 2 way with interactions and I desire the SS of each one of the contrasts. Thanks, Jos Cludio Faria UESC/DCET Brasil 73-634.2779 [EMAIL PROTECTED] [EMAIL PROTECTED]

Re: [R] axis.POSIXct: Datetime data and plotting

2004-05-28 Thread Don MacQueen
At 2:10 PM -0400 5/27/04, Barrowman, Nick wrote: But is there any easy way to avoid this problem? It is confusing and inconvenient that axis.POSIXct(1,d,format=f) and axis.POSIXct(3,at=d,format=f) give seemingly contradictory scales. In one case you provide the locations of the tick marks, in

[R] simulate dependent probabilities

2004-05-28 Thread Chris Wallace
I need to simulate from a random process and am not sure how to go about it. The process is the probability of an event occuring between a pair of points on a line. (This probability is between 0 and 0.5). I have estimates of these probabilities for a series of points, their standard errors and

[R] RE: SCO R

2004-05-28 Thread Siddique, Amer
SCO is grasping for straws. and is now gasping. vexatious lawsuits will only drive so far before puttering out. eventually the curtain gets pulled back on glorified attempts at racketeering. cheers. Message: 76 Date: Thu, 27 May 2004 21:08:12 -0500 From: Deepayan Sarkar [EMAIL PROTECTED]

Re: [R] Build R-1.9.0 with static libraries ?

2004-05-28 Thread Mélanie PELEGRINI-ISSAC
I am still unable to build a binary for R (linux Red Hat 9 distribution) which would be self-consistent (without shared objects and use of dynamics librairies). I would be very grateful if someone could explain and/or help me... so please set MAIN_LD as it suggests (and that's why I quoted

[R] How to insert a bitmap in a grph device

2004-05-28 Thread Luis Rideau Cruz
If I have a plot and I want to insert a bitmap in the same devicehow to do it? Thank you Luis Ridao Cruz Fiskirannsóknarstovan Nóatún 1 P.O. Box 3051 FR-110 Tórshavn Faroe Islands Phone: +298 353900 Phone(direct): +298 353912 Mobile: +298 580800 Fax:

Re: [R] RE: SCO R

2004-05-28 Thread Berton Gunter
My $.02 (sorry, couldn't resist). This is but one skirmish in the whole vexatious battle on intellectual property rights, which extends from software to many other venues, one of the most interesting of which is patents on human genes. An example I heard on the radio this morning was eHarmony.com

[R] how to pass defined function with more than one arguments to apply?

2004-05-28 Thread anoly
Dear all: I meet a problem of apply function. I have a matrix called tb tb V1 V2 V3 1 0 3 1 2 1 4 0 3 0 3 0 4 0 4 0 5 0 3 1 6 1 4 1 7 1 1 0 8 1 3 0 9 0 1 1 10 0 3 1 I hope to get the number of row that match c(0,3,1) I do this way: length(apply(t(tb) = =

Re: [R] how to pass defined function with more than one arguments to apply?

2004-05-28 Thread Wolski
Hi! ?apply ...: optional arguments to 'FUN'. you can pass more arguments to your function. length(apply(tb,comare,c(0,3,1))) Sincerely Eryk. *** REPLY SEPARATOR *** On 5/28/2004 at 12:05 PM anoly wrote: Dear all: I meet a problem of apply function. I have a matrix called tb

Re: [R] How to insert a bitmap in a grph device

2004-05-28 Thread Dirk Eddelbuettel
On Fri, May 28, 2004 at 04:58:58PM +0100, Luis Rideau Cruz wrote: If I have a plot and I want to insert a bitmap in the same devicehow to do it? IIRC the pixmap package on CRAN helps with that. Dirk -- FEATURE: VW Beetle license plate in California

Re: [R] How to insert a bitmap in a grph device

2004-05-28 Thread Thomas Petzoldt
Luis Rideau Cruz wrote: If I have a plot and I want to insert a bitmap in the same devicehow to do it? Yesterday we have learned, how to load and plot jpeg images (see: Is it possible to read jpeg files into R?) Now you need a way to insert this into another figure. For this purpose the grid

Re: [R] How to insert a bitmap in a grph device

2004-05-28 Thread Thomas Petzoldt
Dirk Eddelbuettel wrote: IIRC the pixmap package on CRAN helps with that. Ah, I see, ?addlogo is much easier than grid in this case, as one sees in a (slightly modified) version of the help example: x - read.pnm(system.file(pictures/logo.ppm, package=pixmap)[1]) plot(sample(1:100)) for (i in

RE: [R] Generate a sequence of random integer values

2004-05-28 Thread Ted Harding
On 28-May-04 Brian Davenhall wrote: I'm trying to generate a sequence of random integer values. I've tried to combine the random (r) and the sequence (seq) functions but this approach does not work. For example, if I use the following command: a - seq(1:100) a [1] 1 2 3 4 5

[R] dotchart questions

2004-05-28 Thread Richard Valliant
I am trying to put 3 dotcharts side-by-side with minimal space between each. Each chart is for a different variable, but the vertical axes are the same. I want to have vertical axis labels on the lefthand chart but no vertical axis labels on the other two. Plus, I would like very little space

[R] Negative binomial glm and dispersion

2004-05-28 Thread Dan Kehler
Using R 1.8.1, and the negative binomial glm implemented in MASS, the default when using anova and a chi-square test is to divide the deviance by the estimated dispersion. Using my UNIX version of S-plus (v 3.4), and the same MASS functions, the deviances are *not* divided by the estimated

[R] Is there a way to represent the ... argument in signature?

2004-05-28 Thread Wolski
Hi! I guess that it cant work. but maybee I am wrong. I would like to define a function rbind. rbind function (..., deparse.level = 1) { which works only for objects of class Myclass Is it possible to use something like setMethod(rbind,signature(...=Myclass,deparse.level=numeric) This gives an

Re: [R] Is there a way to represent the ... argument in signature?

2004-05-28 Thread Gabor Grothendieck
Have a look at rbind.zoo in the zoo package for an example. Wolski wolski at molgen.mpg.de writes: : : Hi! : : I guess that it cant work. but maybee I am wrong. : I would like to define a function rbind. : rbind : function (..., deparse.level = 1) : { : which works only for objects of

Re: [R] orca binary?

2004-05-28 Thread A.J. Rossini
I cross compiled it. If someone wants to donate a pre-configured Windows computer to my research group, I could test it. Else, someone else will have to. I can't. It is http://www.analytics.washington.edu/~rossini/rorca/rorca.zip best, -tony [EMAIL PROTECTED] (A.J. Rossini) writes: Uwe

Re: [R] Sorting Data?

2004-05-28 Thread Martin Klaffenboeck
Am 28.05.2004 15:18:58 schrieb(en) Jonathan Baron: Also the other things. I have to test the hypothesis: Does a daugther answer the questions (semantic differential) more equal the own mother and more different to the mothers of the other daugthers. I don't think this is a trivial problem at

Re: [R] RE: SCO R

2004-05-28 Thread Spencer Graves
The intellectual property questions raised by the GNU General Public License (GPL) are under active discussion in PDF Solutions, where Sundar Dorai-Raj and I work. Sundar and I have recommended linking PDF commercial software to R. However, we have to be careful how we do this, because

[R] How could I find R.exp?

2004-05-28 Thread Rui
Hi all, I try to create an import library for using dll. Following the instruction in the readme.package, I used the command “lib /def:R.exp /out:Rdll.lib”, however, the error message showed “Cannot open file R.exp”. I downloaded the source code R-1.9.0.tgz, but I can not find R.exp. Thanks!

Re: [R] How could I find R.exp?

2004-05-28 Thread Duncan Murdoch
On Fri, 28 May 2004 11:52:58 -0600, Rui [EMAIL PROTECTED] wrote : Hi all, I try to create an import library for using dll. Following the instruction in the readme.package, I used the command “lib /def:R.exp /out:Rdll.lib”, however, the error message showed “Cannot open file R.exp”. I downloaded

Re: [R] gauss.hermite?

2004-05-28 Thread Spencer Graves
Thanks to Andy Liaw, Marc Schwartz, Dimitris Rizopoulos and Tamas Papp for quick replies. Before I sent the email, I found that it was written by Jim Lindsey, but I could not find further information on him, so I appreciate the links. Also, the Advanced Statistical Computing course

[R] Problem: creating shared objects using lapack and blas

2004-05-28 Thread Zhu Wang
I am having trouble with creating shared objects with Fortran files, which use numerical libraries Lapack and Blas. I have read the section in Writing R Extensions but could not find what I needed and I am not at the stage to create a library. What I did was $R CMD SHLIB --output=car file1.f

[R] intro statistics course in San Francisco bay area

2004-05-28 Thread John David Duncan
I'm in San Francisco looking for a basic statistics course... using R. I'll be using R for real work later on, so I'd prefer not to have to do the coursework with some other software. All leads are appreciated... Thanks, - JD __ [EMAIL PROTECTED]

[R] Converting data frame to array?

2004-05-28 Thread TAPO (Thomas Agersten Poulsen)
Dear List, Please bear with a poor newbee, who might be doing everything backwards (I was brought up in pure math). I want to make a simple multi-linear regression on a set of data. I did some expreiments, and if X is a 4 by 2 array and Y is a 4 by 1 array, I can do a linear

Re: [R] vector normal to a plane

2004-05-28 Thread Wolski
Hi! In which library is this function? I cant find it using help.search() on my R installation. Eryk *** REPLY SEPARATOR *** On 5/28/2004 at 12:22 PM [EMAIL PROTECTED] wrote: Hi All, (I have a degree in math, but I am too embarassed to ask my colleagues, so here goes:) I

Re: [R] (Correction) vector normal to a plane

2004-05-28 Thread wwsprague
if I do R cross(v1, v2), I get .94. Huh? Meant to say R crossprod(v1, v2) Sorry __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] vector normal to a plane

2004-05-28 Thread wwsprague
Hi All, (I have a degree in math, but I am too embarassed to ask my colleagues, so here goes:) I would like to get a vector normal (orthogonal) to a plane formed by two other vectors. In matlab I do this: v1 = [.4, .6, .8]; v2 = [.9, .7, .2]; nn = cross(v1,v2) (gives ~[-.48, .65, -.24] if I

Re: [R] vector normal to a plane

2004-05-28 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Hi All, (I have a degree in math, but I am too embarassed to ask my colleagues, so here goes:) I would like to get a vector normal (orthogonal) to a plane formed by two other vectors. In matlab I do this: v1 = [.4, .6, .8]; v2 = [.9, .7, .2]; nn = cross(v1,v2) (gives

Re: [R] vector normal to a plane

2004-05-28 Thread Uwe Ligges
Uwe Ligges wrote: [EMAIL PROTECTED] wrote: Hi All, (I have a degree in math, but I am too embarassed to ask my colleagues, so here goes:) I would like to get a vector normal (orthogonal) to a plane formed by two other vectors. In matlab I do this: v1 = [.4, .6, .8]; v2 = [.9, .7, .2]; nn =

Re: [R] vector normal to a plane

2004-05-28 Thread Thomas Lumley
W is looking for the vector cross product (a specifically three-dimensional object important in physics and engineering). The crossproduct() function provides a matrix product so that crossprod(x,y) is t(x)%*%y, something completely different. For the three-dimensional case you could define the

RE: [R] Converting data frame to array?

2004-05-28 Thread John Fox
Dear Thomas, In fact, the more common way to fit a linear regression in R is to use variables in a data frame (or list) along with a model formula specifying the model. All of this is explained in the Introduction to R manual that is distributed with R: see, in particular, Sec. 6.3 on data

Re: [R] vector normal to a plane

2004-05-28 Thread Gabor Grothendieck
Thomas Lumley tlumley at u.washington.edu writes: : : W is looking for the vector cross product (a specifically : three-dimensional object important in physics and engineering). The : crossproduct() function provides a matrix product so that crossprod(x,y) : is t(x)%*%y, something completely

[R] GLMM error in ..1?

2004-05-28 Thread Spencer Graves
I'm trying to use GLMM in library(lme4), R 1.9.0pat, updated just now. I get an error message I can't decipher: library(lme4) set.seed(1) n - 10 N - 1000 DF - data.frame(yield=rbinom(n, N, .99)/N, nest=1:n) fit - GLMM(yield~1, random=~1|nest, family=binomial, data=DF,

Re: [R] how to pass defined function with more than one arguments to apply?

2004-05-28 Thread Gabor Grothendieck
Not sure what you intend with regard to length but to get a logical vector indicating which rows equal a particular vector: f1 - function(tb, row) apply(tb,1,function(x)all(x==row)) # or without using apply: f2 - function(tb, row) colSums( abs( (t(tb) - row) ) ) == 0 # or in terms of a

[R] glmmPQL:

2004-05-28 Thread Spencer Graves
I'm getting a strange error from glmmPQL. Consider the following sample code: set.seed(8) N. - 1000 z - rnorm(N.) pr.good - exp(-1e-4*exp(2+2*z)) quantile(pr.good) DF. - data.frame(yield=rbinom(N., N., pr.good)/N., Offset=rep(-10, N.), nest=1:N.) fit -

Re: [R] GLMM error in ..1?

2004-05-28 Thread Deepayan Sarkar
Is this the new experimental lme4 (version 0.6-x) ? If so, this is due to an error in our use of method dispatch. It has been fixed in the development version, and there should be a new release in a few days. On Friday 28 May 2004 19:32, Spencer Graves wrote: I'm trying to use GLMM in

[R] Re: Problem: creating shared objects using lapack and blas

2004-05-28 Thread Zhu Wang
To solve my own problem, it seems better to follow instructions in Writing R Extensions to build a library. I am making progress. Thanks, Zhu Wang Statistical Science Department Southern Methodist University Dallas, TX 75275-0332 On Fri, 2004-05-28 at 13:11, Zhu Wang wrote: I am having

[R] distance in the function kmeans

2004-05-28 Thread Jari Oksanen
My thread broke as I write this at home and there were no new messages on this subject after I got home. I hope this still reaches interested parties. There are several methods that find centroids (means) from distance data. Centroid clustering methods do so, and so does classic scaling