[R] [R-pkgs] scuba 1.1-8

s usage calculations. New features in version 1.1-8: . support for dive profiles uploaded from a dive computer . new dataset: dive profile from a wreck dive on nitrox . bug fix in oxygen toxicity calculations Adrian Baddeley ___ R-packages

[R] [R-pkgs] QCA version 0.4-5

minimization function, using a shortcut instead of the classical complete and exhaustive algorithm. This new function is faster and uses significantly less memory (50 MB compared to 1.5 GB for large datasets). It should appear soon on CRAN, feedback is welcome. -- Adrian Dusa Romanian Social Data

[R] filehash package: error in unserialize(con)

r reading from connection > How do I get to touch the data again? It's really slow to recreate the filehash again, and it seems to be very fragile. Any suggestions much appreciated. Thank you, Adrian Dragulescu __ R-help@stat.math.eth

[R] C code generators

ld like to learn whether there's an R package or other open source program that already does it. Any pointers would be welcome thanks Adrian Baddeley __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] Problem Connecting to Oracle with R from Windows XP

You also need the ROracle package. On 8/15/07, Song, Alex <[EMAIL PROTECTED]> wrote: > > Hello, > > > > I installed RGui 2.5.1 and package DBI on Windows XP and tried to connect > to Oracle database which is on a Linux server. When I tried to use > dbDriver("Oracle"), I got an error as follows: >

[R] (no subject)

"Lawrence D. Brenninkmeyer" <[EMAIL PROTECTED]> writes: > I am trying to find a way to diffuse GIS data on a European map. I have a > dataset consisting of particular locations scattered across Europe, > along with magnitude and value information. The dataset is a spatial point pattern, where the

[R] vertically oriented color key in heatmaps

Rajarshi Guha <[EMAIL PROTECTED]> writes: > Hi, I have some data which I was plotting using image(). I wanted to > add a vertical color key to the plot and I found that heatmap.2 in > gplots does let me add a color key. However, I was thinking of a > vertical bar with the color range rather than

Re: [R] BDS test - results unclear to me

eries to be sure. Best regards Adrian > Subject: > [R] BDS test - results unclear to me > From: > Nicolas Navet <[EMAIL PROTECTED]> > Date: > Mon, 13 Aug 2007 22:14:40 +0200 > > To: > r-help@stat.math.ethz.ch > > > Hello, > > I would like to use t

Re: [R] bars' values on barplot

plot(my.values, ylim=c(0,11)) text(x, my.values, "wibble", pos=3) # always does what you want, whereas: text(x, 0.4+my.values, "wibble") # doesn't look very nice HTH, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5

Re: [R] the large dataset problem

$ operator. I think filehash is more general than R.huge. R.huge works very well with numerical 2D data only. Adrian Dragulescu On 7/31/07, Eric Doviak <[EMAIL PROTECTED]> wrote: > > > Just a note of thanks for all the help I have received. I haven't gotten a > chance to implement

Re: [R] Looping through all possible combinations of cases

sub.aa <- aa[rowSums(aa) == combinations, ] result <- apply(sub.aa, 1, function(x) sum(myvalues[x == 1])) names(result) <- apply(sub.aa, 1, function(x) paste(names(myvalues)[x == 1], collapse="")) HTH, Adrian On Friday 27 July 2007, Dimi

[R] LSD, HSD,...

Hi, I'm designing a experiment in order to compare the growing of several clones of a tree specie. It will be a complete randomized block design. How can I decide what model of mean comparision to choose? LSD, HSD,TukeyHSD, Duncan,...? Thanks in advance ___

Re: [R] Levene Test with R

the homogeneity of variances test, and the recommended ones are: Ansari-Bradley for two groups (i.e. for t.test) Fligner-Killeen for three or more groups (i.e. for ANOVA) hth, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel.

[R] [R-pkgs] spatstat version 2

