Re: [R] How to measure mode (central tendancy)

Julian Burgos wrote: mode()... of course! Wrong. See ?mode Get or set the type or storage mode of an object. Thomas P. Raymond Balise wrote: What is the name of the function to give me the mode (central tendancy) of a numeric variable that can be negative? [[alternative HTML version

Re: [R] How to measure mode (central tendancy)

Raymond Balise wrote: What is the name of the function to give me the mode (central tendancy) of a numeric variable that can be negative? Similar questions appear occasionally on this list. Please see the following help page for possible solutions:

Re: [R] Using odesolve to produce non-negative solutions

Dear Jeremy, a few notes about your model: The equations of your derivatives are designed in a way that can lead to negative state variables with certain parameter combinations. In order to avoid this, you are using if constructions which are intended to correct this. This method is however

Re: [R] question about formula for lm

Why not using: lm(X[[Ytext]]~Xvar,data=X) ThPe Pedro Mardones wrote: Dear all; Is there any way to make this to work?: .x-rnorm(50,10,3) .y-.x+rnorm(50,0,1) X-data.frame(.x,.y) colnames(X)-c(Xvar,Yvar) Ytext-Yvar lm(Ytext~Xvar,data=X) # doesn't run lm(Yvar~Xvar,data=X) #

Re: [R] JSS Volume 18 (www.jstatsoft.org)

Thomas Petzoldt and Thomas Kneib, guest editors -- Thomas Petzoldt Technische Universitaet Dresden Institut fuer Hydrobiologie[EMAIL PROTECTED] 01062 Dresden http://tu-dresden.de/hydrobiologie/ GERMANY __ R-help

[R] Reminder: JSS Special volume R in ecology

Petzoldt, guest editors The original CFP was posted on Fri, 6 Oct 2006: https://listserv.umd.edu/cgi-bin/wa?A2=ind0610aL=ecolog-lP=11582 -- Thomas Petzoldt Technische Universitaet Dresden Institut fuer Hydrobiologie[EMAIL PROTECTED] 01062 Dresden http://tu-dresden.de

Re: [R] Lattice / Trellis analog of axis(graphics) ?

Hi Derek, see ?xyplot and ?panel.axis Hint: RSiteSearch(panel.axis) will point you to examples. Thomas Derek Eder wrote: My question is so basic that I am (almost too) embarrassed to admit that I could not find an answer after an hour's worth of homework. What is the Trellis / Lattice

Re: [R] ML fit of gamma distribution to grouped data

Dear Augusto, thank you for your example. Your solutions are the in fact the usual methods, but these do not apply to my case because I have grouped data. The good news is, that the solution of Prof. Brian Ripley works perfectly -- of course :-) Thank you both for your help. Thomas P.

Re: [R] Vertical line in densityplot?

Hi, lattice graphics work by utilizing so called panel functions. Here is a working version of your example: library(lattice) x-rnorm(100) plot.new() densityplot(x, panel=function(x, ...){ panel.densityplot(x, ...) panel.abline(v=0) } ) For mor information, please look into

[R] ML fit of gamma distribution to grouped data

Hello, we have a set of biological cell-size data, which are only available as frequencies of discrete size classes, because of the high effort of manual microscopic measurements. The lengths are approximately gamma distributed, however the shape of the distribution is relatively variable

[R] [R-pkgs] package simecol uploaded to CRAN

of incompatibility. Comments are welcome, Thomas Petzoldt -- Abstract: = The simecol package is intended to give users (scientists and students) and interactive environment to implement, distribute, simulate and document

Re: [R] Problems Creating an R package

Hi Aarti src/ if you have source files in C or Fortran. R/for the R sources. If you have no C or Fortran files, you should delete the src/ directory. The syntax error in your R function is simply that the , is on a new line so R thinks that the line above is complete. Hope it helps

[R] assign / environment side effect on R 2.4.0

other stuff } L - list(test = function() 1 + 2) e1 - environment(L$test) solver(L) e2 - environment(L$test) print(e1) # environment: R_GlobalEnv print(e2) # environment: 0x01d0b088 -- Thomas PetzoldtTel. +49-351-463 3 4954 Technische Universitaet DresdenFax +49-351-463 3

