search for: 0.825

I have a file: x y z 0.025 0.025 1.65775 0.025 0.050 1.62602 0.025 0.075 1.63683 0.025 0.100 1.91847 0.025 0.125 2.00913 0.025 0.150 1.82222 0.025 0.175 1.70901 0.025 0.200 1.39759 0.025 0.225 1.39089 0.025 0.250 1.04762 If I read the file like this: data<-read.table("file.dat") How do I access the vectors x,y,z that are inside the dataframe data? I studied Venables and

hi, I have 2 files example 1 and example 2 and would like to know what is in example2 and not in example1 (attached) V1 contain data which could be in duplicated which I am using as identifiers I used setdiff(example2$V1,example1$V1) to find the identifiers which are specific to example2: [1] "rs2276598" "rs17253672" I am looking for a way to get an output with all

Hi All, simple question: do you know how to graph the following object/matrix in a 'surface manner': [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.154 -0.065 0.129 0.637 0.780 0.221 [2,] 0.236 0.580 0.448 0.729 0.859 0.475 [3,] 0.401 0.506 0.310 0.650 0.822 0.448 [4,] 0.548 0.625 0.883 0.825 0.945 0.637 [5,] 0.544 0.746 0.823 0.877 0.861 0.642 [6,] 0.262 0.399

Hi, I hope somebody can help me with this issue: I am doing population pyramids using the barplot command, so in the left side I have male age structure and in the right side the female age structure. To plot the male age structure I put the data in negative numbers. Now, I want to change the sign in the bar plot in such way that I have no-sign numbers, both in left and right side of the graph. I

Hello, I’m trying to follow the syntax of a script from a journal website. In order to create a regression formula used later in the script, the regression matrix must have column names “X1”, “X2”, etc. I have tried to assign these column names to my matrix ScoutRSM.mat using a for loop, but I don’t know how to interpret the error message. Suggestions? Thanks, Paul

Hi all, what indicates the parameter theta in the summary of a random effect panel model estimated with the plm function? example: data("Produc", package = "plm") zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, model="random", data = Produc, index = c("state","year")) summary(zz) Effects: var std.dev

Kevin Wright wrote: > library(nlme) > m2 <- gnls(conc ~ t1*(1-t2*exp(-k*time)), > data = df.Chloride, > start = list( > t1 = 35, > t2 = 0.91, > k = 0.22)) So my error was to use nls instead that gnls. Thanks a lot, Kevin. > summary(m2) > plot(m2) > lag.plot(resid(m2), do.lines=FALSE) >

Dear all, I?m having some problems getting optim with method="Nelder-Mead" to work properly. It seems like there is no way of controlling the step size, and the step size seems to depend on the *difference* between the initial values, which makes no sense. Example: f=function(xy, mu1, mu2) { print(xy) dnorm(xy[1]-mu1)*dnorm(xy[2]-mu2) } f1=function(xy) -f(xy, 0,

Dear list members, First of all thank you for your helpful advices. After your answeres to my firt mail I studied a lot (R-News n?5) and I tried to perform my analysis: First, to fit a GLM with a nested design I decided to use the function "lmer" in package "lme4" as suggested by Spencer Graves and Filippo Piro. I remember to you that my data were: land use classes, 3 levels

Dear Thomas, The head of my dataset > head(wsuv) parcel sp time censo treatment species 1 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1 2 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1 3 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1 4 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1

Hi all, I have a package that is specific to a task I was repetitively using a few years ago. I now needed to run it again with new data. However I am told it was built with an older version or R and will not work. How can I tweak the package so it will run on 11.1? It was a one-off product and has not been maintained. Is there a way to "unpackage" it and repackage it to work? I

Hi, When I run the survfit function, I want to get the restricted mean value and the standard error also. I found out using the "print" function to do so, as shown below, print(km.fit,print.rmean=TRUE) Call: survfit(formula = Surv(diff, status) ~ 1, type = "kaplan-meier") records n.max n.start events *rmean *se(rmean) median 200.000

Dear All, What is the preferred way to get Bayesian analysis results (such as those from MCMCpacki, MCMCglmm, and DPpackage) into LaTeX table automatically? I have been using the "apsrtable" package and similar functions in "memisc" package, but neither seems to handle MCMC output directly. Many thanks. Shige

Hello, I'm trying to adjust a non linear model in which the biological response variable (ratio of germinated fungus spores) is dependent on 2 covariates (temperature and time). The response to temperature is modeled by a kind of beta function with 2 parameters (optimal and maximum temperatures) and the time function is a 2-parameter Weibull. Adjustments with nls or gnls work, but I need to

Dear list members, First of all thank you for your helpful advices. After your answeres to my firt mail I studied a lot (R-News n?5) and I tried to perform my analysis: First, to fit a GLM with a nested design I decided to use the function "lmer" in package "lme4" as suggested by Spencer Graves and Filippo Piro. I remember you that my data were: land use classes, 3 levels

Howdy, I have done many searches and can't seem to find a way around this. I am reading in a .csv file where each row is a dataset and each column represents a position. The values are sparse (there are 2003 positions but usually only 100-200 values) and the idea is to be able to plot each dataset (row) and overlay them in different combinations. What I would like to do is have a plot

Dear list, I am trying to replicate some Stata results but having a tough time doing it in R. The goal is to obtain a difference-in-difference estimate in a model with simple state fixed effects. The "state" variable is a factor, but some levels are missing. It appears that Stata automatically recognizes this and works around it. It also automatically pick "ME" as the base

Hi all, First I dont have much experience with R so be gentle. OK, I am dealing with a dataset (~ tens of thousand lines, each line ~ 10 columns of data). I have to create some subset of this data based on some certain conditions (for example, same first column with another dataset etc...). Here is how I did it: # import data dat <- read.table( "test.txt", header=TRUE,

The time appears to be perfect inside dom0, however all the domU''s tend to have a slightly faster date which gets further out of sync every day. I''m currently using Xen 3.0.3 with Gentoo Linux, under 3.0.2 I had no problems with domU clocks. Are there any known issues which could cause this? I''d strongly prefer not to run ntpd in every domU, having all domU clocks in

I am trying to use the SeasonalMannKendall function in the Kendall package. My dataset (alb_data) is in the same format as the example dataset (manaus) in the package. > class(manaus) [1] "ts" > is.ts(manaus) [1] TRUE > typeof(manaus) [1] "double" > alb_data=read.table("R:/albemarle_manken.txt", header=T) > head(alb_data) year Jan Feb