SPATSTAT version 2 Spatstat is an R package for the statistical analysis of spatial data. A preliminary announcement about the forthcoming Version 2 of spatstat is available here: ___ R-packages mailing list [EMAIL PROTECTED] https:/

Re: [R] Odp: Odp: pie initial angle

tial angle, I just needed the position of the slices the way they are. My geometry seems to be poor towards innexistent :) All the best, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21

[R] pie initial angle

ion, where: pie(c(60, 40)) pie(c(80, 20)) present the slices in different directions. I read everything I could about init.angle argument, I even played with various formulas to compute it, but I just can't figure it out. How can I preserve the desired *direction* of the slices? Many thanks in

Re: [R] learning lattice graphics

Check the documentation link from http://cm.bell-labs.com/cm/ms/departments/sia/project/trellis/software.writing.html Adrian On 5/26/07, Tyler Smith <[EMAIL PROTECTED]> wrote: > > Hi, > > I've just produced my first lattice plot - the graphic is very > impressive, bu

[R] data in packages... a list?

Dear all, Is it possible to add a list in the data folder when creating a new package? In other words, is data in packages restricted to data.frame only? Thank you, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21

[R] how to avoid infinite loops when sourcing R files

"I add two numbers.\n"); flush.console() } Because of the source command in f1, I get into an infinite loop. This must be a common situation but I don't know how to avoid it. I need the source(set.parms) in f1, f2, f3, etc. because I want use a different combination of them in o

[R] R² in a non-linear regresion analisys

Can anybody explain me how do i get Correlation Coefficient R² in a non-linear regresion analisys performed with nls()?. Thanks in advance. -- Adrián J. Montero Calvo [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailm

[R] Quadratcount() plotting in R spatstat

d by plot(qX, add=TRUE, entries=NULL) The code in example(quadratcount) shows how to get thicker lines and change colours, etc. Adrian Baddeley __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

Re: [R] read.spss (package foreign) and SPSS 15.0 files

e multiple Excel fiels but it seems rather an overkill. In the mean time, I use portable files or export in different formats using a comercial software like StatTransfer (this one is really good). Thanks, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest s

Re: [R] read.spss (package foreign) and SPSS 15.0 files

don't have access to an SPSS 15 (only to version 14 for the moment) but I cannot find the option to save as CSV. Is it a version 15 feature? Thanks, Adrian __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do re

[R] sampling from the uniform distribution over a convex hull

points, library(spatstat) W <- convexhull.xy(x, y) P <- runifpoint(42, W) will compute the convex hull and generate 42 independent uniformly-distributed points in the convex hull. The result can be plotted by plot(P) and the coordinates of the simulated points can be extracte

Re: [R] Sorting rows of a binary matrix

func(3, 2) # For 3 columns in base 3 myfunc(3, 3) hth, Adrian On Thursday 22 February 2007 15:00, Adrian Dusa wrote: > Hello Serguei, > > Is this what you need? > > myfunc <- function(x) { > create <- function(idx) { > rep.int(c(rep.int

Re: [R] Sorting rows of a binary matrix

trix() in package QCA. HTH, Adrian On Thursday 22 February 2007 12:50, Serguei Kaniovski wrote: > Hallo, > > The command: > > x <- 3 > mat <- as.matrix(expand.grid(rep(list(0:1), x))) > > generates a matrix with 2^x columns containing the binary representations > of the

[R] R in a production environment

isibility and use of R beyond academia and personal use. Thanks, Adrian __ 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 and provide commented

Re: [R] jump in sequence

On Tuesday 30 January 2007 16:38, Peter Dalgaard wrote: > >[...snip...] > Is this it? > > > as.vector(outer(0:2,seq(4,22,9),"+")) > > [1] 4 5 6 13 14 15 22 23 24 Indeed it is :)) Thanks, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magu

[R] jump in sequence

, 3) Then see which are TRUE: which(myvec) [1] 4 5 6 13 14 15 22 23 24 I'd like to avoid creating the whole vector if possible; for very large ones it can be time consuming. There should be a way to only create the proper indexes... Thanks for any hint, Adrian -- Adrian Dusa Romanian Soc