Re: [R] assign / environment side effect on R 2.4.0

Sorry, the posted example had the side effect on all platforms (correctly: R 2.2.1/Windows, 2.3.1/Linux, 2.4.0/Windows), but in the following corrected example the behavior of 2.4.0 differs from the older versions. The only difference between the wrong and the new example is L[[test]] vs. L$test

[R] list of interdependent functions

Hello, I discussed the following problem on the great useR conference with several people and wonder if someone of you knows a more elegant (or more common ?) solution than the one below. The problem: I have several sets of interrelated functions which should be compared. The

Re: [R] list of interdependent functions

Martin Morgan wrote: Here's another way: makeSolver - function() { f1 - function(x, K) K - x f2 - function(x, r, K) r * x * f1(x, K) function() f1(3,4) + f2(1,2,3) } solverB - makeSolver() solverB() makeSolver (implicitly) creates an environment, installs f1 and f2 into it,

Re: [R] Sampling in R

Does the following what you want: x - 1:100 s - matrix(0, nrow=200, ncol=9) for (i in 1:200) { s[i, ] - sample(x, 9) } m - rowMeans(s) hist(m) The default behavior of sample is without replacement. Thomas Noam Friedman wrote: Hi all! I was wondering if you can help me with a

[R] generating multivariate autocorrelated time series

Hello expeRts, for an application in hydrology I need to generate multivariate (log)normally distributed time series with given auto- and cross-correlations. While this is simple for the univariate case (e.g. with conditional normal sampling) it seems to be not so trivial for multivariate

Re: [R] Running DOS command prompt from R 2.0.1 Windows?

Mr Toby Daniel Ahrens wrote: I am trying to use R (windows version 2.0.1) to manage runs of a program that is run from a DOS command prompt. Is R able to call a DOS prompt? I am hoping that there is something analogous to the spawn command in IDL, but I can't see to find any help in the R

Re: [R] remedial stats education

Doran, Harold wrote: There is a Springer publication All of Statistics: a concise course in statistical inference by Larry Wasserman that might be what you are looking for. The book also has an emphasis on R and his web site has code and data sets for analysis of the examples used throughout.

Re: [R] oma and sub-title

Thomas Steiner schrieb: I want to add an outer subtitle to my 2x3-plot, but it's distance to the lowest plots is very high: 6 lines or more?! If I choose oma=c(5,0,2,0), it lies out of the plotting region and disapprears. Obviously I make something wrong here. Help appreciated, Thomas

Re: [R] Win32 network drive install

[EMAIL PROTECTED] wrote: I want to install r on a windows network drive so that users who have their own Win2000 machines can run something like x:\bin\rterm.exe CMD BATCH x:\url\prog.r c:\out\prog.Rout I do not want to make n users install their own versions of R; I want to

Re: [R] shapefiles manipulations ??

justin bem wrote: cheers , I need help in shapefile manipulations. I have two shapefiles of my country Cameroon. The first contain 10 provinces and each province contains a certain number of administratives units. I dont have ESRI Arc view. I want to add a admistrative unit of the center

Re: [R] Interpolating / smoothing missing time series data

Francisco J. Zagmutt wrote: I don't have much experience in the subject but it seems that library(akima) should be useful for your problem. Try library(help=akima) to see a list of the functions available in the library. I hope this helps Francisco Yes, function aspline() of package

Re: [R] Using R map data to determine associated state for a coordinate?

the example of inside.owin() from the spatstat package. Hope that helps Thomas Petzoldt ## library(maptools) library(spatstat) ger - read.shape(germany.shp) plot(ger) pger - Map2poly(ger) sx- pger[[13]] lines(sx, type=l, col=red) # Saxony

Re: [R] using system()

Omar Lakkis schrieb: Using system() is theer a way to make the R interpreter not wait for the command to finish? system(cmd, wait=FALSE) see ?system in online help. Thomas P. __ R-help@stat.math.ethz.ch mailing list

Re: [R] model selection vs. H0 based testing

