similar to: help with matrix creation

> Date: Thu, 10 Aug 2000 15:33:39 -0700 > From: Jesus Maria Frias Celayeta <iosu at bureau.ucc.ie> > > hi R-help! > > I would like to know a simple and easy way (if posible) to do the following > data manipulation. I have a matrix of experimental data (with replicae in time) actually, they seem to be data frames. > >

I'm trying to create a function which will allow me to subset a data set based on values of various specified variables. I also want to then apply some other function(s) (e.g., summary). This is what I've tried so far.... > test.fx <- function(dta, expvar, expval) { + newdta <- subset(dta, eval(expvar)>expval) + summary(newdta$eval(expvar)) + } > >

I haven't been an R lister for a bit, but I hope to enlist someone's help here. I think this is a simple question, so I hope the answer is not much trouble. Can you please respond directly to this email address in addition to the list (if responding to the list is warranted)? I have a histogram and I want to see if the data fit a Poisson distribution. How do I do this? It is

Hello, i only got a small problem. i try to create automatic new dataframes, or png?s. the main problem i got is: how can i create automatic a new name for a file (read out by simply "for") - i tried to use "(paste...) but theres an errormessage, about a wrong declination. R told it is as.character, but need as.Real. Should i use another method than "paste"? i tried as

I am trying to run this loop but it keeps coming up with the following error: 'Error: no function to return from, jumping to top level' n.simul <- 1000 A <- {for(i in 1:12){{ nT <- rpois(n.simul, T[i]) nP <- rpois(n.simul, P[i]) nD <- rpois(n.simul, D[i]) nC <- rbinom(n.simul, nT, C[i]) Score <- (5*nT + 2*nC + 3*nP + 3*nD) } return(Score)}} My aim is to calculate

Hello, I have a question on export and import data. Thank you for any suggestions. data 'simul' is generated as follows: N <- 20 n <- N/2 nsets <- 10 simul <- matrix(0,nsets,N) th <- c(0,1, 1) for(i in 1:nsets){ simul[i,] <- rnorm(N,mean= rep(th[1:2],N/2),sd=th[3]) } I exported data as follows: write.table(simul,

I have a question about the lme function in R. My question is: After I got the object from function lme, why the numIter value of the object is always NULL? Following is my code: jjww<-lme(y~x*zz,data=simul,random=~x|group, control=lmeControl(returnObject=TRUE)) attributes(jjww) jjww$numIter the first 20 observation of data simul are: > simul y

Hello, I have a data frame with different groups (Time, Target, Conc) and each entry has a triplicate value of the measurements OD and ODnorm. How can I merge the triplicates into a single mean value? I tried the following: ``` df = data.frame(Time=rep(1, 9), Well=paste("A", 1:9, sep=""), OD=c(666, 815, 815, 702, 739, 795, 657, 705, 663),

Hi listers, I've been looking for a procedure, but I am not succeding... I have a function that give multiple results... Then, I would like to simulate this function n times, so I need to save/keep the n multiple results, in order to calculate my desired statistics... I have already tried with the RETURN and LIST FUNCTION, but I am not getting it right... An example of what I am looking for

Dear R-users, I would like to ask whether it's possible (for sure it would be), to plot each rows (or columns) in different graphs and in the same figure region without using the function "par" and then struggling around with "axes" and labels etc. Luckily, I would always have "rows + columns = even number" and the same "ylim". The next one could be a

Hello, I would like to add to axis labels special caracters. Instead of writing :plot(simul, xlab="beta", ylab="sigma11") It would be great if I could write something as in LaTeX :plot(simul, xlab="\beta", ylab="\sigma_{11}") Is there a way to do that ? Thank you _________________________________________________________________ ues clics pour retrouver tout

Hello, I wanted to fit a linear mixed model to a data that is similar in terms of design to the 'Machines' data in 'nlme' package except that each worker (with triplicates) only operates one machine. I created a subset of observations from 'Machines' data such that it looks the same as the data I wanted to fit the model with (see code below). I fitted a model in

Hey gang, I've been able to use sipp to produce some call volume on our asterisk server. The server has no problems handling 50 simul calls. But then again, no RTP is being done. I tried to use the rtp echo ability of sipp but that doesn't seem to work right. I also setup a fake number in asterisk that when called by sipp, would dial another number via PRI, hoping that some 729

Hi, I think you're misunderstanding which set of variables go on either side of the formula. Is this what you're looking for? > aggregate(OD ~ Time + Target + Conc, data = df, FUN = "mean") Time Target Conc OD 1 1 BACT 1 765.3333 2 1 BACT 2 745.3333 3 1 BACT 3 675.0000 > aggregate(ODnorm ~ Time + Target + Conc, data = df, FUN =

Hi listers, I am working on a simulation... But I am having some troubles... Suppose I have a function A which produces two results (mean and variance)... Then I would like to simulate this function A with a function B many times using the results from function A For example: #Function A boot<-function(a,b,c){ mean_boot<-(a+b)/2 var_boot<-c #list(a=a,b=b,c=c) return(a) } Then I would

My recent migration from Courier-IMAP to Dovecot has been marred by these errors: Oct 24 12:12:40 server dovecot: IMAP(juha): Fixed index file ~user/Maildir/.Geekzone/dovecot.index: first_recent_uid_lowwater 431 -> 430 Oct 24 12:12:40 server dovecot: IMAP(juha): Corrupted index file ~user/Maildir/.Geekzone/dovecot.index: Duplicate header extension keywords Oct 24 12:12:40 server dovecot: child

Hi All, I have the following dataset (triplicates values of 5 independent measurements and duplicate vaues of a control): 1 2 3 4 5 C 181.8 58.2 288.9 273.2 290.9 53.9 120.3 116.8 108.9 281.3 446 39.6 86.1 148.5 52.9 126 150.3 My aim was to find if mean values of Samples 1 - 5 are

Dear useRs, I have written a function that implements a Bayesian method to compare a patient's score on two tasks with that of a small control group, as described in Crawford, J. and Garthwaite, P. (2007). Comparison of a single case to a control or normative sample in neuropsychology: Development of a bayesian approach. Cognitive Neuropsychology, 24(4):343?372. The function (see

Also, > aggregate(cbind(OD, ODnorm) ~ Time + Target + Conc, data = df, FUN = "mean") Time Target Conc OD ODnorm 1 1 BACT 1 765.3333 108.33333 2 1 BACT 2 745.3333 88.33333 3 1 BACT 3 675.0000 18.00000 (You might wish for "cbind(OD,ODnorm) ~ . - Well", but aggregate.formula is not smart enough for that.) -pd > On 24 Oct 2023, at

Hi I've two servers : 1. UK 2. Pakistan Pakistan * server has ISDN30. Pakistan(ISDN30) <====> UK ===> User Im planning to setup an IAX2 trunk between these two server ? so , how much bandwidth I need for 30 simul. calls ? Im planning to use G729 on both my server ? to support 30 calls over IAX2 do I've to change some setting during compile time or not ? pls suggest.