[R] [R-pkgs] version 0.3 of QCA

ithmic one. This approach openes the way for _exact_ multi-value minimizations, and an even better (and faster) approach is searched for the future versions. Best, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40

Re: [R] multiple bases to decimal (was: comparing two matrices)

cbind(mat1, conv=colSums(apply(mat1, 1, function(x) x*mbase))) YES! Thank you so much Jim, this made my day :)) Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 3120210 / int.101

[R] multiple bases to decimal (was: comparing two matrices)

. Thanks for any hint, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 3120210 / int.101 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] comparing two matrices

quot;) > ind2 <- apply(mat2, 1, paste, collapse = "/") > match(ind2, ind1) Oh yes, I thought about that too. It works fast enough for small matrices, but I deal with very large ones. Using paste() on such matrices decreases the speed dramatically. Thanks again, Adrian -- Adrian

Re: [R] comparing two matrices

1[c(19, 16, 13, 24, 8), ] mylines <- mat2[, 1] for (i in 2:ncol(mat2)) {mylines <- 3*mylines + mat2[, i]} mylines + 1 [1] 19 16 13 24 8 I was still hoping for a direct matrix function to avoid the for() loop. Thanks, Adrian On Sunday 21 January 2007 01:06, Marc Schwartz wrote: >

Re: [R] comparing two matrices

Hi Christos, It's... more or less the same thing. I was looking for a matrix.to.matrix.run.me() function :) Cheers, Adrian On Sunday 21 January 2007 00:56, Christos Hatzis wrote: > Here is a slightly more compact version of your function which might run > faster (I did not test tim

[R] comparing two matrices

which(apply(mat1, 1, function(y) { sum(x == y) }) == ncol(mat1)) }) The code is vectorized, but I wonder if there is a simpler (hence faster) matrix computation that I miss. Thank you, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Buc

[R] [R-pkgs] spatstat 1.11-0

ll handled). More stringent checking of function arguments. Improved handling of pixellation effects. Extensions to variance-covariance calculations for fitted models. Adrian Baddeley and Rolf Turner ___ R-packages mailing list

Re: [R] a question of substitute

want more complicated things. I'm looking forward to the release of the 2.5.0 version. Best regards, Adrian On Thursday 11 January 2007 03:08, Prof Brian Ripley wrote: > The 'Right Thing' is for oneway.test() to allow a variable for the first > argument, and I have altered

Re: [R] a question of substitute

ithout > the workarounds I posted: >fun.2(values ~ group) Brilliant :) Super fast change, this is why I love R. Cheers, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21

Re: [R] a question of substitute

h fligner.test(), printing the p.value and based on that changing the var.equal argument in the oneway.test() It's just for convenience, but they do like having it all-in-one. Best regards, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5

Re: [R] a question of substitute

On Tuesday 09 January 2007 15:14, Gabor Grothendieck wrote: > oneway.test is using substitute on its arguments so its literally > getting formula rather than the value of formula. Ah-haa... I understand now. Thanks for the tips, they both work as expected. Best, Adrian > Try these:

[R] a question of substitute

o vectors in a data.frame, with no avail. I did find a hack, creating two new vectors inside the function and creating a fresh formula, so I presume this has something to do with environments. Could anybody give me a hint on this? Thank you, Adrian -- Adrian Dusa Romanian Social Data Archi

[R] dashed lines and SVG files

example is: library(RSvgDevice) devSVG("myplot.svg", width=10, height=10) plot(1:10) abline(v=5, lty=¨dashed¨) dev.off() Is there anything more (or different) I should do? Many thanks in advance, Adrian __ R-help@stat.math.ethz.ch mailing

[R] dashed lines and SVG files devSVG("/folderul/unde/salvez/myplot.svg", width=10, height=10) plot(1:10, 1:10) dev.off()

