similar to: Random Sample with Unique function

Hi, If any one has time I need some help understanding the P-values given in the lmer output. Using AIC for model selection I find my minimal model is FOLLOW~MOVERSTATUS+DISTANCE however it appears DISTANCE is not significant at 95% confidence, see output quoted below. However, removing DISTANCE gives a higher AIC=433.5, therefore I will keep it in, but am confused as to what is adds to the

Dear R users, Simple question. Can anyone help with the code that would allow me to view only the variables who's correlation output is >0.8? This is the code I'm using to date >cor(data, method="spearman") Kind regards Krys -------------------------------------------- _________________________________________________________________ Save time by

Dear all, I have to select 122 stratified random samples from a population of >3900 cells. I have 41 strata and I have to draw a different number of samples from them(between 2 and 8). I have tried to apply the funtion strata following the instruction in the manual: strata(dataframe, stratanames=NULL, size, method=c("srswor"), pik,description=TRUE) but I get the error

OK, now i've to start to move the big part of my users from my old NT-like domains to my new AD domain. I've setup roaming profile in the new domain following the wiki (https://wiki.samba.org/index.php/Roaming_Windows_User_Profiles, 'using windows ACL') and for new profiles works like a charm. But i've tried to move/copy old profile to the new domain, and seems work, with

Hi R users. I have a problem with function strata in sampling packages. > st0 = strata(dom, stratanames="stratas", size=sample.size, method="systematic",pik, FALSE) Error in sort.list(y) : 'x' must be atomic for 'sort.list' Have you called 'sort' on a list? In previous version of R 2.9.1 and previous version of package sampling this code worked

Suppose I have a long format for a longitudinal data id time x 1 1 10 1 2 11 1 3 23 1 4 23 2 2 12 2 3 13 2 4 14 3 1 11 3 3 15 3 4 18 3 5 21 4 2 22 4 3 27 4 6 29 I want to select the x values for each ID when time is equal to 3. When that observation is not observed, then I want to replace it with the obervation at time equal to 4. otherwise just use NA. How can I implement this with a quick

Dear all, I want to draw a barplot with the following data: [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 110.0675 118.4167 117.16 109.735416 101.6312 104.0312 101.8263 99.38541670 114.2613 [2,] 0.0000 0.0000 0.00 1.658333 0.0000 0.6250 0.0000 0.05208333 0.0000 [,10] [,11] [,12] [,13] [,14] [,15] [,16]

Hi, > Say we have a large-sized vector of observation A and a time vector with > the same size. The latter vector basically is the time flag when ith > obervation of A occurs, hence A[i] occurs at T[i]. > > I have no trouble with doing: T[A <= 10 | A >=20] for example > > Now, what I am confused: I tried to pull out when the largest value of A > with: > >

Hi Guys, I have this part of a program: library(survival) Gastric <- cbind.data.frame(Gp=c(rep(1,45),rep(0,45)), ### 2nd gp 0 time=c(1,63,105,129,182,216,250,262,301,301,342,354,356,358, 380,383, 383,388,394,408,460,489,499,523,524,535,562,569,675,676, 748,778,786,797,955,968,1000,1245,1271,1420,1551,1694,2363,2754,2950,

Hi everyone, I want to have a group and dial multiple phones/lines simultaneously. If I use this Dial command: exten => 222,2,Dial(${TRUNKBP}/246&SIP/258&${TRUNKBP}/243,20,tTr) ... all phones ring just once, after that only the first one continues ringing and only that one can answer. Can anyone tell me why? thanks! Tomica -------------- next part -------------- An HTML

Dear R-users, Following is part of my data, where slide has 36 levels and block 48 levels. I have done boxplot for each slide on the same graph. There are outliers for each slide and I tried to use indentify functtion to identify outliers in such a way that when I click on an outlier or point, the points will be labelled by either their block or ID or by both but without success. How can I make

This ought to work: resultdt <- lapply(PGWide[, 240:389], function(x, ...) try(tdt(x, ...))) You can then check the class of each component to see which one failed. Andy From: Farrel Buchinsky > > "Berton Gunter" <gunter.berton at gene.com> wrote in message > news:008601c67097$de1b46e0$5bc4fea9 at gne.windows.gene.com... > > ?try > > > > as in

Hello ! I am attempting to switch from being a long time SAS user to R, and would really appreciate a bit of help ! The first thing I do in getting a large dataset (thousands of obervations and hundreds of variables) is to run a SAS command PROC CONTENTS VARNUM command - this provides me a table with the name of each variable, its type and length; then I run a PROC MEANS - for numerical

Hi All.. I am new to ruby on rails.I want set the unquie on two fields. Means in one model three fileds are there. 1.Id, 2.C-id 3.Q-id I want set the unquie on C-Id,Q-Id combindely on the model validations Can u give me advice? I tryed with statemwnt validates_uniqueness_of [:c_id, :qid] But it now giving the proper results Regards, Ramya -- Posted via http://www.ruby-forum.com/.

unique(c(1:10,1)) gives 1:10 (i.e. unique values), is there any method to get only 2:10 (i.e. values that are unique) ? _________________________________________________________________ Easily edit your photos like a pro with Photo Gallery. [[alternative HTML version deleted]]

Hi, The man page for base::unique.matrix() and base::unique.array() says that MARGIN is expected to be a single integer. OTOH the code in charge of checking the user supplied MARGIN is: if (length(MARGIN) > ndim || any(MARGIN > ndim)) stop(gettextf("MARGIN = %d is invalid for dim = %d", MARGIN, dx), domain = NA) which doesn't really make sense. As

I'm very new to R, so please forgive me if I just missed the answer in existing documentation... I have a data set with at least three columns, X, Y, and Z. I want to produce a chart where one axis shows all the unique values of X, and the other axis shows all the unique values of Y. Each cell within the chart should contain the result of applying an aggregate function (such as mean(), for

hello! my question conserns with use of "panel" package (written by R.C.Gentlman) (unfortunately the manual and help sites are very short) 1. is it possible to do analysis just without a(ny) covariate? i suggest do it by introducing a covariate with level=0 in all obervations, this because of Q(z)=Q_o exp(beta*z), but it seemingly doesn't work 2. in the option gamma in the

I have 1000 oberved binary reponses Y, and I also have 1000 predictited probabilities predicted_prob (between 0 and 1) for those observed reponses. How can I plot a gain chart in R? Thanks. -- View this message in context: http://www.nabble.com/how-to-plot-a-gains-chart-in-R--tp25083628p25083628.html Sent from the R help mailing list archive at Nabble.com.

Colleagues, I am trying to de-duplicate a large (long) database (approx 1mil records) of diagnostic tests. Individuals in the database can have up-to 25 observations, but most will have only one. IDs for de-duplication (names, sex, lab number...) are patchy. In a first step, I am using Andreas Borg's excellent record linkage package (), that leaves me with a list of 'pairs' looking