similar to: Reading CSV file with unequal record length

Hello, I am working on the frailty model using coxph functions. I am running some simulations and want to store the variance of frailty (theta) values from each simulation result. Can anyone help me how to extract the theta values from the results. I appreciate any help. Thanks Shankar Viswanathan

Dear R helpers   Suppose e.g. I have a csv file having three variables defined and each of these variables have data items of say 40, 50, 45 length. When I open this csv file in 'R', I get 10 trailing 'NA's under first column and 5 'NA' s in case of 3rd column.   How do I get rid of these NA's s.t. when I read the first column, there should be only 40 data items, 2nd

I am new to R. I have a data set like this (given below is a fictional dataset): AgeCat FINWT 1 98 2 62 1 75 3 39 4 28 2 47 2 66 4 83 1 19 3 50 I need to calculate the weighted distribution of the variable AgeCat. In SAS i can do: proc freq data=ageval; tables agecat; weight finwt; run; What or is there an equivalent in R? TIA, Krishnan -- Krishnan Viswanathan 1101 High Meadow Dr

I have the following sample dataset (CSV input here:http://goo.gl/YR8LP. CSV output here: http://goo.gl/EFCC8) which I want to transform as follows. For each person in a household I want to create two new variables OrigTAZ and DestTAZ. It should take the value in TripendTAZ and put that in DestTAZ. For OrigTAZ it should put value of TripendTAZ from the previous row. For the first trip of every

Hi, I was wondering if multiple remote agents can access the same libvirtd. How do we prevent conflicting commands from the different agents. Can different access permission be set for the agents? (The documentation usually goes over the case where one agent can access multiple libvirtd and not this case) Thanks Arvind -------------- next part -------------- An HTML attachment was scrubbed... URL:

I am trying to recode a variable into another variable and while the package 'car' works well when it is only recoding A into B, I am not sure I can do the same with recoding (A or C) into B. If i can use recode please advise on how to. So i am using an if/else if conditions. My sample dataset is below along with the code and the warning and results i get. TIA Krishnan #****Code****

Hello- I have run an ANOVA on 4 treatments with unequal sample sizes (n=9,7,10 and 10). I want to determine where my sig. differences are between treatments using a Bonferroni test, and have run the code: pairwise.t.test(Wk16, Treatment, p.adf="bonf") I receive an error message stating that my arguments are of unequal length: Error in tapply(x, g, mean, na.rm = TRUE) :

Dear Douglas and R-help, Does lme assume normal distribution AND equal variance among groups like anova() does? If it does, is there any method like unequal variance T-test (Welch T) in lme when each group has unequal variance in my data? Thanks, Shirley

Hi, I am doing multiple regression, and have several X variables that are categorical. I read that I can use dummy or contrast codes for that, but are there any special rules when there're unequal #observations in each groups (4 females vs 7 males in a "gender" variable)? Also, can R generate these codes for me? THanks.

Hello, I have some patient data for my masters thesis with three groups (n=16, 19 & 20) I have completed compiling the results of 7 tests, for which one of these tests the variances are unequal. I wish to perform an ANOVA between the three groups but for the one test with unequal variance (<0.001 by both bartlett and levene's test) I am not sure what to do. I thought i would run

I am using levelplot() from lattice with grids that have unequal cell sizes. This means that the boundary between two cells is not always half-way between nodes, as levelplot() assumes. The result is that some cell sizes are rendered incorrectly, which can be painfully obvious if using relatively large cells. Is there any work-around? I am using the conditioning capability of lattice and

Dear R users, I need to merge two data frames based on both equal and unequal comparisons. The "sqldf" package used to work well , but today, I cannot resolve the following error by reinstallation of the sqldf package. Can anyone suggest a different way to perform this kind of merge function? Thank you, Ding > DMRlog2pbde47DMS <- sqldf("select * from DMR_log2pbde47 as a

Hi, i'm looking for an R function to fit a one-way ANOVA with one factor containing 10 levels. The factor levels have different numbers of observations (varying between 20 to 40). For most of the dependent variables i'm testing there are unequal variances among the factor levels. I see the function oneway.test: oneway.test(variable ~ factor, data=dataset) which by default does not

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110828/ac62948a/attachment.pl>

Hello All, Would anyone be able to help me understand how R computes a quantile-quantile plot for comparing two data samples with unequal sample sizes? Normally, the procedure should be to rearrange the larger data sample into n equally-spaced parts using interpolation, where n is the sample size of the smaller sample, and then plot the matching data pairs. I tried using different plotting

This is in followup to a thread on R-help with subject "Sampling". I claim that R does the wrong thing by default when sampling with unequal probabilities without replacement - the selection probabilities are not proportional to 'prob', for any draw after the first: I suggest that R do what S-PLUS now does (though you're free to choose a better implementation). What S-PLUS

Hi, I am having problems getting cohen's kappa to work. I have been using the function: ><-ckappa(x,y) from the psy package. I am trying to test for inter-observer reliability, I have 2 observers and 26 categories, however, the two observers might not necessarily have the same range of categories (I have unequal score ranges). However, I thought R could cope with this. Each time I

I am getting a repeated error when I try to run a logistic regression in R 2.8.1 >(glm(prop1~x1,data=glm1,family=binomial("logit"),weights=nt1)) Error in model.frame.default(formula = prop1 ~ x1, data = glm1, weights = nt1, : invalid type (list) for variable 'x1' x1 is multistate categorical (3 categories). 2 of the categories have 12 observation, one has 9. Is this what

Well fellow R users, I throw myself on your mercy. Help me, the unworthy, satisfy my employer, the ungrateful. My feeble ramblings follow... I've searched R-Help, the R Website and done a GOOGLE without success for a one way ANOVA procedure to analyse data that are both non-normal in nature and which exhibit unequal variances and unequal sample sizes across the 4 treatment levels. My

Hello, I am hoping someone could shed some light on power calculations for me. I have two small data sets of unequal sample size after NA removal (m = 5, f = 7). m <- c(2.0863, 2.1340, 2.1008, 1.9565, 2.0413, NA, NA) f <- c(1.8938, 1.9709, 1.8613, 2.0836, 1.9485, 2.0630, 1.9143) In a R help message/reply from Sep 30, 2001, it was noted that the "power.t.test" function