10) plot(1:10) abline(v=5, lty=¨dashed¨) dev.off() Is there anything more (or different) I should do? Many thanks in advance, Adrian __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide ht

[R] autofill for a tkentry field

try1.f <- tkentry(ppw, width=30, textvariable=f1, validate="key", validatecommand=autofill) tkgrid(entry1.f) tkgrid(ppw) tkgrid(pw) I don't know how to pass arguments to the autofill function, so I can manipulate the textvariable. Any

Re: [R] how to replace some objects?

t;b" with "2", and "c" with "3". > > Like this: 1,2,1,3,2,2,1,1,3. > > let <- c('a', 'b', 'a', 'c', 'b', 'b', 'a', 'a', 'c') > library(car) > num <- recod

Re: [R] Count cases by indicator

uot; argument: truthTable(df, outcome="OUT", inside=TRUE) I hope it helps, Adrian On Monday 04 December 2006 15:24, Liaw, Andy wrote: > I might be missing something, but the data you showed don't seem to > match your expectation. Firstly, 1 in binary is 511 in decimal

Re: [R] Diebold Mariano Test

Hello Graham Pls find attached some old R code for the Diebold Mariano test. The code is at least 6 years old and was not used in the meantime. Pls check first before using it. Best regards Adrian Dear List Has anyone used R to distnguish between alternative forecasting models? In

[R] Installing package rpvm under Windows

nfo' rpvm_ser.o:rpvm_ser.c:(.text+0x25): undefined reference to `pvm_pkbyte' rpvm_ser.o:rpvm_ser.c:(.text+0x64): undefined reference to `pvm_pkbyte' rpvm_ser.o:rpvm_ser.c:(.text+0x9f): undefined reference to `pvm_upkbyte' rpvm_ser.o:rpvm_ser.c:(.text+0xe4): undefined reference to

Re: [R] Suppress blanks/spaces in character

> c <- "abcdef ghi" > gsub(" +"," ", c) [1] "abc def ghi" > On Fri, 27 Oct 2006 [EMAIL PROTECTED] wrote: > Hi all > > I'm have a character vector and would like to suppress the blanks if there > are more than one after the other. > > Example: > > Character value is: "abc def ghi" > T

[R] ROracle error in Windows. Memory could not be read.

quot;0x4a280ae2". The memory could not be read. And R freezes. Something is working, because I get to establish an Oracle driver. Any ideas would be much appreciated. Thank you, Adrian __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.c

Re: [R] weight cases?

On Saturday 14 October 2006 16:52, Gabor Grothendieck wrote: > Try this (and round the result to make to it comparable to your > calculation): > > xtabs(weight ~ var1 + var2, my.data) Oh yes... :) It was so simple. Thanks for the cov.wt() as well. Regards, Adrian -- Adrian Dusa Rom

Re: [R] weight cases?

;- table(my.data$var1, my.data$var2)) round(unweighted*total$weight, 0) Yet another question: how would the weight variable be applied to correlate two numerical variables? Best, Adrian On Saturday 14 October 2006 16:00, Gabor Grothendieck wrote: > Try this: > > table(lapply(my.

[R] weight cases?

.data$var1, my.data$var2) A B C D E a 1 0 0 0 0 b 0 0 0 1 0 c 0 0 0 0 1 d 0 0 1 0 0 e 0 1 0 0 0 Applying the weight variable, the table should yield a value of 2 for the "eB" combination: > table(my.data$var1, my.data$var2) A B C D E a 1 0 0 0 0 b 0 0 0 1 0

Re: [R] How to repeat vectors ?

Maybe this one? > MyMatrix <- matrix(1:4, nrow=2) > MyMatrix [,1] [,2] [1,]13 [2,]24 > MyMatrix[rep(seq(nrow(MyMatrix)), each=2), ] [,1] [,2] [1,]13 [2,]13 [3,]24 [4,]2 4 HTH, Adrian On Saturday 30 September 2006 09:33, Ton

Re: [R] delete a entire vector of a dataframe

dat) %in% "Var_7") # now don't want Var_7 > dat <- dat[, -not.want] > dat > > This can be extended to many variables: > > not.want <- which(names(dat) %in% c("Var_10", "Var_2", "Var_8")) > dat <- dat[, -not.want] >

