search for: ay

...e real dataset is of cause much larger, app. 100.000 samples I need to adjust each value to remove any tendency of some years generally having higher values and others lower, since this is an unwanted artifact from different measuring traditions. My plan is to generate an average for each year Ay, as well as a global average Ag. Then each value should be multiplied by Ay/Ag. I can make the averages like this: > Ag <- mean(dafSamp[,2]) > Ag [1] 95.11111 > Ay <- aggregate(x=dafSamp[,2], by=list(dafSamp[,1]), FUN='mean') > Ay Group.1 x 1...

Dear list members, Any help on this efficiency issue would be greatly appreciated. I would like to find the most efficient way to run a non-vectorized function (here: fisher exact test p-value) iteratively using 4 matrices with identical dimensions. And as a result I aim for an array with identical dimensions containing the corresponding p-values. Please consider some code using a trivial example with 3x4 arrays below. Eve...

...of several textbooks and the R list serve. I have three treatments of interest (A,B,C), each with 2 treatment levels. One level of treatment A (Y or N) was applied to 6 blocks (Block, labelled 1,2,3,4,5,6). Each block was subdivided into 4 subplots (Positon, labelled 1,2,3,4, so implicit nesting may be a problem), to which treatment B (Y or N) was randomly assigned (2 replicates of each level of treatment B per block). Each of these subplots was divided in half and randomly assigned treatment C (B or M). My response variable is the number of predetermined spots (out of 132) occupied by a cer...

Hi, I have a large data set that I input as a matrix, where I have 1:x rows with names AX, then x+1: y rows named AY, etc. The idea is that I have to count how many rows exactly I have with name AX and how many I have with name AY (or find which row numbers have names AX). Is there any way in R to count a number of rows with a name matching a required pattern? Thanks in advance, Katie -- View this message i...

Finally got a chance to install Centos 6.x x86_64 version and having an issue with the latest version of freenx. This is a desktop install on a esxi 5 with all updates and only thing I have installed is the latest vmware tools. When I enable the extras repo and install the nx-3.5.0-1.el6.ay.x86_64.rpm and freenx-0.7.3-8.el6.ay.x86_64.rpm. When I try to logon the desktop pops up but I get the little busy circle and it never completes, I am unable to close the session or do anything to the freenx window. Here is the odd, part, if I install the earlier versions from toracat, freenx-0.7...

...he +# search path than package:base, base methods appeared to +# preferentially call the original base versions, rather than the +# patches versions that I wanted. So, big hammer it, override +# everything - which will probably make it easier to contribute these +# improvements back to stock R anyway. +# +# --atp at piskorski.com, 2005/10/01 02:29 EDT +# +# $Id: all.equal.R,v 1.4 2005/10/01 13:10:25 andy Exp $ + + +# In S-Plus, all.equal() prefers to index objects by name, while in +# stock R, it prefers to index by position. IMO, *NEITHER* of those +# behaviors are fully correct. What we rea...