Hello, I wish to thank Douglas Bates very much for clarification and pointing me to the MCMC simulation method to get p values even for cases where Wald tests are inappropriate. One question however remains when publishing statistical results: does it help readers if we combine both, - AIC

Re: [R] r: chinese installation of r

. If your system runs on Windows, define a variable LANGUAGE in the systems settings (environment) and set it to EN. Hope it helps Thomas Petzoldt __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] model comparison and Wald-tests (e.g. in lmer)

Dear expeRts, there is obviously a general trend to use model comparisons, LRT and AIC instead of Wald-test-based significance, at least in the R community. I personally like this approach. And, when using LME's, it seems to be the preferred way (concluded from postings of Brian Ripley and

Re: [R] TeX distribution on Windows

reports using R figures, - with SWeave and - in the R package building process. You can get the web installer from: http://sourceforge.net/projects/miktex Miktex can be used with WinEDT, Emacs, Texniccenter and others as editor. HTH Thomas Petzoldt __ R

Re: [R] png scaling problem - solved :-)

Yet another Windows solution without winfig: 1. Create a postscript image in R 2. Open this image in Ghostscript 3. Select a reasonable resolution using Display Settings in ghostscript 4. Copy the image via clipboard into your favorite image viewer (e.g. IrfanView) 5. Save the image in the

Re: [R] R binaries

Nam-Ky Nguyen schrieb: Dear Rexperts, I intend to burn some R CDs to colleagues in Vietnam. I want to put all binary files for base as well as contributed packages (for both Windows and Linux). It is very time consuming if I download files by files. Is there a place a can buy a CD with all

Re: [R] odesolve/lsoda differences on Windows and Mac

On 27 Jul 2005, Peter Dalgaard wrote: One thought: Integrating across input pulses is a known source of turbulence in lsoda. You might have better luck integrating over intervals in which the input function is continuous. Tweaking the lsoda tolerances is another thing to try. Yes, that's

Re: [R] odesolve/lsoda differences on Windows and Mac

. Sorry that I can't do more in the moment Thomas Petzoldt Thomas Petzoldt, Institute of Hydrobiology, Dresden University of Technology [EMAIL PROTECTED] http://www.tu-dresden.de/fghhihb/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] Help with Mahalanobis

Hello, a proposed solution of Bill Venables is archieved on the S-News mailing list: http://www.biostat.wustl.edu/archives/html/s-news/2001-07/msg00035.html and if I remember it correctly (and if the variance matrix is estimated from the data), another similar way is simply to use the Euclidean

Re: [R] How to change the value of a class slot

Ross Boylan wrote: I believe your example assumes that foo is updating the outer a by cheating and directly modifying enclosing environments. (I figure it also needs to get the name of its actual argument to do this, which would also involve slightly dirty tricks.) This does seem to be the

Re: [R] How to change the value of a class slot

Ross Boylan wrote: The paper also notes that call-by-value vs call-by-reference, which is the root of the slot update problem, is really orthogonal to the FOOP/COOP distinction. It's easy to imagine FOOP with call-by-reference. R.oo uses references. Yes, it can be found on

Re: [R] How to change the value of a class slot

