search for: poulsen

...boratory equipment (e.g. spectrophotometers) to have a linear response in a certain interval and then go into saturation. I wonder if there is an R-function that models this; for instance by estimating the breakpoint and fitting a line below the breakpoint and a constant above. Best regards Thomas Poulsen -- Thomas Poulsen Research Scientist. PhD, MSc Novozymes A/S Protein Design / Bioinformatics Brudelysvej 26, 1US.24 Phone: +45 44 42 27 23 DK-2880 Bagsv??rd.

...ments and an other file with 10 rows of one value each. When I read in those files using read.delim(), I get data frames, and apparently I can no longer do the multi-linear regression. Is there a way to convert the data frames into arrays, or am I going the wrong way about this? Sincerely Thomas Poulsen

The Register is carrying a article written by Kevin Poulsen of Securtiy Focus, calling asterisk "..the most powerful tool for manipulating and accessing CPN data.." > http://www.theregister.co.uk/2004/07/07/hackers_gut_voip/ I hope NuFone doesn't drop asterisk-set-able callerid's after this article; i've been wanting that feature...

...P/66.207.32.30-40700600' This indicates that I am passing the wrong context to Directory, however I am able to do mailbox lookups. It is just when i try to say 'yes this is the mailbox I want' that I have a problem. Any pointers are greatly appreciated. Thank you for your time. Dan Poulsen

Hello R-Users, I am trying to label points on an xYplot. I can label the points, but then I lose the error bars. I suspect that I should be using panel.xYplot rather than panel.xyplot, but that didn't work either... Can anybody help me out? Example code is below. First xYplot has error bars, second tries to label the points and the bars disappear. Thanks, John a=c(1:4)

...l pllan so I know I have a working system. 2. Once I am in the ACMEInc mailbox, if I dial 9 to reach the directory, it doesn't do the transfer. How can I set some kind of global paramater that says if I dial 9, go to Directory(${EXTEN}) ? Your input is greatly appreciated. Thank you. Dan Poulsen

Hello, I am trying to run a logistic regression with random effects on proportional data in glmmBUGS. I am a newcomer to this package, and wondered if anyone could help me specify the model correctly. I am trying to specify the response variable, /yseed/, as # of successes out of total observations... but I suspect that given the error below, that is not correct. Also, Newsect should be a

...roblem surfaced runs samba-3.0.28a-0.fc8 and I think yum updated this quite recently. It it likely that the behavior changed an a recent Samba update? Is there a good source of documentation that I just plain overlooked? Should I be using a different mechanism to set up the access rights? / Lars Poulsen

Dear list, I would like to call a Kerberos-authenticated web-service from within R. Curl can do it: $ curl --negotiate -u : "http://my.web.service/" so I would expect that RCurl also has the capability, but I have not been able to find the correct options to set. listCurlOptions() does not return anything with negotiate, and searching the source of RCurl, the only thing I found was

Hello R-users, I have code that gives me the important variables from an analysis. I need to input these variables into a different analysis. To do this, I need to modify them slightly... 1) remove all numbers at the end of the variables, 2) remove all periods. I tried to do it with the awkward code below. It works to remove all the numbers, but when I try to remove the period everything

Hello, I want to build a function to call up a column of a data.frame by the names of the columns. I have column names that are sequentially named (col1, col2, etc.). How do I change a character expression into something that will be understood as a data.frame column. For example: example<-data.frame(cbind(col1=1:10, col2=21:30, col3=41:50)) call.fun<-function(t){

Hello, I am trying to plot a curve over points plotted with se's in xYplot (see example below). I can get Figure 1 below to plot the data with error. However, I keep getting a the error message "Error using packet 1 object "y" not found" Can anyone see what I am doing wrong? Thanks! John ------------------------------- a=0.002; b=31.7; c=0.51

Hello, I am trying to make a xyplot plot with points that are different symbols. I want to call the symbol type (pch) from a column in my dataframe. Here is a simplified example. In my real example I also have groups, which I have not included here. This example doesn't change the symbols or colors. Any help you can provide would be appreciated. Thanks, John x<-c(1:12)

Hello, I have been using mgcv to run GAM hurdle models, analyzing presence/absence data with GAM logistic regressions, and then analyzing the data conditional on presence (e.g. without samples with no zeros) with GAMs with a negative binomial distribution. It occurs to me that using the negative binomial distribution on data with no zeros is not right, as the negative binomial allows zeros.

Hello, I apologize in advance as this example seems really elementary. Below I have created a simple scatterplot with lines. I would like to label each line with the name of the village, instead of using a legend. Can someone please let me know how to do it. Thanks in advance -- John Pop=c(1406,1862,2333,2664,2685,2521,2726,0,0,0,0,214,289,365,900,952,1002,1205,1041,942,1004,

Dear list, using the packages Hmisc and lattice i produced some nice xYplots. However, since the data range of the conditioning variable is very big, i need to define more than one y-scale for the plot (in some panels you just see a flat line of data points very close to the x-axis), e.g. different y-axis scales for the different rows of the plot. Is there a way to do so? Thanks, Henning

Hello, The simple example below plots three species by year and site. However, I would like the points to occur in groups horizontally, rather than on top of each other. I thought jitter might work, but haven't been successful in getting it to do what I want. Any ideas? Thanks in advance. John wt=rnorm(27,5,0.6) yr=rep(c(2000,2001,2002),each=3)

Hello, I am using glmer() from lmer(lme4) to run generalized linear mixed models. However, I am having a problem extracting the standard errors for the fixed effects. I have used: summary(model)$coef fixed.effects(model) coef(model) to get out the parameter estimates, but do not seem able to extract the se's. Anybody have a solution? Thanks, John

Hello, I have a dataset that has counts, but I need to expand the dataset so that each of the counts has its own line in the dataset (row) and is given and id. It looks something like: Site Type Cnt 1 "A" 3 1 "B" 0 2 "C" 2 I want the dataset to look like: Site Type ID 1 "A" 1 1 "A" 2 1 "A" 3 1 "B" 0 2 "C" 1 2

Hello, I am looking for some general advice on which graphics package to use to make a figure demonstrating my experimental design. I want to design a circle with 7 sectors inside. Then I will want to shade the sectors depending on densities of observations in the sectors. I will also want to draw horizontal lines at increments along the sectors to demonstrate different distances out to