The function of crossprod in R puzzled me. I would like to calculate the cross product of two vectors. According to my text book, it defines like this: a = (ax, ay, az) b = (bx, by, bz) then, the cross product of a and b is: a X b = (ay*bz-az*by, az*bx-ax*bz, ax*by-ay*bz) It can also write in a determinant format. But the crossprod or tcrossprod function in R appeared not calculate the cross product. Suppose I enter this: > a = c(1, 2, 3) > b = c(2, 3...

Hello! I've been looking for a way about how to change the emulated 'windows' standard terminal font now quite for a while. I have a windows program which normally uses the standard terminal-font of Windows as the font for its 'text-output'. If have tried to copy the 'Windows/Fonts'-folder of my original wi...

...are independent.   Suppose A represents (7-7.50), B represents (7.50-8.50) and C represents (8.50 – 10.00). Also let X be (11-12), Y is (12-14) and Z is (14-21).                         These two groups are independent i.e. Prob(A and Y) = P(A) * P(Y)   So there are 9 combinations possible as (AX, AY, AZ, BX, BY, BZ, CX, CY and CZ) respectively with the joint probabilities (0.059, 0.235, 0.126, 0.031, 0.123, 0.066, 0.05, 0.202, 0.108) respectively.   My problem is (i)                               How to obtain these 9 combinations of probabilities in the sense how do I obtain the various combi...

Dear list, I'm seeking advice to extract some numeric values from a log file created by an external program. Consider the following example, input <- readLines(textConnection( "some text <ax> = 1.3770E-03 <bx> = 3.4644E-07 <ay> = 1.9412E-04 <by> = 4.8840E-08 other text <aax> = 1.3770E-03 <bbx> = 3.4644E-07 <aay> = 1.9412E-04 <bby> = 4.8840E-08")) ## this is what I want results <- c(as.numeric(strsplit(grep("<ax>", input,val=T),...

I renamed my volume with vgrename however I didn't complete the other steps. Mainly update fstab and intiramfs. Once I booted, I was dropped on the Dracut shell. From here I can see the newly rename VG and I can lvm lvscan as well as activate it, lvm vgchange -ay. However I can't figure out what to do next, I'm assuming I need to regenerate the initramfs and then boot to change grub? Could someone point me in the right direction to recovering a FS from Dracut, or other means, once the volume group name was changed. TIA, Paul

I have a degraded raid array (originally raid-10, now only two drives) that contains an LVM volume. I can see in the appended text that the Xen domains are there but I don't see how to mount them. No doubt this is just ignorance on my part but I wonder if anyone would care to direct me? I want to be able to retriev...

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi I know there is a function - I have used it before - but I always forget what it is called... I need the combination of two character vectors, i.e: x <- c("a", "b") y <- c("x", "y") z <- THEFUNCTION(x, y) z == c("ax", "ay", "bx", "by") I promise I will write the name of...

Dear all, I'm looking for a way to replace NA's with 0 for a number of variables which share the same ending and which constitute a subset of a data frame. Say, for example, there is > data <- data.frame(ax=c(1,2,3,NA,0) , bx=c(3,5,NA,5,1) , ay=c(5,NA,2,NA,0) , by=c(4,2,NA,2,1)) > data ax bx ay by 1 1 3 5 4 2...

...problem with Raid controller Tekram TR-822... I'm using dmraid software: root at ns1: / $ dmraid -r /dev/sda: sil, "sil_ajacddafagff", mirror, ok, 156299440 sectors, data@ 0 /dev/sdb: sil, "sil_ajacddafagff", mirror, ok, 156299440 sectors, data@ 0 When I'm start dmraid -ay: root at ns1: / $ dmraid -ay RAID set "sil_ajacddafagff" was not activated And I see from dmesg: device-mapper: table: 253:0: mirror: Device lookup failure device-mapper: ioctl: error adding target to table My system centos My kernel 2.6.18-128.1.10.el5.centos.plus

Dear R helpers,   I am working on the scenario analysis pertaining to various interest rates. In this connection I need to form the various combinations as under :   Suppose I have two sets A = (a, b, c) and B = (x,y,z)   Then I can easily form the cominations as (ax, ay, az, bx, by, bz, cx, cy, cz)   However, if I have say 5 variables, then total no of possible combinations will be 3^5 = 243. Thus, A = (a,b,c), B = (x, y, z), C = (l, m, n), D = (p,q,r), E = (s, t, u). Then may be my possble combination will start as (a, x, l, p, s), then next combination may be (a...

Hi, I'm working on the codes below however every time I run them when they get to OpenBUGS I keep getting the error message: array index is greater than array upper bound for hab. Any help would be greatly appreciated, Suzie Codes: ungulate <- read.csv(file.choose ()) #ungulate ungulate <- as.matrix(ungulate);colnames(ungulate)<-NULL;rownames(ungulate)<-NULL habitat <- read.csv(file.choose ())...

Is there code for bivariate kernel density estimation? For bivariate kernels there is kde2d in MASS kde2d.g in GRASS KernSur in GenKern (list probably incomplete) but none of them seems to accept a weight parameter (like density does since R 2.2.0) -- Erich Neuwirth, University of Vienna Faculty of Computer Science Computer Supported Didactics Working Group Visit our SunSITE at

Hi, I am trying to recover data from my old system which had LVM. The disk had two partitions - /dev/sda1 (boot, Linux) and /dev/sda2 (Linux LVM). I had taken a backup of both partitions using dd. Now I am booting of CentOS live cd for system restore. I recreated partitions like previous system using fdisk and then used dd to dump all the data onto it. I would like to mount sda2 as LVM, but I

Hi, I have always been wondering why S-Plus/R can not fit a linear model to an arbitrary large data set given that, I thought, it should be pretty straightforward. Sometime ago I came across a reference to LM package, http://www.econ.uiuc.edu/~anovo/LM.html, by Roger Koenker and Alvaro Novo. So I thought here it is...