Ross Boylan wrote: I defined an S4 class with a slot i. Then I wrote a regular function that attempted to increment i. [... details deleted ...] What do I need to do to update slot values? Here are some possibly relevant code fragments setClass(CompletePathMaker,

Re: [R] Time series indexes

Fernando Saldanha schrieb: I tried to assign values to specific elements of a time series and got in trouble. The code below should be almost self-explanatory. I wanted to assign 0 to the first element of x, but instead I assigned zero to the second element of x, which is not what I wanted. Is

[R] lme: error message with random=~1

Hello, I have an unbalanced mixed model design with two fixed effects site (2 levels) and timeOfDay (4 levels) and two random effects day (3 consecutive days) and trap (6 unique traps, 3 per site). The dependent variable is the body length (BL) of insect larvae from 7 to 29 individuals per trap

Re: [R] Converting integers to chars i.e 1 to 01

Gregor GORJANC wrote: Hello! I am producing a set of images and I would like them to be sorted by names I give. I was able to produce my names and add integer to them. That is easy. But my problem lies in sort of file from this process: figure_10.png figure_11.png figure_12.png ... figure_1.png

Re: [R] lme: error message with random=~1

Dimitris Rizopoulos wrote: Hi Thomas, random=~1 works if your data frame is in groupedData format, check this: # Orthodont is in groupedData format fm1 - lme(distance~age+Sex, data=Orthodont, random=~1) # dat - as.data.frame(Orthodont) fm2.1 - lme(distance~age+Sex, data=dat, random=~1) `dat'

Re: [R] lme: error message with random=~1

Douglas Bates wrote: I'm not sure what model you want to fit here. To specify a random effect in lme you need both a grouping factor and a model matrix. The error message indicates that lme is unable to determine a grouping factor. It would be correct syntax if you added a single level

Re: [R] find parameters for a gamma distribution

Andrew Collier wrote: hello, i have just started exploring R as an alternative to matlab for data analysis. so +far everything is _very_ promising. i have a question though regarding parameter +estimation. i have some data which, from a histogram plot, appears to arise from +a gamma

Re: [R] find parameters for a gamma distribution

Andrew Collier wrote: hello, i have just started exploring R as an alternative to matlab for data analysis. so +far everything is _very_ promising. i have a question though regarding parameter +estimation. i have some data which, from a histogram plot, appears to arise from +a gamma

Re: [R] find parameters for a gamma distribution

The help page of ?dgamma says: The mean and variance are E(X) = a*s and Var(X) = a*s^2. So, to estimate the parameters in your example, try the following: d - rgamma(10, 20, scale = 2) var(d)/mean(d) [1] 1.992091 mean(d)^2/var(d) [1] 20.09559 ... you may use these as start values for

Re: [R] Detecting incomplete commands

Philippe Grosjean wrote: Hello, I need a similar behaviour as with the prompt: asking to complete incomplete R command with eval(parse(text = )) Is it a way to make the difference between an illegal and an incomplete R command in a string? For instance: parse(text=ls()) expression(ls()) This

Re: [R] bootstrap package

mismatches) and at runtime (multi-argument returns are deprecated...) and it *still* works. Thomas Petzoldt __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting

[R] lattice graphics empty in Sweave

Hello, I have problems using lattice graphics together with Sweave under Windows XP and R 2.0.0; 2.0.1 Patched; 2.1.0 devel respectively. Base graphics are o.k., but for lattice graphics empty files (.eps and .pdf) are generated. The only workaround I found was to write and read the respective

[R] one dimensional AKIMA in R?

Hello, I am looking for an interpolation method similar to the one-dimensional AKIMA interpolation as in Akima (1970). Is there already such an algorithm in R which I may have overlooked? Thank you in advance Thomas P. H. Akima, A new method of interpolation and smooth curve fitting based on

Re: [R] surface fitting

m p wrote: Hello, I am looking for a R routine to fit a function/surface to my data 3d data. I'd like to use the function in a model so something like splines is not applicable. The data are smooth and can provide a plot if helpful. You may try surf.ls or surf.gls from package spatial. Thomas P.

Re: [R] akima

Erin Hodgess wrote: There is a contributed package called Akima which has the akima function. R Version 2.0.1 Windows The akima is from October 2004, so it is up to date. Yes I know. But, unfortunately package akima contains a different algorithm meant for bivariate interpolation and smooth

[R] decision about random effects in lme

Hello, in an experimental field study my collegue made a design with samples on two manipulated sampling sites (site: control, treatment). Within each site she sampled 3 traps (trap) at day and night (light: light, dark) at 3 consecutive days (day). We applied lme models with abundance as

Re: [R] R-(wiki)-pedia?

Tony Plate wrote: At Thursday 11:29 AM 10/7/2004, Dan Bolser wrote: [snip] I just added some pages... I think it would be great if people could get motivated to contribute to something like this. Its one of those cases of just getting the ball rolling... Do you think you can dump the existing

[R] R 2.0.0 (Windows): slow startup over the network

