similar to: lme convergence

Dear All: I would like to center the data in 'x' by 'group'. The following code scale the data and I have not been able to figure out how to change it so I get the centered data. x <- c(1, 2, 3, 4, 5, 6, 7, 8) group <- c(1,1,1,2,2,2,2,2) unsplit(lapply(split(x,group),scale),group) I would appreciate your help. Ashraf

Hi All: I have to compute bootstrap confidence intervals, the statistic (incremental cost effectiveness ratio) is computed from two samples (intervention and control) of different sizes. All the bootstrap functions that I have seen use one dataset as argument. I may go ahead and get the desired number of bootstrap replications separately. I would appreciate if you could point me to a source of a

I am trying to find out if R can recognize specific criteria for removing rows (i.e. a prexisting function) I have a matrix myMatrix that is 12000 by 20 I would like to remove rows from myMatrix that have: -999 across all columns -999 across all columns but one -999 across all columns but two -999 across all columns but three -999 across all columns but four -999 across all columns but five

Hi all: I get the following error message that I am not able to resolve. Error in if (const(t, min(1e-08, mean(t)/1e+06))) { : missing value where TRUE/FALSE needed It appears right before the last data.frame statement. Below is the program that simulates data from one way random effects model and then computes normality and bootstrap confidence interval for

Dear All: I have been struggling to run the conditional logistic regression on a dataset. I am using clogit function in survival package. Here is how the data is generated. I tried to shorted the data by using the second count variable but still the clogit did not run. Why I am not able to run clogit on my dataset. I appreciate any input. Regards, Ashraf age <- c(rep(0,8),rep(1,8)) scc <-

I am producing 2X2 correlation matrices by a class variable. I need to extract a vector of correlation coefficients only. I am doing that in a loop (see below) but I am sure there would be a simpler way. Please help! > by(d1[,c(2,3)],d1[,1],cor) d1[, 1]: 1 c e c 1.0000000 0.1972309 e 0.1972309 1.0000000

Hi All: I would like to generate a binomial random variable that correlates with a normal random variables with a specified correlation. Off course, the correlation coefficient would not be same at each run because of randomness. I greatly appreciate your input. Ashraf

On Thu, Jan 31, 2013 at 10:29 AM, Weijia Wang <wwang.nyu@gmail.com> wrote: > Hi, > > > > I have a new question about subsetting in R. > > > > Say we have this data frame: > > > > PT_ID Blood_Pressure OBS_TYPE > > 92 1900 90.0 DBP > > 94 1900 90.0 DBP > > 174 2900 140.0 SBP > > 176 2900

Hello, I do have two different plots from LAD and OLS regression objects like this: # LAD regression and related plot rq(formula = SBP ~ Age) f = coef(rq(SBP ~ Age)) pred = f[1] + f[2]*Age plot (Age, SBP) lines (Age, pred) # OLS regression and related plot Pred = lm(SBP ~ Age) plot (Age, SBP) lines (Age,fitted(Pred)) Well for comparatative reason I would would like to subsume both plots

>>> Jeff, Patrick: Could you try my code from my github repo branch efi-multinic?? It's derived from Patrick's code and I finally see good responses with a VMware VM's e1000e NIC (never saw ANYTHING good from it until now). git://github.com/geneC/syslinux.git https://github.com/geneC/syslinux.git -- -Gene <<< Hi there I think in the case of a particular

Dear R communtiy, I do not understand why this does not work...: > betaS$SBP [1] 0.03274 -0.04216 -0.08986 -0.45980 0.60320 -0.63070 -0.05682 0.20130 > t<-c("betaS$SBP") > t [1] "betaS$SBP" > get(t) Error in get(t) : object 'betaS$SBP' not found [I am trying to use the variable "t" in a loop to call many different objects, but the

>>> > I think these changes would solve the thing. > > ... > -EFI_SERVICE_BINDING *sbp; > +EFI_SERVICE_BINDING *sbp =NULL; > ... > > -? ? if (status != EFI_SUCCESS) > +? ? if (status != EFI_SUCCESS || sbp == NULL) >? ? ? ? goto free_binding; > ... The presumption was that we'd see an error rather than non-error from LibLocateHandle()

Dear all, I have found some difference of the results between multinom() function in R and multinomial logistic regression in SPSS software. The input data, model and parameters are below: choles <- c(94, 158, 133, 164, 162, 182, 140, 157, 146, 182); sbp <- c(105, 121, 128, 149, 132, 103, 97, 128, 114, 129); case <- c(1, 3, 3, 2, 1, 2, 3, 1, 2, 2); result <- multinom(case ~ choles

Dear R help, I have the following data frame: structure(list(prochi = c("ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1", "ind_1"), date_1st_event = structure(c(14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784, 14784 ), class = "Date"),

For an analysis of data from twins, I'd like to do some rearranging in my dataframe. Currently, each twin has his or her own row. I'd like to add several columns of the co-twin's data to each row. In other words, this is the current situation: FamilyID ParticipantID IQ Digit_span 1 1 95 6 1 2 93 7 2 3 102 8 2

On Sat, Aug 15, 2015 at 9:50 AM, Gene Cumm <gene.cumm at gmail.com> wrote: > On Jul 27, 2015 12:30 PM, "Patrick Masotta" <masottaus at yahoo.com> wrote: >> >> >>> >> > I think these changes would solve the thing. >> > >> > ... >> > -EFI_SERVICE_BINDING *sbp; >> > +EFI_SERVICE_BINDING *sbp =NULL;

Hi All; The following question is directed more to statisticians on the list. Suppose I have a normal random variable. Is there a transformation that I can apply so that the new variable is slightly positively skewed. I greatly appreciate your help. Ashraf

# uname -a FreeBSD ck.com 7.0-STABLE FreeBSD 7.0-STABLE #1: Tue Jul 15 11:38:41 EDT 2008 ck@friendlybearonskates.com:/usr/obj/usr/src/sys/CK i386 # Make buildkernel kernconf=CK fails here : linking kernel.debug sbp.o(.text+0xd9d): In function `sbp_free_ocb': /usr/src/sys/dev/firewire/sbp.c:2905: undefined reference to `xpt_release_devq' sbp.o(.text+0x123e): In function

Dear all, When calling rlm with the following data, I get an error. (R v.1.8.1, WinXP Pro 2002 with service pack 1.) > d <- na.omit(data.frame(CPRATIO, HEIGHTZ, FAMILYID)) > c <- tapply(d$CPRATIO, d$FAMILYID, mean) > h <- tapply(d$HEIGHTZ, d$FAMILYID, mean) > c 1 2 3 6 7 9 10 11 6.000000 2.500000 3.250000

Hi, I am trying to get FreeBSD 4.7 to connect to a Maxtor 5000XT drive, but I get the following when connecting it -> Apr 26 09:19:43 chowder /kernel: fwohci0: node_id = 0xc000ffc1, CYCLEMASTER mode Apr 26 09:19:43 chowder /kernel: firewire0: 2 nodes, maxhop <= 1, cable IRM = 1 (me) Apr 26 09:19:50 chowder /kernel: fwohci0: BUS reset Apr 26 09:19:50 chowder /kernel: fwohci0: node_id =