similar to: unequal number of observations for longitudinal data

I am wondering if there is a function which will do a join between 2 data.frames by minimum distance, as it is done in ArcGIS for example. For people who are not familiar with ArcGIS here it is an explanation: Suppose you have a data.frame with x, y, coordinates called track, and a second data frame with different x, y coordinates and some other attributes called classif. The track data.frame has

Dear all, I have used 'read.table' to create a data frame of 720 columns and 360 rows (and assigned this to 'Jan'). The row and column names are numeric: > columnnames <- sprintf("%.2f", seq(from = -179.75, to = 179.75, length = 720)). > rnames <- sprintf("%.2f", seq(from = -89.75, to = 89.75, length = 360)) > colnames(Jan) <- columnnames

Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 --------------------------------- Colleagues Is there a function to calculate distances between adjacent latitude/ longitude pairs in a matrix? It is basically a spatial diff() function that I have in mind. Data: long1, lat1 long2, lat2 long3, lat3 looking for: diff(data$long, data$lat)

All - I have a large data frame that looks like ID p1 p2 p3...p20 Lat1 Lat2 Lat3...Lat20 Long1 Long2 Long3...Long20 1 0 0 1 0 NA NA 29.xx NA NA NA -89.xx NA 2 1 0 0 1 27.xx NA NA 29.00 -88.00 NA NA -89.xx 3 0 0 0 0 NA

Deemed colleagues I would appreciate your help with a sentence to calculate the linear distance between 2 geographical points (coordinates in UTM). In advance thnks for your attention, -- Ricardo Bandin Llanos rbandin@udec.cl Estudiante - Magíster Cs. m. Pesquerías Universidad de Concepción, Región del Bio-Bio, Chile Celular: (0056-41) 97949957 [[alternative HTML version deleted]]

Apologies for troubling the list, but it is a test that needs to be carried out. I am resending the message that I sent earlier on behalf of Paul Miller, but with a certain word used in the variables names of the SAS code changed to a different word. With thanks for your tolerance ... Ted. [*** PLEASE NOTE: I am sending this message on behalf of Paul Miller: Paul Miller <pjmiller_57 at

Hello list, I have defined a new MoH-class in musiconhold.conf : [default] mode=files directory=/var/lib/asterisk/moh random=yes ; *[106002] mode=files directory=/var/lib/asterisk/moh/106002 random=yes* In sip.conf I have this commented out : ;mohinterpret=default ;mohsuggest=default Asterisk sees these moh-classes and files : vps2301*CLI> moh show classes Class: default Mode: files

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.

I have the following longitudinal data: id time y 1 1 10 1 2 12 1 3 15 1 6 18 2 1 8 2 3 9 2 4 11 2 5 12 3 1 8 3 4 16 4 1 9 4 5 13 5 1 7 5 2 9 5 6 11 .... I want to select the observations at time 4. if the observation at time 4 is missing, then i want to slect the observation at time 3. if the observation at time 3 is also missing, then i want to select observation at time 5. otherwise i will

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

Just out of curiosity, I took the default "iris" example in the RF helpfile... but seeing the admonition against using the formula interface for large data sets, I wanted to play around a bit to see how the various options affected the output. Found something interesting I couldn't find documentation for... Just like the example... > set.seed(12) # to be sure I have

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

I ask for a histogram of counts with slightly uneven breaks but histogram() ignores the 'type' argument and provides density scaling instead. x = sample(1:3, 100, replace=TRUE) histogram( ~ x, breaks=c(0,1.5,2.5,3.5), type='count') My real application has time data with month boundaries for the breaks and I DO want counts on the y-axis. A work-around is not difficult but

Dear R folks - Sorry to be so dense, but can you help me with two programming problems? 1. Randomizing within factors. I have factored the values in an array according to a series of ranked breakpoints. Now I want to randomize the values that fall within each of the breakpoint boundaries, without replacement. I want to end up with an array in which the values within each factored category

Hi I am trying to run an ANOVA on data with unequal variance. I am new to nlme, but to my understanding I need to use the gls function. I have single response variable (distance which is continuous) and the explanatory variable is individual ID (class variable: individuals differ in the variance in their distance values hence the need to using nlme). So I would create a model

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>

Hi, i need to test the equality of proportions in a paired case, like here: after before + - + 10 14 - 5 53 So usually i use the mcnemar.test(stats). Due to mortality, I now got the problem that my sample sizes in both factors are not equal any more -- there are less cases in the "after" condition. I already learned that Ekbohm (1982) and Marascuilo adressed the problem

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

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

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