, with no effect. Does anyone an have idea, what I can do next? Thank you in advance Thomas Petzoldt __ [EMAIL PROTECTED] 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] Beginners problem

Rolf Wester wrote: Hi, I'm new to R and have a problem with a little test program (see below). Why doesn't - in function rk4 The reason is lexical scoping (see FAQ), and I suggest to implement rk4 as function with a return value (see # !!!). BTW there are already a rk4 (and lsoda) functions in

Re: [R] Indexing dataframe

Jacques VESLOT wrote: I am sorry to ask such question, but I can't find a solution... I have a dataframe 'd2004' and I want to remove two columns: 'd2004$concentration' and 'd2004$stade. d2004$concentration - NULL d2004$stade - NULL Hope it helps! Thomas P.

[R] Sweave echoing comments (again)

Hello, I try to document some R scripts for my collegues and observed the problem, that Sweave strips comment lines away. As a small example I write in an Rtex file: \begin{Scode} ## a small example test() # line comment \end{Scode} ... the .tex file generated by Sweave only contains:

Re: [R] newsgroup on R

Bin Jiang wrote: HI, by newsgroup, I mean a kind of user forum like MATLABs http://newsreader.mathworks.com/WebX?14@@/comp.soft-sys.matlab where anyone can post a question, without receiving unnecessary emails daily Hi, gmane.org IS a gateway from and to a newsgroup system in the original (NNTP)

Re: [R] importing data in excel

Uwe Ligges wrote: Gerardo Prieto Blanco wrote: Hello, I need to care excel data to be used in R,..., how do I make it? Thank you and greetings, Gerardo Prieto Please read the R Data Import/Export manual! Uwe Ligges Hello Gerardo, I completely agree with Uwe, and the following simple example may

Re: [R] text() with text, variables and math HOWTO?

Johannes Graumann wrote: Hello, One more question from the 'abusing R for blotting - particularly anally' department: How can I in the expression below make the '%~~%' show up as the aprrox-sign I want it to be? Thanks for any hint, Your code does not work because %~~% is a character string and

Re: [R] rbind for similar data frames

Gladys Castillo Jordán wrote: Hi all: I have a problem when I try to concatenate two similar data frames with different number of rows using rbind. I did something like this: d-data.frame(a=1:10,b=2:11,c=3:12) e-data.frame(a=101:105,b=102:106,c=103:107) data=rbind(d,e) The resulting row

Re: [R] Modalwert

Sonja Dornieden wrote: Hai - kann mir jemand sagen, wie ich den Modalwert in R berechne?! IRgendwie finde ich den Befehl nicht greetz und herzlichen Dank Sonja Hi Sonja, this question seems to appear in regular intervals and you find something about it in the R-Help archives, e.g.:

Re: [R] Population simulation.

[EMAIL PROTECTED] wrote: Hello, can anyone tell me if R has any special function for simulating the structure of human populations? Something like the genetic algorithm? I need to simulate a sample of a population with a specific structure. Is there something on R that can help me? Please be

Re: [R] questions about principle component analysis (princomp)

Yunfeng Hu wrote: Hi, I am a new R user and am currently using princomp to conduct a PCA. I have read the help(princomp) and still do not quite understand everything in the help. Basically I want to get the covariance matrix, and eigenvector/eigenvalues (loadings()?) so that I can find the

Re: [R] Vector to complex scalar

Laura Quinn wrote: Is there an R function to convert vectors into complex scalars? Thanks Is this because of the wind problem in one of the last questions? For conversion of two dimensional data into complex numbers and vice-versa you may look for ?complex in the online help, however I do not know

Re: [R] Histogram without common borders

Ross Darnell wrote: Is it possible to produce a histogram directly using the hist() function with the common borders removed? It can be done by plotting the histogram object using type 's'teps. my.hist - hist(x,plot=FALSE) plot(my.hist$breaks,c(0,my.hist$counts),type='s') Hello Ross, I think, your

[R] strucchange: breakpoints in inequally spaced data

Hello, we want to identify breakpoints (different phases) in environmental data, algae cell counts of three years with intervals between 7 and 30 days (N=40). We found that breakpoints(cells ~1) works great and identifies 5 very good breaks, however we are uncertain about these, because the

Re: [R] distance in the function kmeans

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

[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

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

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

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] adress the current index

Mag. Ferri Leberl wrote: How can I adress the current index of a vector? I want to work with time series and therefore give the n-th element of a vector some value dependent on the value of the n-1th element. Sorry, your question is not really clear, but possibly the following may help: 1) If

