similar to: REPLACE and REMOVE fns

I am calculating batting averages using the attached code ("dye.02_R_code.doc" and "giambi.03_R_code.doc") and the attached data. I've run this code for several players' data but always get the same results (see attached "MLB R CONSOLE.doc). It's obvious that the results wouldn't be exactly the same. Any help would be greatly appreciated. *I've

Hi, I cant understand where I am going wrong.Below is my code.I would really appreciate your help. Thanks. > genfile<-read.table("c:/tina/phd/bs871/hw/genfile.txt",skip=1) > > #read in SNP data > snp.dat <- as.matrix(genfile) > snp.name <- scan("c:/tina/phd/bs871/hw/genfile.txt",nline=1,what="character") Read 100 items

Dear R helpers, I need to find out maximized log likelihoods, parameters estimates and standard errors (in parentheses) of r largest-order statistics model, with different values of r by using the function rlarg.fit. I want to specify required number of order statistics to the model. I attached my data file with this mail.please help me. Ruposh --- r-help-request at stat.math.ethz.ch wrote:

Saludos, estoy intentando usar la función reprojectHDF() ( http://r-gis.net/?q=ModisDownload). #### source('ModisDownload.R') library(raster) library(RCurl) path<-"~/R/MODIS/Data/Test" input<-list.files(path,"*.hdf",all.files=T,recursive=T,full.names=T) input outname<-paste(substr(input[1],40,45),".tif",sep='') outname

Hi, How to make an if condition work, if the condition inside if() is created dynamically ,and that is a string .If i type teh dynamically created string the if works fine but when dynamically created,it is a string and going inside the if() ,an error is thrown saying : rgument is not logical ..I even tried changing teh mode of the string to logical,but it doesn't work Say my dynamically

Hi all, I have the following code: set.seed(1) x1 <- matrix(sample(1:12), ncol=3) x2 <- matrix(sample(1:12), ncol=3) x3 <- matrix(sample(1:12), ncol=3) X <- list(x1,x2,x3) tt <- matrix(round(runif(5*4),2), ncol=4) Is there a way I can construct a new list where newlist[[i]] = tt[i,] %*% X[[i]] without using a for loop? Each element of newlist will be 3 x 1 vector. Thanks -- Tina

On Wed, Jan 11, 2012 at 04:52:36PM +0800, Zhenzhong Duan wrote: > add polling interface to xen-netfront device to support netconsole > Ian, any thoughts on the spinlock changes? > Signed-off-by: Tina.Yang <tina.yang@oracle.com> > Cc: Konrad Rzeszutek Wilk <konrad.wilk@oracle.com> > Cc: Jeremy Fitzhardinge <jeremy@goop.org> > Signed-off-by: Zhenzhong.Duan

We have an occasional problem which manifests with multiple processes being created for a particular user. For example, for a user "xy004": xy004 8463 0.0 0.0 6.02M 0K ?? IW 11:48:59 0:02.02 smbd xy004 9426 0.0 0.0 5.99M 0K ?? IW 11:52:18 0:03.68 smbd xy004 10433 0.0 0.0 5.81M 0K ?? IW 12:17:20 0:00.85 smbd xy004

I am trying to calculate Bayesian Credible Intervals for a proportion (disease prevalence values to be more specific) and am having trouble using R to do this. I am working with ncredint() function but have not had success with it. Please help! Example: Positive samples = 3 Total sampled = 10 Prevalence = 0.3 pvec <- seq(1,10,by=1) npost = dbinom(pvec,10,prob=0.3, log=FALSE) ncredint(pvec,

I have a list 'ans' from the following code: tt <- rnorm(50) rr <- rnorm(50) ans <- lm(rr~tt) ans[1] is "$coefficients", ans[2] is "$residuals", ans[3] is "$effects", ... and so on up to ans[12]. Is there an easy way to display just these names and not the data they contain? I thought I saw my advisor type "ans$" and they were displayed,

I was writing a simple plugin a while ago which just sets the paint values on inactive windows. The plugin is very simple and the main part is below. When I load this plugin it blocks dbus and fuse plugins so that they only reply to requests when the active window changes. Is there anything obvious here that would cause that? static Bool inactivePaintWindow (CompWindow *w, const

Hello, I'm trying to convert a character column in several dataframes to lower case. ### # # Sample data and 'spp' column summaries: # dput(ban.ovs.1993[sample(row.names(ban.ovs.1993), 20), 1:4]) ban.ovs.93 <- structure(list(oplt = c(43L, 43L, 38L, 26L, 35L, 8L, 39L, 1L, 34L, 50L, 10L, 29L, 31L, 24L, 18L, 12L, 27L, 49L, 28L, 51L), rplt = c(NA_integer_, NA_integer_,

Hello Everyone Both the MCMCpack and the bayesm libraries allow us to make draws from the Inverse Wishart distribution. But I wanted to find out how exactly is the Inverse Wishart distribution parameterized in these libraries. The reason I ask is the following: Now its generally standard to express Inverse Wishart as IW(0.5 * DOF,0.5* Scale). (DOF-> Degree of freedom, Scale -> Scale

I tried toy problems and there doesn't seem to be a basic problem between lm() and ts objects: X = data.frame(x=c(1,2,7,9), y=c(7,2,3,1)) lm(y ~ x, X) X <- lapply(X, function(x) ts(x, frequency=12, start=c(1994,7))) lm(y ~ x, X) and this works fine - whether you do an lm() before or after making ts objects, it's okay. But I have a situation where things aren't okay.

Hello R users, I am trying to run a model with a binary response variable (nesting success: 0 failure, 1 success) and 8 fixed terms. Nesting success was examined in 72 cases in 34 territories (TER) during a 6 study years. Territories are nested within 14 patches (PATCH). I want to run a model taking into account these nested factors and repeated observation. To do this, I assume that the best

Dear R list, I have the following data: set.seed(1) n <- 5 # number of subjects tt <- 3 # number of repeated observation per subject numco <- 2 # number of covariates x <- matrix(round(rnorm(n*numco),2), ncol=numco) # the actual covariates x > x [,1] [,2] [1,] -0.63 -0.82 [2,] 0.18 0.49 [3,] -0.84 0.74 [4,] 1.60 0.58 [5,] 0.33 -0.31 I need to form a matrix

Samba not responding to broadcast calls from networked PCs. Short: My Redhat Linux 8.0 box is not replying to udp broadcasts, and I can't figure out why. Network: The offending Redhat 8.0 Linux Box (//DS) attempting to run Samba 2.2.7 Firewall Redhat 7.3 (//IW) successfully running Samba 2.2.4 Win98 (//MS) Win2k (//SL) <Not running at the time of these tests> Very long: I've

Hola, paso el mini programita q estoy viendo, lo q me llama la atencion es una parte donde se definen las funciones. Probe primero meter adentro del boots la estadistica a estimar usando directamente gini(varible, pesos) pero no me dejo. Vi q en el ej del manual de boots, siempre define antes la funcion, entonces probe definir antes una funcion haciendo grini<-function(x) {gini(variable,

I can't install freebsd 6.0 STABLE amd64 on my notebook: HP Pavilion 2420la AMD Turion ML 32 , 1024 RAM, 80 GB HD, ATI Radeon XPRESS 200M. I try all modes , default, whitout acpi, safe mode. In all the cases after the next message the installation doesn't work: timecounters tick every 1000 msec Linux ELF exec handler install lo0.bpf attached I installed linux with parameter noapic

chan_sip.c6524 reload_config= unable to get ip address from asterisk, sip disabled The ip address is working fine, Internet works great. Can anyone help...Thanks