[R] [R-pkgs] The rpanel package

tcltk variables are managed behind the scenes so that users need not be concerned with any technicalities. A descriptive paper and some simple example scripts are available at www.stats.gla.ac.uk/~adrian/rpanel/ The package is now available from CRAN. Best wishes, Adrian Bowman (on behalf

[R] [R-pkgs] spatstat 1.9-5

h over 270 functions. For more information: See the Release Notes for version 1.9-5 at Adrian Baddeley and Rolf Turner ___ R-packages mailing list [EMAIL PROTECTED] https://stat.ethz.ch/mailman/listinfo/

[R] test the tcltk package

fine? I have R 2.3.1, tcl and tk version 8.4 (dev packages installed as well). Thanks in advance, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 31

Re: [R] R Site Search directly from Firefox's address bar

On Friday 18 August 2006 10:08, Romain Francois wrote: > Le 17.08.2006 20:56, Adrian Dusa a écrit : > [...] > > It breaks also every usage of the google feeling lucky default behaviour > which is really useful I think. > There are R related firefox search plugins on mycroft.

Re: [R] R Site Search directly from Firefox's address bar

"r". Now, everytime I type "r" in the address bar it takes me to www.r-project.org :) -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 3120210 / int.101 __

[R] R Site Search directly from Firefox's address bar

ions&idxname=docs&idxname=Rhelp02a&query="; >From now on, every keyword(s) you type in the address bar will take you directly to the first page of hits at http://finzi.psych.upenn.edu I found this very helpful. Best, Adrian -- Adrian Dusa Romanian Social Data Archive 1,

Re: [R] RE : tcltk library on linux

.R-project.org/bin/linux/ubuntu/ dapper/ Then: $ sudo apt-get update $ sudo apt-get install r-cran-rcmdr (as Dirk Eddelbuettel advised). It _should_ work flawlessly. HTH, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 312

Re: [R] objects and environments

is not > nearly so ugly so you might reconsider whether its ok > for you to just pass bb. Aah-aaa!! :) So that's the way to do it... I don't know how many times I read the help from apply and I missed it every time. Well, I learned many things today, I feel much better now

Re: [R] objects and environments

- apply(mymatrix, 1, fun2) > } Beautiful :) Thanks very much Dimitris, I was out of energy after several hours of struggling with this. Best, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Te

[R] objects and environments

4, e2) fun2() } fun2 <- function(idx) { get("bb", e2) } > fun1() Error in get("bb", e2) : object "e2" not found Any hint would be highly appreciated, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5

Re: [R] missing documentation entries

