Re: [R] creating a certain type of matrix

2006-02-07 Thread Petr Pikal
Hi as only you know perfectly which halves and other portions of your matrices contain zeroes and which contain random numbers you has to finalize the function yourself. Here are few ideas. n-20 mat-matrix(0,n,(n/2)/5+1) #matrix of zeroes dd-dim(mat) # actual dimensions mat[1:(dd[1]/2),1]-1

Re: [R] R loop

2006-02-07 Thread Petr Pikal
Hi you are correct as usually. It is about 10 versions 1.2.0 - 2.2.1 and as there are 2 versions per year it is about 4-5 years. But nevertheless it seems to me as if I used R forever (and hopefully it will continue) Anyway subsetting data frames **was** similar in the old versions as it is

[R] install package in Mac

2006-02-07 Thread Rhett Eckstein
Hi, all I have problem about installing package (built in Windows) in Mac. In the R for MacOS X FAQ, it mentioned that Mac required the prebuilt package named as XXX.tgz. However, the package built under Windows named as XXX.zip. So how should I do to let the package (built in Windows) to

Re: [R] R loop

2006-02-07 Thread Uwe Ligges
Petr Pikal wrote: Hi you are correct as usually. It is about 10 versions 1.2.0 - 2.2.1 and as there are 2 versions per year it is about 4-5 years. But nevertheless it seems to me as if I used R forever (and hopefully it will continue) Anyway subsetting data frames **was** similar in

Re: [R] install package in Mac

2006-02-07 Thread Uwe Ligges
Rhett Eckstein wrote: Hi, all I have problem about installing package (built in Windows) in Mac. In the R for MacOS X FAQ, it mentioned that Mac required the prebuilt package named as XXX.tgz. However, the package built under Windows named as XXX.zip. So how should I do to let the

Re: [R] cdf of multivariate normal

2006-02-07 Thread Uwe Ligges
Juan Pablo Lewinger wrote: I was wondering if anybody has written R code to compute the cdf of a multivariate (or at least a bivariate) normal distribution with given covariance structure. See ?pmvnorm in package mvtnorm. Uwe Ligges __

Re: [R] R loop

