[R] 10 superscript x lables in log-plot

Hi ListMembers, I created a log plot with nums = c(1:10, seq(20,100, 10), seq(200,1000, 100), seq(2000,1, 1000), seq(2,10, 1),20) labl = rep(, length(nums)) labl[1] = 1; labl[10] = 10; labl[19] = 100; labl[28] = 1.000; labl[37] = 10.000; labl[46] = 100.000

Re: [R] error bars in lattice xyplot *with groups*

On 6/5/06, Mike Lawrence [EMAIL PROTECTED] wrote: Hi all, I'm trying to plot error bars in a lattice plot generated with xyplot. Deepayan Sarkar has provided a very useful solution for simple circumstances (https://stat.ethz.ch/pipermail/r-help/2005-October/081571.html), yet I am having

Re: [R] 10 superscript x lables in log-plot

maybe the first example of ?plotmath could be of help. 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:

[R] about string

Hello sir: There are 2 questions about string. 1 How to calculate the width of a string? e.g string abc's width is 3; 2 How can I get the substring in such kind of condition: f:\\JPCS_signal.txt f:\\PC1_signal.txt f:\\PC2_signal.txt What I wanna get is JPCS PC1 PC2.How can I achieve them by

Re: [R] about string

try the following: strg - abc nchar(strg) strg - c(f:\\JPCS_signal.txt, f:\\PC1_signal.txt, f:\\PC2_signal.txt) strg. - sapply(strsplit(strg, ), [, 2) sapply(strsplit(strg., _), [, 1) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student

Re: [R] time series clustering

Spencer Graves wrote: I know of no software for time series clustering in R. Google produced some interesting hits for time series clustering. If you find an algorithm you like, the author might have software. Alternatively, the algorithm might be a modification of something

Re: [R] about string

?nchar sapply(strsplit(f:\\JPCS_signal.txt, _), function(x) substring(x[1], 4)) --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33 (0)3.20.87.10.44 Fax : 33

Re: [R] How to call a value labels attribute?

Thank you, Richard. As soon as I find time I will carefully look at your solution and your book. Heinz At 10:01 05.06.2006 -0400, Richard M. Heiberger wrote: Aha! Thank you for the more detailed example. My solution for that situation is an attribute position and function as.position(). I use

[R] SPSS variable lables import

Hi, I try to get the variable labels of a SPSS data file into R but don't find this mentioned in the help file for foreign. Is there another way to get them ? BTW: An SPSS variable name is like: VAR001, whereas the variable label might be 'Identification no.' Thanks in advance, F. Thomas --

Re: [R] SPSS variable lables import

data1 - read.spss(file1.sav, F, T) # works ! --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33 (0)3.20.87.10.44 Fax : 33 (0)3.20.87.10.31 http://www-good.ibl.fr

Re: [R] about string

On Tue, 6 Jun 2006, Dimitris Rizopoulos wrote: try the following: strg - abc nchar(strg) strg - c(f:\\JPCS_signal.txt, f:\\PC1_signal.txt, f:\\PC2_signal.txt) strg. - sapply(strsplit(strg, ), [, 2) sapply(strsplit(strg., _), [, 1) or (a bit simpler)

[R] Problems using quadprog for solving quadratic programming problem

Hi, I'm using the package quadprog to solve the following quadratic programming problem. I want to minimize the function (b_1-b_2)^2+(b_3-b_4)^2 by the following constraints b_i, i=1,...,4: b_1+b_3=1 b_2+b_4=1 0.1=b_1=0.2 0.2=b_2=0.4 0.8=b_3=0.9 0.6=b_4=0.8 In my opinion the solution

Re: [R] Problems using quadprog for solving quadratic programming problem

G'day Fabian, FB == Fabian Barth [EMAIL PROTECTED] writes: FB I'm using the package quadprog to solve the following FB quadratic programming problem. FB I want to minimize the function FB (b_1-b_2)^2+(b_3-b_4)^2 FB by the following constraints b_i, i=1,...,4: FB

[R] Error in inherits(x, data.frame) : object Dataset not found

I have been trying to run a logistic regression using a number of studies. Below is the syntax, error message data. Any advice regarding what I am doing wrong or solutions are appreciated, regards Bob Green logreg - read.csv(c:\\logregtest.csv,header=T) attach(logreg) names(logreg)

Re: [R] Error in inherits(x, data.frame) : object Dataset not found

you probably want to use: model - glm(cbind(successes, failures) ~ medyear + age + sex + where + who + dxbroad + firstep + standard, family = binomial, data = logreg) since you store the data you imported in the data.frame 'logreg' not 'Dataset'. I hope it helps. Best, Dimitris

Re: [R] Error in inherits(x, data.frame) : object Dataset not found

On Tue, 6 Jun 2006, Bob Green wrote: I have been trying to run a logistic regression using a number of studies. Below is the syntax, error message data. Any advice regarding what I am doing wrong or solutions are appreciated, I suspect you wanted logreg - read.csv(c:\\logregtest.csv,

Re: [R] Multivariate skew-t cdf

Thanks to Spencer Graves for provinding this clarification about pmst. As the author of pmst, I was really the one expected to answer the query, but I was on travel in the last two weeks and did not read this query. As a complement to what already explained, the reason of the sharp change

Re: [R] Returned mail: see transcript for details

Dear Colleague, SAWN 2005 submission is now closed. For further details on the workshop program please see http://dna.engr.uconn.edu/SAWN2005/ SAWN 2005 Organizers __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] parameters in fitdistr

Hello, today I've updated from Version 2.1.1 to 2.3.1. The function you can find below for MLE with double tuncated Lognormal Distr. is working very well with 2.1.1 but now with the same function I get the error: object u_right not found ... although I've declared standard-values. I also tryed

[R] (no subject)

-- This e-mail and any files transmitted with it may contain privileged or confidential information. It is solely for use by the individual for whom it is intended, even if addressed incorrectly. If you received this e-mail in error, please

Re: [R] SPSS variable lables import

On Tue, 6 Jun 2006, Frank Thomas wrote: Hi, I try to get the variable labels of a SPSS data file into R but don't find this mentioned in the help file for foreign. Is there another way to get them ? BTW: An SPSS variable name is like: VAR001, whereas the variable label might be

[R] How to create list of objects?

Hi I am doing several mle and want to store them in a list (or whatever is the right construct) to be able to analyse them later. at the moment I am doing: f - list() f$IP - mle(...) f$NE - mle(...) but when I say: summary(f) I get: Length Class Mode IP 0 mle list NE 0

Re: [R] How to create list of objects?

lapply(f, summary) sapply(f, AIC) --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33 (0)3.20.87.10.44 Fax : 33 (0)3.20.87.10.31 http://www-good.ibl.fr

Re: [R] How to create list of objects?

try something like: lapply(f, summary) sapply(f, function(x) AIC(logLik(x))) 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:

[R] xYplot, lmline and abline

Dear All, I need to produce a multi-panel plot where: 1) groups within panels receive distinct symbols; 2) a linear regression is fit to the entire panel, not the individual groups; 3) a common abline is plotted in each panel. Essentially I would like to merge the two plots below, but I can't

[R] [R-pkgs] misc3d_0.4-0 now available on CRAN

Release 0.4-0 of package misc3d is now available from CRAN. This package provides a small collection of functions for 3D plotting, including 'contour3d' for computing and rendering 3D contours or isosurfaces and 'parametric3d' for rendering parameterized surfaces. A major change in this release

[R] use of R in big companies (references) R-support esp in Germany

Dear R users, sorry for this general email and I am sure it has been asked way too many times. IT departements in big companies only want to support the big standards. Whatever big standards means apart from being expensive. We are in the process of trying to get a risk management project

[R] Blocked File Type

Your message to r-help@stat.math.ethz.ch was blocked because it contained a file type prohibited by our server. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] [OFF] The best tool for a space-temporal analyses?

Hi, I try to make an analyses to discover what is the time that an area begin to have spacial autocorrelation. And after, what is the number of individuals responsible for this autocorrelation. The main idea is to discover if exist a contamination of a quadrat from others quadrats and how is

[R] memory.limit function not found

I have installed R 2.2.1 in Solaris 10 and am trying to increase the memory capacity (the system has 16G RAM) to 3 or 4G, but I keep getting: memory.limit(size=3000) Error: couldn't find function memory.limit Am I missing anything? I do that all the time under Windows. Any help would be

Re: [R] coloring leaves in a hclust or dendrogram plot [solved]

On Saturday 11 March 2006 07:37, Romain Francois wrote: Le 11.03.2006 15:56, Gabor Grothendieck a écrit : Where does one find A2Rplot that is called in that example? A2R does not appear to be on CRAN. I did find this: http://addictedtor.free.fr/Download/A2R.zip which is a zip

Re: [R] memory.limit function not found

On Tue, 6 Jun 2006, Priscila Darakjian wrote: I have installed R 2.2.1 in Solaris 10 and am trying to increase the memory capacity (the system has 16G RAM) to 3 or 4G, but I keep getting: memory.limit(size=3000) Error: couldn't find function memory.limit Am I missing anything? I do that

Re: [R] SPSS variable lables import

Thomas Lumley wrote: On Tue, 6 Jun 2006, Frank Thomas wrote: Hi, I try to get the variable labels of a SPSS data file into R but don't find this mentioned in the help file for foreign. Is there another way to get them ? BTW: An SPSS variable name is like: VAR001, whereas the variable label

Re: [R] memory.limit function not found

Priscila Darakjian wrote: I have installed R 2.2.1 in Solaris 10 and am trying to increase the memory capacity (the system has 16G RAM) to 3 or 4G, but I keep getting: memory.limit(size=3000) Error: couldn't find function memory.limit The function is only required on Windows. No need to

Re: [R] memory.limit function not found

Thanks, Prof Ripley. Now I understand. And , regarding the version, R was installed a while ago in this machine, when v. 2.3.1 was not available. Thanks again. Prof Brian Ripley [EMAIL PROTECTED] 6/6/2006 9:50 AM On Tue, 6 Jun 2006, Priscila Darakjian wrote: I have installed R 2.2.1 in

[R] Suspicious behaviour of sort on POSIXct vectors in R-2.3.0 and R-2.3.1

Hi , When I sort a vector of POSIXct values in R-2.3.0 and R-2.3.1, I get a vector of numeric values and this gets some of my code to crash (class object creation). Is that a R bug? In the mean time. I'll try to override the sort function for the POSIXct objects. Thanks, Patrick Guével Head of

[R] bootstrap data from groups

Hi R-friends. I have a mammal´s dataset looking like: Region Species Sex Bodysize 1 Sp1 M 10.2 1 Sp1 M 12.1 1 Sp1 M 9.1 ... I have three regions, four species and the body size

[R] Accessing lme source code

Dear all; This an FAQ. I tried to access lme source script so I can step into it to debug the problems resulting from a lme() call. I used getAnywhere(lme) or nlme:::lme, both produced only the function definition and UseMethod(lme). Any idea how to list the source code?

Re: [R] Accessing lme source code

1. You need to learn about S3 methods. ?UseMethod will tell you what you need to know. 2. methods(lme) will tell you the available methods. 3. nlme:::lme.formula will give you the code. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The business of the statistician

Re: [R] Suspicious behaviour of sort on POSIXct vectors in R-2.3.0 and R-2.3.1

Use order instead: x - as.POSIXct(Sys.Date() + 10:1) x[order(x)] On 6/6/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Hi , When I sort a vector of POSIXct values in R-2.3.0 and R-2.3.1, I get a vector of numeric values and this gets some of my code to crash (class object creation). Is

Re: [R] Fastest way to do HWE.exact test on 100K SNP data?

Anna Pluzhnikov [EMAIL PROTECTED] wrote: Hi everyone, I'm using the function 'HWE.exact' of 'genetics' package to compute p-values of the HWE test. My data set consists of ~600 subjects (cases and controls) typed at ~ 10K SNP markers; the test is applied separately to cases and controls.

[R] Score test to evalutate the proportional odds assumption.

To the list: What R commands will perform the score test on an ordered multinomial logit model to evaluate the proportional odds assumption? Many thanks. Mtichell Wachtel [[alternative HTML version deleted]] __

Re: [R] Suspicious behaviour of sort on POSIXct vectors in R-2.3.0 and R-2.3.1

On Tue, 6 Jun 2006, [EMAIL PROTECTED] wrote: Hi , When I sort a vector of POSIXct values in R-2.3.0 and R-2.3.1, I get a vector of numeric values and this gets some of my code to crash (class object creation). Is that a R bug? No, it is as documented: see ?sort As from R 2.3.0, all

[R] Strange behaviour of cbind

Hi, Is this intended behaviour of cbind? a-c(0,1,2,3) a [1] 0 1 2 3 a-as.ordered(a) a [1] 0 1 2 3 Levels: 0 1 2 3 a-a[a!=0] #remove the zero from a a [1] 1 2 3 Levels: 0 1 2 3 cbind(a) a [1,] 2 [2,] 3 [3,] 4 #cbind adds +1 to each element a-as.ordered(as.vector(a)) a [1] 1

[R] build R with Visual Studio

Hi, Has anyone had success in building R source with Visual Studio? I followed the instructions in README.packages, but failed on the very first step, where it's looking for R.dll. I looked through R source and couldn't find the file. Can someone point me to where this file is located or

[R] vague errors on R CMD check for very minimal S4-style package

Hello, I have a very minimal package simplepkg (DESCRIPTION, NAMESPACE, and R) with S4 classes/methods (defines a class foo and a show method for that class - both the class and show method are exported). I can seemingly install the package, then load and use it: sessionInfo() Version 2.3.1

Re: [R] build R with Visual Studio

Hi, R.ll is in my C:\Program Files\R\R-2.3.0\bin. My OS is Win XP. On 6/6/06, Jennifer Lai [EMAIL PROTECTED] wrote: Hi, Has anyone had success in building R source with Visual Studio? I followed the instructions in README.packages, but failed on the very first step, where it's looking

[R] Subset data in long format

I have data in a long format where each row is a student and each student occupies multiple rows with multiple observations. I need to subset these data based on a condition which I am having difficulty defining. The dataset I am working with is large, but here is a simple data structure to

Re: [R] Strange behaviour of cbind

On Tue, 2006-06-06 at 13:44 -0700, Dirk Vandekerckhove wrote: Hi, Is this intended behaviour of cbind? a-c(0,1,2,3) a [1] 0 1 2 3 a-as.ordered(a) a [1] 0 1 2 3 Levels: 0 1 2 3 a-a[a!=0] #remove the zero from a a [1] 1 2 3 Levels: 0 1 2 3 cbind(a) a [1,] 2

Re: [R] Subset data in long format

Apologies, but there were some word wrap issues in the prior email it seems. So, here is code for the sample data to avoid confusion tmp - data.frame(id = 1:3, matrix(rnorm(30), ncol=10) ) long - reshape(tmp, idvar='id', varying=list(names(tmp)[2:11]), v.names=('item'),timevar='position' ,

Re: [R] Strange behaviour of cbind

Hi, It doesn't have anything to with cbind, but rather with as.ordered - you converted it to a factor in that step. In the cbind step, you are actually getting the position of that ordered factor, rather than anything to do with the values themselves. a - c(1,7,5,3) a - as.ordered(a) a [1] 1

Re: [R] build R with Visual Studio

actually what I really want is to build the entire R source using Visual Studio. I know the current supported way is to use MinGW to build R on Windows. But since there is no 64-bit MinGW, I thought I'll experimenting with Visual Studio on 64-bit Windows. If anyone has success in building R

Re: [R] vague errors on R CMD check for very minimal S4-style package

Roels, Steven [EMAIL PROTECTED] writes: Hello, I have a very minimal package simplepkg (DESCRIPTION, NAMESPACE, and R) with S4 classes/methods (defines a class foo and a show method for that class - both the class and show method are exported). I can seemingly install the package, then

[R] [R-pkgs] zoo: new version 1.1-0

Dear useRs, I've just uploaded the new version 1.1-0 of the zoo package for dealing with regular and irregular time series data to the main CRAN site. Source and binary packages should be available from the mirrors in the next days. There are two changes that are nut fully backwards compatible:

Re: [R] Subset data in long format

Try this: subset(long, seq(id) - match(id,id) 6) On 6/6/06, Doran, Harold [EMAIL PROTECTED] wrote: I have data in a long format where each row is a student and each student occupies multiple rows with multiple observations. I need to subset these data based on a condition which I am having

[R] gamm error message

Hello, Why would I get an error message with the following code for gamm? I want to fit the a gam with different variances per stratum. library(mgcv) library(nlme) Y-rnorm(100) X-rnorm(100,sd=2) Z-rep(c(T,F),each=50) test-gamm(Y~s(X),weights=varIdent(form=~1|Z)) summary(test$lme) #ok

[R] spatial corStruct in lme

Hi, I'm fitting a relatively simple growth model to some forest plot data. Two species of trees were planted in different mixtures in 10 (nearly-adjacent) plots and measured on four occasions over 10 years. The model is constructed in terms of the diameter increments (per year; DI) in the 3

Re: [R] spatial corStruct in lme

Hello Lyndon, a couple of things you might try - for starters, the use of spatial correlation structures is documented in Pinheiro and Bates (2002) Mixed Effects Models in S and S-Plus from Springer-Verlag. So, you should try that. Secondly, it's difficult to provide useful advice in the

[R] knn - 10 fold cross validation

Hi, I was trying to get the optimal 'k' for the knn. To do this I was using the following function : knn.cvk - function(datmat, cl, k = 2:9) { datmatT - (datmat) cv.err - cl.pred - c() for (i in k) { newpre - as.vector(knn.cv(datmatT, cl, k = i)) cl.pred -

Re: [R] Nested and repeated effects together?

1. Am I correct that you are asking for a 4-factor interaction without a full set of 2- and 3-factor interactions? SAS may let you do that, but I'm not certain what that means. The standard treatment of a factor with k levels is to code it as (k-1) linearly independent contrasts

[R] S4 Methods with large numbers of args using memory

Thanks all to the people on the r-help list for the backtick help. I've run into another problem. With my SWIG C++ wrapper for R, I'm finding that S4 generic functions with large numbers of argument (such as 10 or 11) seem to be taking up a lot of memory especially when a lot of the arguments

Re: [R] knn - 10 fold cross validation

You might want to check out the function tune.knn() in the e1071 package. Andy _ From: [EMAIL PROTECTED] on behalf of Tim Smith Sent: Tue 6/6/2006 8:29 PM To: r-help@stat.math.ethz.ch Subject: [R] knn - 10 fold cross validation [Broadcast] Hi, I was trying to get the optimal

Re: [R] making matrix monotonous

I agree it would be great to sort the variables in a correlation matrix to make it easier to read and see patterns. I don't know any functions for doing that. If it were my problem, I might order the variables by their first principal component. There may also be some cluster

[R] x y averaging

Hello, I am trying to average a number of data sets where the x vector contains times and the Y data are instrument readings. The vectors all have different numbers of values, but many X time values match. For example: A fragment of the first data set: x vectvalue 14:56:10 0.325 14:62:11

Re: [R] x y averaging

Modifying the invalid times 14:62:11 and 14:68:11 to other values try the following: Lines1 - x y 14:56:10 0.325 14:57:11 0.111 14:58:11 0.214 16:05:18 0.245 16:11:09 0.266 16:17:05 0.271 16:33:00 0.304 16:39:05 0.300 Lines2 - x y 15:59:08 0.255 16:05:44 0.281 16:11:25 0.249 16:17:39

[R] Building packages in R - 'private' functions

Hello. I am creating an R package that I'd like to submit to CRAN (OS Windows XP). How do I distinguish among 'public' functions, e.g., those that are intended to be called by users of the package and for which I am providing documentation examples, and 'private' functions, which are used