uestion for R-devel: please do study the posting guide. > > `R-devel is intended for questions and discussion about code development > in R.' Thank you very much for your reply, I'll post to R-devel from now on. It seems to me that name spaces are the solution for my problem. Adr

[R] missing documentation entries

d have two questions: - is it mandatory to document _all_ functions (will the source package be rejected by CRAN if otherwise)? - if not, is there a way to tell R which are the functions that I don't want to document? Thanks, Adrian -- Adrian Dusa Romanian Social Data Archive 1, Schitu Magu

Re: [R] Tcltk package

.org/bin/linux/ubuntu/ dapper/ This repository has lots of other packages compiled for Ubuntu, feel free to take a look. HTH, Adrian -- Adrian DUSA Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 312

[R] [R-pkgs] spatstat 1.9-4

Version 1.9-4 of package 'spatstat' has been sent to CRAN. It includes new code for perfect simulation of point processes and various improvements. The release notes are available at <http://www.spatstat.org/spatstat/current/spatstatRELEASE-NOTES-1.9-4>

Re: [R] Extracting from a matrix w/o for-loop

volving matrix indexing (M1[M2]) is simpler but just for the sake of it, since we're dealing with matrices it is not a case of sapply but of _apply_: apply(M2, 1, function(x) M1[x[1], x[2]]) My 2c, Adrian -- Adrian DUSA Arhiva Romana de Date Sociale Bd. Schitu Magureanu nr.1 050025 Bucur

[R] greek letters, text, and values in labels

for", theta, "=", eval(theta or using bquote plot(1:10, main=paste("Results for ", bquote(theta == .(theta or using substitute, etc. I could not make it work. This should be easy. I would appreciate your help. Thanks, Adrian

Re: [R] Getting rid of for loops

[2]]))[[1]]$"Pr(>F)"[3])     } p.interaction <- apply(all.pairs, 1, myfun) Actually, you don't need "as.matrix" there, just cbind all your vectors to obtain the final dataframe: finally <- as.data.frame(cbind(marker1, marker2, p.interaction)) Adrian On

Re: [R] Getting rid of for loops

]$"Pr(>F)"[3]) } p.interaction <- as.matrix(apply(all.pairs, 1, myfun) HTH, Adrian On Monday 17 July 2006 05:18, Kevin J Emerson wrote: > Hello R-users! > > I have a style question. I know that for loops are somewhat frowned upon > in R, and I was trying to figure out a

[R] [R-pkgs] new package QCAGUI

ment and advice. Regards, Adrian -- Adrian DUSA Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 3120210 / int.101 ___ R-packages mailing list [EMAIL PROTECTED] https://stat.

[R] Simulate Strauss process in 3D

stat' that will deal with 3D point patterns, including simulation of the Strauss process in 3D. regards Adrian Baddeley __ 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] checking package dependencies

for the answer. The problem seem to have vanished (and I haven't done anything in particular). When it didn't work, I remember I checked the permissions and everything was OK. Really have no idea what went wrong, but now it works. Best, Adrian -- Adrian DUSA Romanian Social Data Archiv

[R] uniform and clumped point plots

clumps. Hope this helps Adrian Baddeley __ 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

[R] checking package dependencies

s specify the path to the installed packages (and all depending packages are installed): R_LIBS=${R_LIBS-'/home/adi/Installed/R/site-library:/usr/local/lib/R/site-library:/usr/local/lib/R/library'} Is there something changed about defining R_LIBS? Thank you in advance, Adrian -- Adri

[R] bivariate weighted kernel density estimator

the moment it only offers the Gaussian kernel but we plan to include other kernels. Adrian Baddeley __ 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

[R] weighted kernel density estimate

Weighted spatial kernel density estimation is available in the function 'density.ppp' in the package 'spatstat'. __ 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-projec

Re: [R] choosing a particular object

Thanks, it's exactly what I want. Adrian On Friday 31 March 2006 12:59, Adaikalavan Ramasamy wrote: > Try > > test.fn <- function(obj.name, var.name="q2"){ > > stopifnot( is.character(obj.name) & is.character(var.name) ) > x <- subs

[R] choosing a particular object

function(obj.name, var.name) { temp <- give.me.the.object.called(Object) table(temp[, var.name]) } This should perfom the same thing as: table(bb$q2) Is this possible? TIA, Adrian -- Adrian DUSA Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5

Re: [R] ifelse problem

On Friday 10 March 2006 16:31, Uwe Ligges wrote: > [...] > aa[!(nchar(aa) < 3)] Thanks very much, I got it now. All the best, Adrian -- Adrian DUSA Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 3120210

[R] ifelse problem

ifelse function return the whole aa vector? Using if and else separately, I get the correct result... > if (any(nchar(aa) < 3)) { aa[-which(nchar(aa) < 3)] } else { aa } [1] "test" "name" Thanks in advance, Adrian -- Adrian DUSA Roma

Re: [R] QCA adn Fuzzy

ions and/or feedback are more than welcome. I hope this helps you, Adrian -- Adrian DUSA Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 3120210 / int.101 __ R-help@s

[R] about"Riply's K function"and "envelope"

explained in Kest.fft. I apologise that the function Kest.fft is not mentioned in the help page for envelope. We will fix that. Also please note that there is a bug in spatstat 1.8-6 which causes some functions to run slower. It is fixed in 1.8-7 which will be released this week. regards Adrian Badd

[R] weird behavior of nsmall in format

dered if the current implementation of nsmall in format is correct. Thanks, Adrian __ 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

[R] how to clear the last error?

Hello, I would like to clear the last error that I get with geterrmessage(). Not even rm(list=ls()) clears it. Can I set it to NULL somehow? Thank you, Adrian __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] adding variable into dataframe by indice

Thank you for your assistance. I find that the by using the list I learn something new everytime. --Adrian Hi not sure if I understand correctly but table() can be used ttt <- table(asubs112$fir) prop <- ttt/nrow(asubs112) asubs112$prop <- prop[match(asubs112$fi

[R] adding variable into dataframe by indice

dce,] } Below is the output that I get from the script above. Notice the proportion for the first NA but not any of the others. Not sure what I am doing wrong, any suggestions are a big help. TIA, Adrian asubs112[1:50,] IND_ID rs1042364 first_drink age_int V5

[R] insert value according to indice

(indce)/nrow(asubs112) asubs112[indce,5] <- prop asubs112[indce,] } Below is the output that I get from the script above. Notice the proportion for the first NA but not any of the others. Not sure what I am doing wrong, any suggestions are a big help. TIA, Adrian

[R] [R-pkgs] new package QCA

ethod"). Future versions of this package will have more functions to address the fuzzy-set minimization problems, as well. Big thanks to the r-help list members, supportive as ever, especially to Gabor Grothendieck and Martin Maechler for excellent ideas in the key parts of the algorithm. -- A

Re: [R] yet another vectorization question

. It seems that creating the necessary row and column indexes to cbind is much slower than copying chunks of data at certain columns: > library(combinat) # for the combn function > system.time(all.expr(LETTERS[1:12])) [1] 6.12 0.39 6.54 0.00 0.00 > system.time(all.expr2(LETTERS[1:12])

Re: [R] yet another vectorization question

aa[, i]] <- c(0, 1, 0) Is there any possibility to vectorize this "for" loop? (sometimes I have hundreds of columns in the "aa" matrix) Many big thanks in advance, Adrian -- Adrian DUSA Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Te

Re: [R] yet another vectorization question

and.grid(lapply(1:3, function(x) c(0, 1, NA))) Best, Adrian -- Adrian DUSA Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 3120210 / int.101 __ R-help@stat.math.ethz.ch m

Re: [R] yet another vectorization question

) and write the contents of "tt" in the "result.matrix"... Best, Adrian -- Adrian DUSA Romanian Social Data Archive 1, Schitu Magureanu Bd 050025 Bucharest sector 5 Romania Tel./Fax: +40 21 3126618 \ +40 21 3120210 / int.101 ___

Re: [R] yet another vectorization question

Adrian DUSA roda.ro> writes: > > I'm trying to develop a function [...snip...] Sorry for the traffic, I forgot to say that I'm using library(combinat) for the "combn" function... Thank you, Adrian __ R-help@stat.ma

[R] yet another vectorization question

ol(idk)) { end.row <- start.row + nrow(tt) - 1 return.matrix[start.row:end.row, idk[ , k]] <- tt start.row <- end.row + 1 } ## How can one modify "return.matrix" using apply on "idk"? } return.matrix[is.na(return.matrix)] <

[R] envelopes of simulations

You can compute the envelope for the K-function and transform it afterwards. To plot envelopes of the L function discrepancy LD(r) = sqrt(K(r)/pi) - r, just type E <- envelope(PPPData, Kest, correction="isotropic")

  1   2   3   >