Re: [R] List of data frames...

[EMAIL PROTECTED] wrote: Hi, I would like to create a list of data frames that I could access via index manipulation. An array pointer of dataframe... for (i in 1:length(InputFilelist)) { # create data.frame temp - read.table (file = InputFilelist[i] , header = T, skip = 4) #

Re: [R] Colouring hclust() trees

Richard A. O'Keefe wrote: I have a data set with 6 variables and 251 cases. The people who supplied me with this data set believe that it falls naturally into three groups, and have given me a rule for determining group number from these 6 variables. One possibility is to extract the coordinates

Re: [R] Modalwert

Sonja Dornieden wrote: Hai - kann mir jemand sagen, wie ich den Modalwert in R berechne?! IRgendwie finde ich den Befehl nicht greetz und herzlichen Dank Sonja Hi, there was already a thread in this list about this question with subject Computing the mode on 24.02.2004. You will find several

Re: [R] Stichprobe bilden

Thomas Lumley wrote: On Thu, 29 Apr 2004, Thomas Petzoldt wrote: If you want to make a subsample, you can select cases (rows) using sub(), see ?sub for details, e.g. somethink like: I think you mean subset(). sub() is a regular expression substitution function. Absolutely sure! I should not try

Re: [R] Stichprobe bilden

[EMAIL PROTECTED] wrote: Liebe Alle, ich habe ien Problemm. ich habe einen Panel-Datensatz ueber die 185 Staaten (nach Code von IMF) und von 1948-1999. Ich habe auch 12 Variablen in diesem Datensatz, die ich fuer die Regression benoetige. ich muss ein sub-sample bilden fuer die 65 Staaten und

Re: [R] tapply

[EMAIL PROTECTED] wrote: Question: is there a function that average in a single step over the 3 columns? You may look for ?aggregate Thomas P. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] imputation of sub-threshold values [off-topic]