2006-02-07 Thread Peter Dalgaard
Uwe Ligges [EMAIL PROTECTED] writes: Petr Pikal wrote: Hi you are correct as usually. It is about 10 versions 1.2.0 - 2.2.1 and as there are 2 versions per year it is about 4-5 years. [That calculus only goes back to 2001, though, previous to that the release rate was higher (and

Re: [R] panel.levelplot() for 2D histograms

2006-02-07 Thread Vermeiren, Hans [VRCBE]
-Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Monday, February 06, 2006 9:52 PM To: Vermeiren, Hans [VRCBE] Cc: [EMAIL PROTECTED] Subject: Re: panel.levelplot() for 2D histograms On 2/6/06, Vermeiren, Hans [VRCBE] [EMAIL PROTECTED] wrote: Dear R-wizards,

Re: [R] Scatterplot with corner at (0, 0)

2006-02-07 Thread Jim Lemon
Thomas L Jones wrote: I have two vectors, x and y, and I want to create a scatterplot. I want the lower left corner to be at (0, 0). I can create a scatterplot just by calling plot (x, y), but I am quite throroughly stymied and frustrated by the problem of putting the lower left corner

Re: [R] prehistoric versions of R -- 1995!

2006-02-07 Thread Martin Maechler
UweL == Uwe Ligges [EMAIL PROTECTED] on Tue, 07 Feb 2006 08:56:10 +0100 writes: UweL Petr Pikal wrote: Hi On 6 Feb 2006 at 15:57, ivo welch wrote: Date sent: Mon, 6 Feb 2006 15:57:30 -0500 From:ivo welch [EMAIL PROTECTED] To:

[R] qqplot

2006-02-07 Thread Vincent Negre
Hello, I would like to use qqplot() to compare two experimental distributions. But I do not understand how qqplot() compute quantiles. In fact, quantile() do not return the same results. Thank you for your help. Vincent. __ R-help@stat.math.ethz.ch

Re: [R] qqplot

2006-02-07 Thread Duncan Murdoch
On 2/6/2006 9:59 AM, Vincent Negre wrote: Hello, I would like to use qqplot() to compare two experimental distributions. But I do not understand how qqplot() compute quantiles. In fact, quantile() do not return the same results. Take a look at the source code. qqplot replaces the longer

Re: [R] R-2.2.1-INSTALL Issue with TCLTK

2006-02-07 Thread Bhola, Gautam
I am sorry but I think I didn't do a reply all, when I replied to Prof. Ripley. I did try what he suggested. I have installed tcl-8.4.12 and tk-8.4.12 in another folder and given explicit path to those lib folders for the tclConfig.sh tkConfig.sh. I also tried using the /usr/lib

Re: [R] cdf of multivariate normal

2006-02-07 Thread Marco Geraci
Hi, you can check ?pmvnorm in package 'mvtnorm' Marco --- Juan Pablo Lewinger [EMAIL PROTECTED] wrote: I was wondering if anybody has written R code to compute the cdf of a multivariate (or at least a bivariate) normal distribution with given covariance structure.

Re: [R] qqplot

2006-02-07 Thread François Pinard
[Vincent Negre] [...] I do not understand how qqplot() compute quantiles. Just type ``qqplot`` (without the parentheses) at the R prompt, to see the source code. ``qqplot`` does not especially compute quantiles, which are rather obtained directly through sorting its arguments. -- François

[R] Using R to process spectroscopic data

2006-02-07 Thread Dirk De Becker
Dear R-users, I would like to process some spectroscopic data with R, and I was hoping some people might have some example code on how to do this. I would like to be able to do the following things: * Detect outlier-spectra - This can be done by using scoreplot from the pls package * Determine

[R] R- License

2006-02-07 Thread Antonio_Paredes
Hello. We are trying to install R on our network, and I wanted to ask if there is a user license agreement. I will be grateful if somebody can send me a link to it; if one exists. Thank you very much *** Antonio Paredes USDA- Center for Veterinary

[R] Application of R

2006-02-07 Thread Andy Wong
information. Thanks for your advice. result 20060207.pdf Description: Adobe PDF document __ 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- License

2006-02-07 Thread Doran, Harold
Type license() for this info -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Tuesday, February 07, 2006 8:50 AM To: r-help@stat.math.ethz.ch Subject: [R] R- License Hello. We are trying to install R on our network, and I wanted to

Re: [R] R- License

2006-02-07 Thread Rau, Roland
Hi, From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Subject: [R] R- License Hello. We are trying to install R on our network, and I wanted to ask if there is a user license agreement. I will be grateful if somebody can send me a link to it; if one

Re: [R] R- License

2006-02-07 Thread Roger Bivand
On Tue, 7 Feb 2006 [EMAIL PROTECTED] wrote: Hello. We are trying to install R on our network, and I wanted to ask if there is a user license agreement. I will be grateful if somebody can send me a link to it; if one exists. In addition to license() from an installed R, you can link through:

[R] Prediction method for lowess,loess,lokerns,lpepa,ksmooth

2006-02-07 Thread Amir Safari
Hi Every Body, I don't know why some regression functions have no related prediction function. For example lowess, loess, lokerns, lpridge, lpepa, and ksmooth. What could help? Is there any global or wrapper function so that can help? Regards, Amir Safari

Re: [R] Prediction method for lowess,loess,lokerns,lpepa,ksmooth

2006-02-07 Thread Gavin Simpson
On Tue, 2006-02-07 at 06:24 -0800, Amir Safari wrote: Hi Every Body, I don't know why some regression functions have no related prediction function. For example lowess, loess, lokerns, lpridge, lpepa, and ksmooth. What could help? Is there any global or wrapper function so that

Re: [R] Prediction method for lowess,loess,lokerns,lpepa,ksmooth

2006-02-07 Thread Frank E Harrell Jr
Amir Safari wrote: Hi Every Body, I don't know why some regression functions have no related prediction function. For example lowess, loess, lokerns, lpridge, lpepa, and ksmooth. What could help? Is there any global or wrapper function so that can help? Regards, Amir Safari

[R] Antigen forwarded attachment

2006-02-07 Thread Antigen_PEGASUS
The entire message [R] install package in Mac, originally sent to you by [EMAIL PROTECTED] ([EMAIL PROTECTED]), has been forwarded to you from the Antigen Quarantine area. This message may have been re-scanned by Antigen and handled according to the appropriate scan job's settings. Entire

[R] Antigen forwarded attachment

2006-02-07 Thread Antigen_PEGASUS
The entire message Re: [R] install package in Mac, originally sent to you by [EMAIL PROTECTED] ([EMAIL PROTECTED]), has been forwarded to you from the Antigen Quarantine area. This message may have been re-scanned by Antigen and handled according to the appropriate scan job's settings.

Re: [R] Prediction method for lowess,loess,lokerns,lpepa,ksmooth

2006-02-07 Thread Prof Brian Ripley
On Tue, 7 Feb 2006, Gavin Simpson wrote: On Tue, 2006-02-07 at 06:24 -0800, Amir Safari wrote: Hi Every Body, I don't know why some regression functions have no related prediction function. For example lowess, loess, lokerns, lpridge, lpepa, and ksmooth. What could help? Is there

[R] displaying Cyrillic in RGui under Windows

2006-02-07 Thread Roger Levy
I have a data frame with Cyrillic text that I would like to be able to view under RGui. Unfortunately I can't figure out from the manuals how to make this happen -- can someone point me in the right direction? Thanks, Roger Levy __

[R] post-hoc comparisons following glmm

2006-02-07 Thread Michaël Coeurdassier
Dear R community, I performed a generalized linear mixed model using glmmPQL (MASS library) to analyse my data i.e : y is the response with a poisson distribution, t and Trait are the independent variables which are continuous and categorical (3 categories C, M and F) respectively, ind is the

Re: [R] prehistoric versions of R -- 1995!

2006-02-07 Thread Paul Gilbert
Martin Maechler wrote: UweL == Uwe Ligges [EMAIL PROTECTED] on Tue, 07 Feb 2006 08:56:10 +0100 writes: UweL Petr Pikal wrote: Hi On 6 Feb 2006 at 15:57, ivo welch wrote: Date sent: Mon, 6 Feb 2006 15:57:30 -0500 From:

Re: [R] displaying Cyrillic in RGui under Windows

2006-02-07 Thread Prof Brian Ripley
On Tue, 7 Feb 2006, Roger Levy wrote: I have a data frame with Cyrillic text that I would like to be able to view under RGui. Unfortunately I can't figure out from the manuals how to make this happen -- can someone point me in the right direction? You can only do this if you run RGui in a

[R] (second round) creating a certain type of matrix

2006-02-07 Thread Taka Matzmoto
Hi R users Here is what I got with help from Petr Pikal (Thanks Petr Pikal). I modified Petr Pikal's code to a little to meet my purpose. I created a function to generate a matrix generate.matrix-function(n.variable) { mat-matrix(0,n.variable,(n.variable/2)/5+1) #matrix of zeroes dd-dim(mat) #

[R] Start problem after package update

2006-02-07 Thread CG Pettersson
Hello all! R2.2.1, W2k. When I used update.packages() today, R wanted to update VR, cluster and nlme packages. I did so using the Danish mirror. Everything worked as usual during update, but after closing R and trying to restart I got an error message saying that the .Rdata could not be

Re: [R] symbol() function in 3d

2006-02-07 Thread Gavin Simpson
On Tue, 2006-02-07 at 16:23 +, Laura Quinn wrote: Hello, I was wondering if there is anything within the 3d rgl library that is similar to the symbol() function? I am hoping to overlay a surface3d() plot with some circles of varying size/color. I tried using the points3d() function, but

[R] elements from 'sem' function

2006-02-07 Thread andre tavares correa dias
Hi, I would like to print elements from sem (structural equation modeling) function (e.g., model-reproduced covariance matrix (C); estimated asymptotic covariance matrix of parameter estimates (cov)). How can I do this? Thanks, André -- André Tavares Corrêa Dias Laboratório de Ecologia

[R] question about binary data file and data.frame

2006-02-07 Thread [EMAIL PROTECTED]
I have a binary file with data sequence in the order [age,weight][age,weight] I know the length of the data and I want to load it into a data.frame. of course a way to do this is to read age and weight seperately and then use cbin(age,weight) to combine them into a dataframe, but is

[R] lme question

2006-02-07 Thread Mahdi Osman
Hi list, I am fitting microarray data (intensity) model using the lme package in R environment. I have 5 fixed variables in the model. One of the fixed variables is genes. I am trying to get p-values for different genes. But I am getting only one p-value for all genes together. I can get a list

[R] [R-pkgs] new package QCA

2006-02-07 Thread Adrian DUSA
Dear list members, I am pleased to let you know that R met with QCA - Qualitative Comparative Analysis. This package has a few functions that implement the Quine-McCluskey algorithm, adapted to social sciences by Charles Ragin (as describes in his book from 1987 The Comparative Method).

Re: [R] Significance of degrees of freedom in nlme

2006-02-07 Thread Spencer Graves
Your question begins, Dear Dr. Bates, yet Bates does not appear as an addressee in the copy of your question that I recieved via r-help. I will offer a couple of suggestions in the hopes that they might help you. Since your question did not include enough detail for anyone to

Re: [R] question about binary data file and data.frame

2006-02-07 Thread Barry Rowlingson
[EMAIL PROTECTED] wrote: I have a binary file with data sequence in the order What do you mean by 'binary file'? [age,weight][age,weight] How are age and weight encoded in this 'binary file'? I know the length of the data and I want to load it into a data.frame. of course a way to

Re: [R] displaying Cyrillic in RGui under Windows

2006-02-07 Thread Roger Levy
Prof Brian Ripley wrote: On Tue, 7 Feb 2006, Roger Levy wrote: I have a data frame with Cyrillic text that I would like to be able to view under RGui. Unfortunately I can't figure out from the manuals how to make this happen -- can someone point me in the right direction? You can only do

Re: [R] writing R shell scripts?

2006-02-07 Thread FWS
Use the all powerful here document feature of bash: R --vanilla EOF # Pipe all subsequent lines into R. ### Here's my R code require(tcltk) tkmessageBox(message=It works) ## the end EOF Mike Miller wrote: I'm new to the list. I've used R and S-PLUS a bit for about 15 years but

[R] March course @ 8 locations *** R/Splus Fundamentals and Programming Techniques

2006-02-07 Thread elvis
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce 2-day R/S-plus Fundamentals and Programming Techniques in San Francisco: www.xlsolutions-corp.com/Rfund.htm Philadelphia, March 9-10 Seattle, March 13-14 San Diego,

[R] Shapley Values

2006-02-07 Thread statistical . model
Hi, I am trying to compute the Shapley Values between a set of p independent variables (x1, ..., xp) and the dependent variable y to study the relationships. I believe that the package kappalab could be an appropriate choice. Is anybody able to give me some hints about the code I should write?

[R] hi frinds

2006-02-07 Thread aleid2001
Hello, i hope that every body are OK there. I would like to calculate the inverse of the cdf of mixture of normal distribution (with 2 components, p,(1-p)) from uniform points (0,1] using R. In this case i have to specify the 2 means and 2 variances and the mixing proportions (p, (1-p)). I am

Re: [R] elements from 'sem' function

2006-02-07 Thread John Fox
Dear André, The components of the sem object returned by sem() are explained on the help page ?sem. There are several methods for sem objects: methods(class=sem) [1] mod.indices.sem*normalized.residuals.sem* [3] path.diagram.sem* print.sem* [5]

[R] Reading in FORTRAN data using R

2006-02-07 Thread Tyler Hayes
Hi There: I was wondering if there is a way to read FORTRAN list data (similar to IDL's readf function). I often use FORTRAN for most of my number crunching, and use something like IDL to visualize and perform statistical analysis on that data. Since the each file is rather large (100 Mb),

[R] matching tables

2006-02-07 Thread Eric Pante
Dear Listers, I am trying to match tables that DO NOT have the same length. The tables result from the function table() so they look like this: table 1 2 3 4 3 5 7 table 2 1 2 3 6 4 5 I need the following output: (NOTICE THE ZEROS) 1 2 3 4 table1 0 3 5 7 table2 6 4 5 0

Re: [R] matching tables

2006-02-07 Thread bogdan romocea
t1 - as.data.frame(table(1:10)) ; colnames(t1)[2] - A t2 - as.data.frame(table(5:20)) ; colnames(t2)[2] - B t3 - merge(t1,t2,all=TRUE) -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Eric Pante Sent: Tuesday, February 07, 2006 4:22 PM To:

Re: [R] prehistoric versions of R -- 1995!

2006-02-07 Thread Uwe Ligges
Paul Gilbert wrote: Martin Maechler wrote: UweL == Uwe Ligges [EMAIL PROTECTED] on Tue, 07 Feb 2006 08:56:10 +0100 writes: UweL Petr Pikal wrote: Hi On 6 Feb 2006 at 15:57, ivo welch wrote: Date sent: Mon, 6 Feb 2006 15:57:30 -0500

Re: [R] question about binary data file and data.frame

2006-02-07 Thread [EMAIL PROTECTED]
Each data (edge,weight) is numerical data with type double (4 bytes) df-file(d:/sim.data,rb) age[1]-readBin(df,double()) weigt[1]-readBin(df,double()) ... age[10]-readBin(df,double()) weigt[10]-readBin(df,double()) it is not asc file. === 2006-02-07 13:03:17 === [EMAIL

[R] Calculating Kuiper/Watson test p-value

2006-02-07 Thread Firas Swidan
Hi, I am trying to use the Kuiper and Watson U^2 tests available in the CircStats or circular packages to test for uniformity on a circle. However, the implementation of these tests gives ranges of p-values that are not dependent on n, the size of sample. Are there any methods

Re: [R] Reading in FORTRAN data using R

2006-02-07 Thread Michael H. Prager
Tyler, One relatively easy way to get structured model data into R is to use output routines I have written with my colleagues. They are called For2R and are available at http://shrimp.ccfhrb.noaa.gov/~mprager/Rinter.html To go that route, you would replace (or supplement) your Fortran write

[R] Subsetting matrices

2006-02-07 Thread Tyler Smith
Dear R-Helpers, I am working on a dataset containing ca. 500 specimens of 17 different species. The data is three text columns, followed by ca 70 columns of binary data. I'm trying to write a function that will allow me to interactively exclude species using the inofrmation in the third text

Re: [R] image() and text

2006-02-07 Thread Paul Murrell
Hi Andrej Kastrin wrote: Dear useRs, I have 4×4 symmetrical matrix ; then I use image(log(my.matrix)) to visualise it. Is there any 'simple' way to add text labels into each cell lie on diagonal of the image plot? Thanks for any pointers... Do you mean something like ...? m -

[R] getting strata/cluster level values with survey package?

2006-02-07 Thread Jeff D. Hamann
First, I appoligise for the rooky question, but... I'm trying to obtain standard errors, confidence intervals, etc. from a sample design and have been trouble getting the results for anything other than the basic total or mean for the overall survey from the survey package. For example, using

[R] getting strata/cluster level values with survey package?

2006-02-07 Thread Jeff D. Hamann
First, I appoligise for the rookie question, but... I'm trying to obtain standard errors, confidence intervals, etc. from a sample design and have been trouble getting the results for anything other than the basic total or mean for the overall survey from the survey package. For example, using

[R] explanation of data sets variables

2006-02-07 Thread Kuba
Dear all, where I can find explanation of data sets variables, for example what does (u,x) variables means in city data? Thank you in advance, Kuba Jennifer Anistron, Mark Rufallo, Shirley Maclaine, Kevin Costner, Kathy Bates w komedii

[R] empty wireframe

2006-02-07 Thread Jeffrey Ross-Ibarra
Perhaps I am missing something incredibly basic, but I can't get the wireframe function to plot anything other than an empty graph with axes (and the correct scale). Cloud() works just fine, as does scatterplot3d(), and wireframe() works fine for other data sets I have. Are there some data sets

Re: [R] explanation of data sets variables

2006-02-07 Thread Francisco J. Zagmutt
library(boot) ?city Francisco From: Kuba [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Subject: [R] explanation of data sets variables Date: Wed, 08 Feb 2006 00:35:51 +0100 Dear all, where I can find explanation of data sets variables, for example what does (u,x) variables means in city

Re: [R] empty wireframe

2006-02-07 Thread Deepayan Sarkar
On 2/7/06, Jeffrey Ross-Ibarra [EMAIL PROTECTED] wrote: Perhaps I am missing something incredibly basic, but I can't get the wireframe function to plot anything other than an empty graph with axes (and the correct scale). Cloud() works just fine, as does scatterplot3d(), and wireframe() works

Re: [R] question about binary data file and data.frame

2006-02-07 Thread jim holtman
Read it into a vector and then convert into a matrix; convert to a data frame: df - readBin(d:/sim.data, what='double', n=1, size=4) # make n greater than file size dim(df) - c(2, length(df) / 2) df - t(df) colnames(df) - c('age', 'weight') df - as.data.frame(df) On 2/7/06, [EMAIL

Re: [R] how to use mle?

2006-02-07 Thread Spencer Graves
I see two problems: First, mle expects its first argument to return (-log(likelihood)); your LL is closer to +log(likelihood), but it seems to be missing something beyond that. In particular, Y is 10 x 3, but fit is only 10 x 2. Second, I'm not sure, but it

[R] plotting lines that break if data break

2006-02-07 Thread Denis Chabot
Hi, Sometimes data series (not necessarily time series) suffer breaks where data were expected, but not collected. Often the regular lines command to add such data to a plot is what I want, but other times I'd like the line to break where the data series is interrupted, instead of the

Re: [R] plotting lines that break if data break

2006-02-07 Thread Gabor Grothendieck
It does break already. Try this: plot(1:10, c(1:5,NA,7:10), type = l) plot(c(1:5,NA,7:10), 1:10, type = l) On 2/7/06, Denis Chabot [EMAIL PROTECTED] wrote: Hi, Sometimes data series (not necessarily time series) suffer breaks where data were expected, but not collected. Often the regular

[R] difference between rnorm(1000, 0, 1) and running rnorm(500, 0, 1) twice

2006-02-07 Thread Taka Matzmoto
Hi R users This looks a simple question Is there any difference between between rnorm(1000,0,1) and running rnorm(500,0,1) twice in terms of outcome ? TM __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] plotting lines that break if data break

2006-02-07 Thread Denis Chabot
I'm sorry, too long without coffee. Lately I have done numerous plots with superimposed non-linear quantile regression fits, and those do not break (which is fine), and I wrongly transposed this to lines in my head. Denis Le 06-02-07 à 22:06, Gabor Grothendieck a écrit : It does break

[R] lme syntax for PB examples

2006-02-07 Thread Paul Cossens
Hi helpeRs, I've been working through some examples in Pinhiero Bates( 2000) trying to understand how to translate to the new Lme4 syntax but without much luck. Below is what I think I should do, but either the answers don't come out the same or I get errors. In the Oxide problems I'm

Re: [R] 15-min mean values

2006-02-07 Thread Augusto.Sanabria
Thanks a lot to Gabor for his help with the solution of this problem. The solution using zoo(aggregate) is easy to implement efficient. I have calculated the 15min mean values of a 1min wind speed file containing 2.9 million records x 16 columns (size 179 MB) in just 144 seconds (R-2.1.1