similar to: calculate slope of line

Dear All, I have data in HDF file format and would like to read it into R. I have tried the package hdf5 without success. Any ideas and suggestions?? Kind regards, Katrin -- Katrin Fleischer Vrije Universiteit Amsterdam Faculty of Earth and Life Sciences Subdepartment Hydrolgy and Geo-Environmental Sciences Room E-360 De Boelelaan 1085 1081 HV AMSTERDAM Tel: +31 20 59 87391

Dear All, I have a raster map of the class 'SpatialPointsDataFrame' and coordinates of the class 'SpatialPoints'. I would like to retrieve the values that are contained in the raster map at the specific locations given by the coordinates. Can anyone help me out? Kind regards, Katrin Fleischer

Dear all, I have built glm models based on presences/absences and a number of predictor maps and would like to compute habitat suitability based on the modelled coefficients. I thought this is pretty straight forward and wanted to use predict() and supply the new data in a data frame, with one column for each predictor. However, I do get an error msg warning me that the number of rows for

Hello everyone, I have the following problem: My analysis includes many predictor variables (>50) in the form of raster maps (asc), but I am trying to avoid having to type all their names over and over again in the analysis (e.g. for vectorisation, for deletion of NA's, etc.) So ideally I would like to store them in some way that their names only have to be typed once and can always

Hi all, Dear Dr. Wheeler, I am trying to use the lmPerm package to perform multiple regression on microarray data with certain empirical variables associated with treatments of the experiment. In order the circumvent the very conservative multiple test corrections such as Bonferroni and BH, I try to use permutated probabilities to assess associations. In addition to mu previous posting I

Dear R-helpers, I would like to do a Spearman rank order test, and used the cor() function with the method "spearman". It gives me a number (correlation coefficient?) , but how can I get the p-value? Thank you for the help in advance! Regards, Anne-Katrin -- [[alternative HTML version deleted]]

Dear Madams and Sirs, because the DOI foundation recommends a new, secure resolver [1], I wanted to suggest the attached patch. It a) updates a static DOI link in the docu, b) the code chunks that generate new ones, and c) a reg-ex that detects DOI links. Hopefully, my first venture into SVN & patch files was done correctly. My apologies, if there is anything amiss. In that case, I'll

Hello, I would like to label outliers (or all dots) in a plot plot(as.matrix(ValAddInd_byYear), as.matrix(Cons_Elec_Ind_byYear)) ideally by the 30 different row.names, but anything else (x or y values, for example) would already be helpful. I've tried various things with pch, but they didn't work and I can't find anything else. Any hint appreciated! Thanks, Katrin

Hello, I've downloaded the package "mass" and I've tried to install it, but it didn't work. Can you help me (I use Windows)? Katrin -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20020707/e560f42b/attachment.html

Dear all, how can I plot a line graph and a bar graph in one single figure? I tried to combine "barplot" and "plot". Even though they both have the same x-values (1 to 55),  it just doesnt look as if they match in their scale (the barplot is much wider than the "plot"....even though I tried to put limits on the x-axis). Here is an example of what I did:

Dear all, a program that worked well for weeks today gave me consistently the error message no lines available in input referring to the lines for (i in (0:(timeintervals-1))) { j=subjectquantity+6+i*(subjectquantity+7) print (j) results<-read.table(file, header=F, skip=j, nrows=subjectquantity) timeintervals have been specified as argument as 12 subjectquantity as 9 the file it should

Dear All, I learned about R during my studies at Cologne University of Applied Science. Now I work at Allianz Dresdner Bauspar AG and I would like to install R here too. Is there any license issues that need to be taken in consideration, any fees to pay by the company? I know there are answers to this in the FAQs but I didn't really understand the legal language. Could you therefore just give

Dear R-helpers, I have a time series analysis problem in R: I want to analyse the output of my simulation model which is proportional cover of shrubs in a savanna plot for each of 500 successive years. I have run the model (which includes stochasticity, especially in the initial conditions) 17 times generating 17 time series of shrub cover. I am interested in a possible periodicity of shrub

Dear R Help-Team, If I start R and type in search() the result is: [1] ".GlobalEnv" "package:ctest" "Autoloads" "package:base" I often need an Add-on-package. I can load it with "> library(package)" or include a require-command in the .First-function. The .First-function is fine, if I always use the same working directory. But I

Hallo, I'm using the R Version 2.3.0 (2006-04-24) on Suse Linux 10.1. With an older R and Linux version I could write a R-function into a file and execute it with the command: R --vanilla < script.r In the file script.r was code like this: postscript(file="results.eps") x<-2 y<-3 plot(x,y) dev.off() It worked fine, but now nothing happens. In the manual page this option

Dear List, How do I pick the p-value out of the urppTest result? For adfTest the p-value can be extracted by A2 at test$p.value following A2 <- adfTest(myData[,i], lags=2, type=c("c")) What do I do for urppTest? The above doesn't seem to work. There is a slot @test with $output, which is a list of various test results that didn't want to give away only the p-value

Hi, I would like to know if there's a way to detect the workgroup of my own server with libsmbclient. With earlier versions than 3.0.20 it was possible to use smbctx->workgroup after the context was initialized. Since 3.0.15pre2 or so its always "WORKGROUP". Am I missing something? Regards, Gerd Fleischer ______________________________________________________________

Hello, I am using maptools for plooting geographical data. The colour of the region indicates some region dependent value (population for example). I pass the colours of the regions to the plot.Map function by defining the foreground colour: jet.colors = colorRampPalette(c("#00007F", "blue", "#007FFF", "cyan", "#7FFF7F", "yellow",

Does anyone know of an R-function or method to compare two C-statistics (Harrells's C - rcorr.cens) obtained from 2 different models in partially paired datasets (i.e. some similar and some different cases), with one continuous independent variable in each separate model? (in a survival analysis context)? I have noticed that the rcorrp.cens function can be used for paired data. Thanks

Dear, I want to analyze an outcome in an RCT using lme but I am not sure that I have chosen the right way for the model. We measured the outcome three times repeatedly in the same patient. One time before intervention and two times after intervention. I wanted to adjust for the correlated data in the repeated measurement and baseline differences in the variable in order to get the treatment