Thomas Lumley wrote: [...] If you can persuade the people measuring the values to give you the numbers (assuming they are just below `limit of detection' rather than genuine non-detects) you will reduce the need for imputation. This is often the most powerful technique -- analytical chemists

Re: [R] creating a ps. file

[EMAIL PROTECTED] wrote: Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I overwrite my results with the

Re: [R] creating a ps. file

[EMAIL PROTECTED] wrote: Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I overwrite my results with the

Re: [R] Simple numeric as.is question

Janet Gannon wrote: I am reading a list of numbers from my clipboard, and have been successful, except I can't make a histogram as R doesn't recognize my variable as numeric. I know I need to use as.is, but the specifics escape me. I have used x-read.table(clipboard, header=F) to import from

Re: [R] still spss

Margarida Júlia Rodrigues Igreja wrote: library(foreign) read.spss(H:\\Desktop\\bd1\\experiencia1) Error in read.spss(H:\\Desktop\\bd1\\experiencia1) : unable to open file I suspect, you make still a path error. Does your file really have no extension (e.g. .sav)? The default behaviour of

Re: [R] Summary: do.call and environments

Hello, f(fx,2) [1] 6 I would have naively expected 14. From whence cometh 6? Also, I prefer to use transportable code wherever feasible. The 2*3=6, which was the intention. It is in fact only a proof of correctness, that 7 is not used here. The proposal of Gabor does exactly,

Re: [R] Changingvalues in data frame

Hello Frank, does the following example, what you want? Thomas P. ## some test data x - data.frame(matrix(rnorm(72*10, mean=50, sd=20), ncol=10)) me - colMeans(x) sd - apply(x, 2, sd) coff - me + 2*sd # see ?t and ?pmax x2 - t(pmin(t(x),coff)) # test it x-x2

Re: [R] Changingvalues in data frame

Hello Frank, does the following example, what you want? Thomas P. ## some test data x - data.frame(matrix(rnorm(72*10, mean=50, sd=20), ncol=10)) me - colMeans(x) sd - apply(x, 2, sd) coff - me + 2*sd # see ?t and ?pmax x2 - t(pmin(t(x),coff)) # test it x-x2

Re: [R] julian.default on Windows (was years from as.POSIXlt)

Prof Brian Ripley wrote: Ok, I set my machine to CET and I get the problem. It is indeed a Microsoft bug: it seems its mktime can only regards as valid times after 1970-01-01 00:00:00 GMT (and not after that time in the current timezone). Hence R's strptime has trouble in the period of 1970

[R] do.call and environments

Hello, I want to call a function fx given by name, where some global variables (in the environment of fx) are passed to the function. For compatibility reasons I cannot modify the parameter list of fx and I want to avoid setting variables in the global environment (e.g. via -) Is there a way,

Re: [R] Center labels on a boxplot

Hello Andy, does mtext(c(SPP1, SPP2, SPP3), side=1, at=c(1,2,3)) complete your script? Thomas P. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] do.call and env

Erin Hodgess wrote: Hi Thomas! Try something like this: fx - function(y) print(x*y) f - function(fun,x) { assign(x,x,env=.GlobalEnv) fxx - function() { do.call(fun,list(y=3)) } fxx() } f(fx,13) [1] 39 Hope this helps! Thank you it works, but sorry, this is not what I want, as it assigns x

[R] Summary: do.call and environments

Dear R users, thank you very much for your suggestions. I've learned much, especially that the problem was not as simple as I thought. There have been several proposals, but most of them solved the problem only partly. The proposal(s) of Gabor and Tony (different versions) seemed to be very

Re: [R] years from as.POSIXlt

Prof Brian Ripley wrote: See ?julian, which says Note: Other components such as the day of the month or the year are very easy to computes: just use 'as.POSIXlt' and extract the relevant component. Hello, unfortunately not all mentioned functions work on all machines. Where

Re: [R] years from as.POSIXlt

as.numeric(format(x, f=%j)) which is the right code, works perfectly, too. Thomas P. __ [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

Re: [R] Sweave and Xemacs on Windows2000?

wolski wrote: Hallo! Trying to configure Xemacs to work with .snw files on windows 2000. Tried to do it how it is described in the FAQ for Sweaves. When starting xemacs with and Snw file *ESS* buffer contains hundrets of lines and the few last ones. [...] And no syntax highlighting in the

Re: [R] Computing the mode

). Thomas P. -- Thomas Petzoldt Dresden University of Technology Institute of Hydrobiology [EMAIL PROTECTED] 01062 Dresden http://www.tu-dresden.de/fghhihb/ __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman

Re: [R] r: plots

, line=2.5) # more about this see ?axis, ?pretty and ?mtext Thomas P. -- Thomas Petzoldt Dresden University of Technology Institute of Hydrobiology [EMAIL PROTECTED] 01062 Dresden http://www.tu-dresden.de/fghhihb/ __ [EMAIL

[R] Re: (read many files)

Claudia Paladini wrote: Dear ladies and gentlmen, I want to import a directory with about 400 files (.dat) in R. I know how to import a single file (with scan...) but I've good no idea how to import 400 at once. Can you help me ? Thanks a lot! Claudia setwd(c:/myfiles/) ## list.files uses

Re: [R] be careful: using attach in R functions

li dongfeng wrote: Hi there, I have just found that the ``attach'' function can get you into trouble when called many times. [..] Below is a demonstration of this performance loss, you will see a linear growth in CPU time usage. Adding a ``detach()'' call at the end of ``f'' will get

Re: [R] library nnet

Perez Martin, Agustin wrote: DeaR useRs: I am looking for a function which fits a multinomial model and in Baron´s page I find the function multinom in package nnet but this package is deprecated. Really? What do you mean with deprecated? I suppose that this function is now in other package but

  1   2   >